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ABSTRACT 


It  is  proposed  that  flat  plates  can  have  their  buckling  strength 
increased  by  prestressing.  The  prestressing  is  accomplished  by  first 
cold  rolling  the  plates  into  cylindrical  form  and  then  opening  the 
plate  by  riveting  onto  or  clanping  into  a flat  frame.  This  pre stressing 
will  produce  membrane  stresses  in  the  plane  of  the  plate  of  such  a 
magnitude  and  distribution  as  to  raise  its  buckling  load  in  the  direction 
of  the  generator  of  the  cylinder.  The  present  investigation  was  restricted 
to  the  case  of  all  four  sides  clamped. 

The  stresses  produced  by  this  process  are  measured  and  their  effect 
calculated  by  the  Rayleigh-Ritz  method.  The  analysis  shows  that  small 
stresses  (1000  psi  maximum)  could  raise  the  buckling  load  almost  50$, 
and  probably  more.  These  stresses  would  also  change  the  buckling  mode 
in  some  of  the  cases  examined  from  antisymraetrical  to  symmetrical  modes. 

A testing  frame  was  constructed  to  produce,  and  measure  the  effect 
of  the  proposed  prestressing.  The  tests  showed  that  the  buckling  load 
was  raised  in  some  cases  over  100$  while  the  average  for  all  tests  was 
38$. 


A theoretical  evaluation  of  the  membrane  stresses  and  their  effect 
was  carried  out  yhich  involved  the  solution  of  the  non-linear  plate 
equations  of  Von  Karman.  Because  the  necessity  of  using  the  rather 
approximate  assumptions,  the  analysis  could  be  subjected  to  questioning. 
Nevertheless,  the  analysis  showed  that  the  buckling  load  at  one  particu- 
lar type  of  prestressed  plate  would  be  raised  approximately  15$. 
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SECTION  I 


Introduction 

Before  the  advent  of  modern  high  speed  aircraft,  wrinkling  of  the 
covering  of  aircraft  components  due  to  local  buckling  under  flight  loads 
did  not  introduce  any  serious  problems  because  such  wrinkling  does  not 
have  a material  effect  on  the  performance  of  the  aircraft,  Present  day 
high  speed  aircraft,  however,  often  operate  under  conditions  in  which 
the  local  airstream  Mach  number  adjacent  to  the  aerodynamic  surfaces  is 
nearly  unity;  consequently  even  slight  surface  irregularities,  such  as 
the  appearance  of  wrinkles  due  to  buckling,  may  cause  the  local  velocity 
to  become  supersonic.  Such  changes  from  subsonic  to  supersonic  flow  are 
accompanied  by  the  formation  of  shock  waves  in  the  flow  which  materially 
increase  the  drag  of  the  aircraft.  It  is  therefore  desirable  for  these 
aircraft  to  employ  a covering  which  will  not  buckle  at  flight  loads  and 
yet  will  not  introduce  any  undue  increase  in  the  weight  of  the  structure. 

One  type  of  such  structure  is  sandwich  construction.  Sandwich 
construction  consists  of  two  thin  external  or  face  layers  of  high-strength 
material,  such  as  aluminum  alloy  sheet,  bonded  to  a thick  internal  layer 
or  core  of  light-weight  material,  such  as  balsa  wood  or  cellular  cellulose 
acetate.  The  core  serves  to  separate  the  strong  faces  a fixed  distance 
apart,  thus  giving  the  structure  a high  bending  and  therefore  buckling 
strength,  but  without  a substantially  increase  in  weight.  This  type  of 
construction  has  been  the  subject  of  intensive  investigation  and  considerable 
work  has  been  done  at  the  Guggenheim  School  of  Aeronautics  of  New  York 
University  (1,2). 

In  this  report  an  attempt  is  made  to  introduce  another  type  of  construc- 
tion where  buckling  is  delayed  by  a method  of  prestressing.  The  method 
is  as  follows.  First,  the  thin  sheet  is  curved  by  cold  rolling  in  the 
direction  perpendicular  to  the  direction  of  compression,  then  it  is  opened 
elastically  and  attached  to  the  stiffeners  or  the  frame  of  the  panel.  By 
such  a process,  it  was  found  that  in-plane  stresses  are  induced  and  the 
buckling  loads  are  increased  materially.  Both  experimental  and  theoretical 
investigations  have  been  carried  out  to  study  this  method  of  prestressing. 

The  fact  that  the  buckling  loads  can  be  appreciably  raised  by  first 
overcurving  the  panel  and  then  elastically  bring  it  to  the  desired  curvature 
was  first  noticed  by  Welter  (3)  in  his  study  on  the  effect  of  imperfection 
on  the  buckling  loads  of  curved  plates.  Welter  gave  no  explanation  on  this 
phenomenom.  Later,  Cox  (4)  indicated  that  favorable  deformations  occurring 
from  the  elastic  bending  accounted  for  the  raising  of  the  buckling  load. 

This  view  was  substantiated  by  Cicala  (5)  who  showed  that  in  certain  cases 
the  buckling  load  could  be  increased  due  to  favorable  deformations. 

Although  favorable  deformations  will  doubtlessly  be  a factor  in  raising 
the  buckling  loads,  an  equally  important  cause,  if  not  the  major  cause,  is 
believed  to  be  the  presence  of  the  induced  in-plane  stresses.  This  is  born 
out  in  the  present  Investigation  where  curved  sheets  were  elastically  opened 
to  flat  ones,  because  in  general  no  favorable  deformations  can  exist  for  a 
flat  plate.  The  term  nin  general*1  is  used  in  the  above  statement  in  deference 
to  Cox’s  (4)  statement  that  a favorable  imperfection  could  be  possible  so  as 
to  develop  a mode  corresponding  to  a sliAtlv  higher  buckling  load. 
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SECTION  n 


OUTLINE  OF  THE  EXPERIMENTAL  PROGRAM 

Before  describing  the  experiments  carried  out,  Welter1  s (3) 
experiments  should  be  reviewed  at  this  point  since  it  was  his 
experiments  which  brought  to  light  this  particular  method  of 
prestressing  thin  plates  to  increase  their  buckling  load.  Welter 
presented  the  results  of  some  experiments  he  performed  with  no 
attempt  to  explain  the  phenomenom  discovered. 

He  performed  his  tests  on  both  17ST  aluminum  alloy  of  ,036  in, 
thickness  and  2S  3/4  H pure  aluminum  sheets  of  ,032  in,  thickness, 
both  of  24  in,  width  and  96  in,  length.  The  cylindrical  plates 
were  tested  with  a radius  of  curvature  of  24  inches.  The  Initial 
radii  of  curvatures  of  the  plates  varied  between  6 in,  and  24  in, 

(the  latter  being  a non-pre a tre s sed  plate)  and  were  tested  with  a 
radius  of  24  in,  by  opening  the  plate  elastically  to  the  greater 
radius.  These  tests  showed  Increases  in  the  buckling  loads  up  to 
50$,  There  were,  however,  a total  of  only  nineteen  tests  performed 
and  of  these  perhaps  only  half  the  results  were  conclusive. 

The  present  experimental  program  can  be  divided  into  three 
parts:  first,  to  determine  whether  objectionable  deflections  were 
produced  during  the  prestressing  process;  secondly,  to  determine 
whether  the  prestressing  process  produced  inplane  (membrane)  stresses 
that  could  be  capable  of  raising  the  buckling  load  of  the  flat 
plate;  and  thirdly,  to  measure  the  buckling  loads  of  prestressed 
plates  and  compere  them  with  that  of  the  corresponding  non-pre  stressed 
flat  plate. 


SECTION  III 

TESTS  TO  DETERMINE  THE  INITIAL  DEFINITIONS  DUE  TO  PRESTRESSING 

Since  the  purpose  of  prestressing  is  to  delay  the  occurrence  of 
wrinkles  due  to  buckling,  it  is  important  to  determine  first  whether 
it  is  possible  to  clamp  an  initially  curved  plate  without  undesirable 
initial  deflections.  By  "undesirable",  it  is  inferred  that  the  maximum 
deflection  of  the  plate  should  be  less  than  the  thickness  of  the  plate. 

Three  "flatness"  tests ’were  performed.  One  of  the  tests  was 
performed  on  a half -hard  brass  plate  while  the  other  two  were  performed 
on  24ST  aluminum  plates.  The  brass  plate  was  19  l/2  by  12  inches  and 
had  a nominal  thickness  of  l/32  inch.  The  aluminum  plates  were  19,25 
by  13,633  inches  and  had  nominal  thicknesses  of  ,032  and  ,040  inches. 
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The  clamping  of  the  plate  was  executed  on  the  compression  base 
of  the  200,000  lb.  Baldwin-Southvark  thiversal  Testing  Machine.  In 
this  base  there  are  several  conveniently  located,  "Tee " Slots,  which 
were* used  in  clamping  the  plates  to  the  base. 

The  clamping  of  the  brass  plates  (see  Fig.  1)  was  accomplished 
by  using  eight— l/2  inch  thick,  1 1/4  inch  wide  steel  bars.  Two  of 
these  bars,  19  l/2  inches  long  were  placed  parallel  to  each  other 
so  that  the  edges  of  the  bars  and  the  "Tee"  slots  coincided.  Two 
bars  9 l/2  inches  long  and  parallel  to  each  other  were  positioned  on 
the  base  so  that  they  were  perpendicular  to  the  19  1/2  inch  bars. 

The  brass  test  plate  was  set  on 
this  frame,  and  corresponding  steel 
bars  were  oriented  in  similar  positions 
on  top  of  the  brass  plate,  thus 
forming  the  clamping  frame. 

Twelve,  1/2  inch  standard 
steel  bolts  whose  heads  were  fitted 
with  lock-washers  were  registered 
into  the  "Tee"  slots  of  the 
compression  base  and  were  positioned 
as  indicated  in  Fig.  1.  The  top 
of  each  bolt  was  fitted  with  a large 
diameter  and  a standard  diameter 
steel  washer.  The  diameter  of  the 
large  steel  washer  was  such  that 
it  extended  over  a considerable 
portion  of  the  width  of  the  steel 
clamping  bars.  On  the  side  of  each 
bolt  opposite  to  that  of  the  steel 
bar  were  placed  three  shims,  two  of 
them  being  l/2  inch  thick,  and  one 
being  l/32  inch  thick.  Thus  the 
total  thickness  of  the  shims  equaled 
that  of  the  steel  bars  and  test  plate. 
By  fitting  a nut  on  the  top  of  each 
bolt  so  that  it  bears  on  the  washer, 
and  tightening  the  nut,  a clamping  of 
Fig.  1 the  spedimen  between  the  steel  bars 

was  obtained. 

After  the  plate  was  clamped  by  tightening  down  on  each  of  the 
twelve  bolts,  the  deflections  of  the  plate  at  several  positions  were 
measured  by  a Starrett  dial  gage  graduated  with  divisions  of  *001 
inches. 


Clamping  of  Brass  Plates  for 
"Flatness"  tests 
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The  damping  of  the  Aluminum  plates  (see  Fig*  2)  was  accomplished 
by  using  eight-  5/8  inch  thick,  1 1/2  inch  wide  steel  bars  in  a manner 
similar  to  the  procedure  used  for  the  brass  plate*  In  this  series  of 
tests  two  holes  were  drilled  in  each  of  the  steel  bars  and  half  the 
bars  were  tapped  and  the  other  half  drilled  clear  to  accomodate  bolts 
that  would  clamp  the  steel  bars  tightly  to  the  edges  of  the  plates. 

Two  pairs  of  the  steel  bars  were  19*25  inches  long  while  the  other 
two  pairs  were  10,75  inches  long.  The  ends  of  the  shorter  bars  which 
were  underneath  the  plate  were  beveled  to  leave  a passageway  for  the 
wires  connected  to  the  strain  gages  on  the  bottom  surface  of  the  plate 
as  required  in  the  determination  of  inplane  strains  described  in  the 
next  section* 

The  bolts  for 
clamping  the  plate  to 
the  base  of  the  testing 
machine  were  specially 
modified  to  fit  the 
"Tee"  slots  thus  making 
clamping  of  the  plate 
easier.  The  hexagonal 
head  of  the  bolt  was 
replaced  by  a square  head 
machined  to  fit  the 
bottom  of  the  "Tee"  slot 
exactly.  Instead  of 
large  washers  being  used 
to  transmit  the  clamping 
force  from  the  bolt  to 
the  steel  bars,  special 

Clamping  of  Aluminum  Plates  for  "Flatness"  "clamps"  were  made  from 
Tests  3/8  x 1 1/4  x 2 inch  bars 

with  holes  drilled  through 
2 the  center  to  accomodate 

the  bolt.  One  end  of 
these  "clamps"  rested  on 

the  steel  bars,  forming  a frame  in  which  the  plate  was  held,  while  the 
other  end  rested  on  bars  1*30  inches  high  (approximately  the  height  of 
the  frame). 

The  dial  gage  was  fastened  by  means  of  rods  to  a gage  "base  plate". 
The  gage  was  adjusted  so  that  its  measuring  rod  was  perpendicular  to 
the  plane  of  the  test  plate.  The  measurements  at  the  various  positions 
were  obtained  by  sliding  the  gage  "base  plate"  along  the  surface  of  the 
compression  base  of  the  testing  machine,  until  the  rod  of  the  gage  was 
situated  directly  on  the  desired  point.  At  each  point,  the  measuring  rod 
was  lifted  a small  distanoe  above  the  plate  and  then  released,  so  as  to 
obtain  the  same  measuring  pressure  at  each  point. 

The  various  measuring  positions  are  illustrated  in  Figs.  3,  4,  and 
5 for  the  brass  and  aluminum  plates*  A particular  point  is  defined  by 
two  numbers,  the  first  indicating  the  position  along  the  short  side  of 
the  frame,  the  second  along  the  long  side* 
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Deflection  Measuring  Points 
Fig.  3 


Deflection  and  Strain  Gage  Measuring 
Points  on  First  Series  of  Aluminum 
Plates 

Fig.  4 


Deflection  and  Strain  Gage  Measuring  Points 
on  Second  .040  inch  Aluminum  Plate 

Fig.  5 

The  reference  point,  or  the  point  at  which  the  dial  gage  was  set  to 
zero,  was  arbitrarily  selected  at  point  5-9*  for  the  brass  plate  and 
1-1  for  the  aluminum  plates. 
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The  results  of  the  tests  are  in  Tables  1 (a,b,  and  c). 

The  comercial  designation  of  the  brass  plate  was  1/32  inch  with 
a standard  commercial  tolerance  of  * .002  inches.  From  the  res  nits, 
it  can  be  observed  that  the  maxima  positive  deflection  is  ♦ ,003 
(position  4-1)  , and  the  maximum  negative  deflection  - .007  inches 
(position  3-3,  3-4,  3-6).  Thus  the  percentage  maximum  deflections 
based  upon  the  plate  nominal  thickness  of  l/32  inch  are  + 9.6$  and 
- 22.456  or  a total  of  3256. 

The  .040  inch  thick  aluminum  plate  had  an  Initial  radius  of 
curvature  of  8.7  inches,  and  when  flattened  had  a maximum  total 
deflection  (the  difference  between  the  highest  and  lowest  points  on 
the  plate)  of  .022  inches  or  55%  of  the  nominal  thickness. 

Similarly  the  .032  inch  thick  aluminum  plate  with  an  initial 
radius  of  curvature  of  7 3/4  inches  and  a maximum  total  deflection 
of  .049  Inches  or  15356  of  the  nominal  thickness.  The  unusually  large 
deflection  of  this  thin  plate  was  probably  caused  by  initial  imperfections 
(dim  to  the  fact  that  this  plate  had  been  flattened  several  times 
previous  to  this  experiment)  and  unavoidable  errors  in  the  measurements. 
However,  it  should  be  noted  that  deflections  of  this  order  might  be 
expected  in  aircraft  due  to  the  pressure  of  the  airforces. 

SECTION  IV 

TESTS  TO  DETERMINE  THE  INPLANE  STRESSES  DUE  TO  PRESTRESSING 

There  were  two  series  of  tests  performed.  In  the  first  series, 
three  plates  were  tested  - two  brass  plates  and  one  aluminum  plate. 

These  plates  were  tested  before  any  theoretical  work  was  carried  out. 

In  the  second  series,  tests  were  performed  on  four  plates.  These 
plates  were  of  the  same  size  but  of  different  thickness.  The  theory 
indicates  that  the  stress  distribution  depends  on  the  parameters 
2 A 

a and  a/b,  where  a is  the  width  of  the  plate,  b the  length,  h the 
Hh 

thickness,  and  R is  the  initial  radius  of  curvature.  The  radii 
of  curvature  of  these  plates  were  rolled  in  such  a way  that  the 
parameter  remains  a constant  value,  337.5. 

Each  of  the  seven  plates  was  clamped  flat  in  the  frames  described 
in  the  preceding  section.  SR-4  strain  gages  were  fastened  to  the 
plates  with  Duoo  Household  cement  in  pans  back  to  back  at  various 
locations.  The  leads  of  the  strain  gages  were  attached  to  a multichannel 
switching  and  balancing  unit.  The  strain  was  measured  by  a standard 
Baldwin-Southwark  Type  K strain  indicator. 
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The  strain  in  a particular  gage  la  obtained  by  subtracting  the 
Initial  reading  taken  when  the  plate  is  standing  free  on  its  curved 
side  and  the  final  reading  taken  when  the  plate  Is  clamped  In  the 
frame.  The  inplane  strain  is  obtained  by  taking  the  average  of  the 
strains  on  both  sides  of  the  plate  at  a particular  point. 

In  each  case  the  radius  of  curvature  vas  determined  by  tracing 
the  outline  of  the  plate  on  a piece  of  paper.  As  the  surface  Is 
approximately  cylindrical,  the  outline  is  approximately  circular.  Thus 
by  drawing  several  chords  of  the  arc  and  finding  the  intersection 
of  the  perpendicular  bisectors  of  these  chords,  the  radius  of  curvature 
is  determined. 

The  tests  on  the  two  brass  plates  did  not  produce  consistent 
results.  A review  of  the  literature  of  the  Baldwin-Soutbwark  Company 
showed  that  bonding  of  paper  backed  strain  gages  on  brass  is  usually 
poor.  This  probably  explains  the  inconsistent  readings.  However, 
the  results  of  these  tests  did  indicate  that  inplane  strains  are 
produced  in  such  a process  and  the  inplane  strains  in  the  direction 
of  the  buckling  are  negligible  compared  to  the  inplane  strains  in 
the  perpendicular  direction. 

The  aluminum  plate  tested  in  the  first  series  was  a 14  inch  by 
19  l/4  Inch  24  - ST  alloy  plate  with  a thickness  of  .020  inch.  Ten 
SR-4  type  A-l,  and  two  type  A-5  strain  gages  were  mounted  at  the 
position  shown  in  Fig.  6.  The  test  results  are  presented  in  Table  2. 

As  indicated,  the  first 
two  tests  were  performed  with 
gages  1 through  8 only.  Since 
the  inplane  strains  are  not 
entirely  as  expected,  it  was 
decided  to  investigate  the 
strain  distribution  over  the 
top  and  bottom  surfaces  as  each 
of  them  should  be  smooth  and 
continuous.  An  investigation 
of  the  data  of  Table  2 reveals 
that  on  the  top  surfaces 
(gages  1,3, 5, 7)  all  gages 
except  7 yields  a smooth 
distribution  of  strain.  A 
similar  conclusion  is  obtained 
for  gage  2 on  the  bottom 
surface.  It  was  therefore 
decided  to  mount  gage  9 
immediately  adjacent  to  gage 
7 and  gage  10  immediately 
adjacent  to  2.  Gages  7-2  and 
8-2  (see  Fig.  6)  were  attached 
for  cheeking  purposes. 


Location  of  Strain  Gages  on  .020 
Aluminum  Plate 

Fig.  6 
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From  the  results  of  tost  3,  it  Is  seen  that  by  employing  gages  9 
and  10  instead  of  7 and  2 a smooth  strain  distribution  is  obtained*  By 
extrapolating,  the  values  shown  in  Table  5 under  "Test  3 (using 
correction)"  are  obtained.  By  comparing  the  axial  strain  at  7-2,  8-2 
with  7,8,  the  rough  extrapolations  are  probably  somewhat  in  error.  On 
a conservative  basis  the  following  inplane  strains  were  assumed.  These 
values  will  be  used  in  calculating  the  buckling  load. 


Position 


Inplane  strains  (xlO  ) 


3,4 

5,6 

1,2 

7,8 


60 

80 

10 

15 


3h  the  second  series,  four  aluminum  plates  of  the  size  of  (and 
including)  the  aluminum  plates  described  in  the  preceding  section  were 
tested.  The  locations  of  the  strain  gages  for  three  plates  (namely 
one  each  of  .032,  .040,  and  .051  inch  thickness),  are  shown  in  Fig.  4, 
while  for  the  remaining  test  (on  a ,040  plate)  are  shown  in  Fig.  5. 


In  this  series  of  tests  the  zero  readings  were  taken  v ' the 
plates  standing  on  its  curved  (shorter)  edge,  and  with  and  without 
the  bars  attached  to  the  longer  edges.  The  reason  for  taking  zero 
readings  with  bars  attached  to  the  longer  edges  was  to  force  these 
edges  to  be  straight,  which  they  would  have  been  had  adequate  means  been 
available  for  cold  rolling  the  plates.  The  data  taken  for  the  five 
tests  are  shown  in  Tables  3 through  10.  Note  that  in  some  cases  the 
tests  were  repeated  on  the  plates  several  times. 


In  the  case  of  the  .040  aluminum  plate  whose  data  is  recorded 
in  Table  10  many  of  the  readings  should  be  identical  and  this  property 
of  symmetry  was  used  in  analyzing  the  data  to  eliminate  any  unreasonable 
data.  This  analysis  is  shown  in  Figs.  7a  and  b where  the  observed  data 
and  the  derived  averages  are  indicated. 


Graph  Used  for  Averaging  the  Values  of  the 
•040  inch  Aluminum  Plate 


Fig.  7a 
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Fig.  7b 


Using  Pig*  7a  and  b to  arrive  at  a "reasonable  "estimate  for  the 
membrane  strains  of  the  plate  in  the  x-direetion  for  the  upper  left 
(of  Fig.  5)  quadrant  we  obtain  the  values  listed  in  Table  11.  Similarly 
in  the  cases  of  the  other  plates  referring  to  Fig.  5 we  find  because 
of  symmetry;  the  membrane  strains  as  listed  in  Table  12  through  14* 


SECTION  V 

DETERMINATION  OF  THE  BUCKLING  LOAD  BT  FINITE  DIFFERENCE  APPROXIMATION 

Assuming  that  the  inplane  stresses  due  to  prestressing  do  not 
change  while  the  buckling  load  is  applied,  the  buckling  load  can  be 
calculated  by  means  of  finite  difference  approximation.  Two  objections 
can  be  raised  against  such  a procedure.  First,  since  there  are  small 
initial  deflections  due  to  prestressing,  the  theoretical  buckling  load, 
which  defines  the  stability  limit,  does  not  exist.  Secondly,  these 
inplane  stresses  due  to  prestressing,  because  of  the  initial  deflections, 
certainly  will  change  while  the  buckling  load  is  applied.  The  excuse 
for  using  such  a procedure  to  calculate  the  buckling  load  is  as  follows. 
Although  the  theoretical  buckling  load  is  not  defined  for  a plate  with 
initial  deflections,  the  load  for  such  a plate  at  which  the  deflections 
suddenly  become  large  is  close  to  the  buckling  load  for  the  corresponding 
plate  without  initial  deflections*  Although  the  change  in  the  inplane 
stresses  due  to  prestressing  while  the  buckling  load  is  applied  may 
introduce  large  errors  in  the  calculation,  in  the  absence  of  a better 
method,  it  is  felt  that  the  present  method  will  at  least  give  some 
qualitative  results  to  compare  with  the  experimental  results  described 
in  the  next  section. 


The  differential  equation  governing  the  deflection  of  a thin  plate 


with  inplane  forces  N , N and  N is 

x y xy 


N 2 N -2  N -2 

H ♦ -*  2-S  ♦ 2 -S  2ji 
D ax  D ay  D 3x3y 


(1) 


Let  the  plate  be  loaded  in  the  y-direction 


Nw„  + N N 
ycr  yO  xy 


V>  where  the  subscript 


inplane  forces  and  the  subscript  cr  refers 

Slues  Nx  - oh,  Hy  - oh, 

experiments  it  was  found  that  c ~ * a « ■ 

yO  xyO 

stress  we  find 


only.  Then  N^  = N^,  Ny  * 

0 refers  to  the  original 

to  the  buckling  value, 
and  from  the  preceding 
0,  from  Hooke's  Law  for  plane 
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X 


Eh 

1-v5 


'xO 


„ _ vEh 

I*  * N + X c^, 

7 ycr  ^ v2  10 


N * 0 
xy 

Substituting  into  equation  (1),  we  have 

2 


Eh" 


2-*  + IN + 


2 xO  . 
v ax 


VEh  \ a‘ 
ycr  1-v2  *° 


(2) 


It  is  convenient  to  put  equation  (2)  in  non-dimensional  form  by  letting 


2*  2 


w*  ■ v/h,  N » ■ N x»  « x/a 


ycr 


Eh 


ycr 


y’  - y/a,  ejo  - 


-xO 


We  have,  therefore, 

^4’  - 12  cLa  Sfsd  - (N  • + 12v  c*  ) 2-4  - 0 

^7  w xO  ax,2  ycr  yO  ^,2 


(3) 


where  is  the  biharmonic  operator  referring  to  the  non-dimensional 
coordinates* 


let  us  now  transform  this  equation  into  finite  difference  equations* 

Let  the  plate  be  divided  into  rectangular  meshes  with  ^ y » k*y,  where 

A x and  Ay  are  the  sides  of 
meshes.  For  the  19  1/4" 
x 14n  *020"  Aluminum  alloy- 
plate  clamped  by  the  1 1/4" 
bars,  the  plate  area  is 
actually  (19  3/4"  - 2 1/2") 
x (14"  - 2 2/2") • If  we 
divide  the  plate  as  shown 
in  Fig.  8,  we  have 

y - b/4,  k - - 1.457,  v 

11.5 


Net  Points  for  Finite  Difference  Solution 
Fig.  8 


■ 
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To  approximate  the  boundary  conditions  more  accurately,  imaginary  net 
points  outside  the  plate  are  introduced.  From  the  boundary  conditions, 
they  are  equal  to  the  values  one  mesh  point  inside  the  boundary  as 
indicated  in  Fig.  8.  The  finite  difference  equations  at  various  points 
are  as  follows:* 

At  vQ: 

(177.7+18.52  +2Ny£r)v0-2(44.4+5.66  c^))w1 

-2(94.2+3.6  c»q  +Ny£r)v2+4(l5.09)v3  - 0 (4) 

At  w^: 

(177.7+18.52  +2Ny*r)w1+2(3.56)w1 

-(44.4+5.66  c^1)w0+2(15.09)w2-2(94.2+3.6  + Ny'r)w3+0 


At  w2: 

(177.7+18.52  +^y£r)w2+2(l6*01)w2 

-(94.2+3.6  c^2  + Ny^r)w0-2 (44.4+5.66  ei2)w3«(X5.°9)w1  (6) 

At  v3: 

(177.7+18.52  c^3  +2Ny^Jw3+2(l6.01)w3+2(3.56)w3 
-(94.2+3.6  c^3  + Ny»r)w1+15.09w0-(44.4+5.66  c^)v3  = 0 


* The  transformation  of  partial  differential  equation  to  finite  difference 
equations  can  be  found,  for  example,  in  "Applied  Elasticity"  by  Chi-Teh 
Wang,  McGraw-Hill  Book  Co.,  1953,  pp.  106-143. 
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where  the  w’s  at  various  mesh  points  are  the  non-dimensional  deflections. 

Using  the  results  for  the  inplane  stresses  in  the  .020  inch  thick 
aluminum  plate  found  in  the  previous  section,  we  have,  using  the  notation 
of  this  section. 


e I * 20  micro-inches  per  inch 


a t b n n n n 

®X1  ° 


o i b 3 n n nr 

X2 


0x3  B 5 " " n M 


The  finite  difference  equations  can  he  put  in  matrix  notation  as  follows. 


where 


(&H®  VJM 


[548 

-315 

-332 

60l 

-196 

648 

30.2  - 

■368 

/-IO3 

30.2 

265  - 

■1231 

Ll5.1 

-112.2 

-72.7 

309] 

~2 

0 

-2  0“ 

0 

2 

0 —2 

-1 

0 

2 0 

L° 

-1 

0 2 

w'  * 


Equation  (8)  can  be  rewritten  in  the  following  form 


(8) 


(9) 
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which  is  the  tisual  form  in  which  a solution,  by  iteration,  for  the 
eigen  values  N^  will  converge  to  the  lowest  eigen  value  Ny^. 

The  inverse  of  A can  be  found  (18)  and  is 


3.396 

2.078 

5.069 

3.832 

1.274 

2.737 

2.370 

3.955 

1.502 

1.038 

6.507 

3.533 

(10  3) 

.6502 

— 

1.137 

2.145 

5.318 

(10) 


Iterating  Eq.  (9)  according  to  Ref*  (18)  the  non-dimensional  buckling 
parameter  NyCr  for  the  prestressed  plate  is  found  to  be 

Nycr  * ~ "*3  (U) 

while  for  the  ordinary  plate  (setting  all  cx  * 0) 


N * - - 47.55 
ycr 


(12) 


Thus  using  the  inplane  stresses  found  in  the  previous  section  a 
solution  by  finite  difference  approximations  shows  that  the  buckling 
load  can  be  increased  about  10956* 


SECTION  VI 

DETERMINATION  OF  THE  BUCKLING  LOAD  BY  ENERGY  METHOD  USING  HARMONIC 

ANALYSIS  OF  THE  TEST  DATA 

The  finite  difference  method  used  in  the  previous  section  loses 
accuracy  when  used  with  only  a few  net  points.  As  an  alternative,  the 
Rayleigh-Ritz  method  can  be  used  in  finding  the  buckling  loads  of  the 
prestressed  and  non-pre stressed  plates.  The  assumed  deflection  will 
be  represented  in  a Fourier  series  and  It  will  be  convenient  to  represent 
the  inplane  stresses  found  in  the  previous  sections  also  by  a Fourier 
Series. 

The  Raylelgh-Rits  method  will  be  shown  in  wore  detail  in  Part  II 
of  this  report  dealing  with  the  theoretical  analysis  and  only  the  essential 
points  of  the  method  will  be  described  in  this  section* 

Let  v b?  the  later*1  deflection  of  the  plate  due  to  buckling*  Is 
the  case  of  a plate,  with  clamped  edges,  the  bending  strain  energy  in 
the  plate:  when  buckling  occur?  is 
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U ■ - f(~  (^-M  + ^-5)  dx  dy  (13) 

2 '/  0x  0y 

A 

where  D Is  the  flexural  rigidity  of  the  plate.  The  work  done  by  the 
inplane  forces  is 

W « - ffU*>  + N (— ) dx  dy  (14) 

2 ''  x dx  7 dy 

where  N ■ N N =5?  N . In  the  above  we  hare  assumed  that  the 
x xu  y ycr 

contribution  to  the  work  by  the  original  stresses  in  the  y-direction 
and  by  the  original  shear  streams  Is  negligi  ble.  At  buckling  we  hare 


U - W « 0 


from  which  we  obtain 


N 


U - W 


ycr 


where 

and 


i rr  to 

^ 'J  dy 


dx  dy 


W - W - N_ 


ycr 


; fj  <-> 


dx 


dx  dy 


(15) 

(16) 


It  must  be  emphasized  that  the  above  relations  are  correct  only 
if  there  were  no  initial  deflections.  With  the  presence  of  the  initial 
deflections,  the  values  of  the  inplane  forces  due  to  the  prestressing 
will  probably  be  changed  when  the  load  N^cji  Is  applied,  and  therefore 

the  above  method  of  calculation  may  induce  an  error  which  is  not  small. 

The  Rayleigh-Ritz  method  can  be  carried  out  as  follows.  First, 
the  deflection  w is  assumed  in  the  form  of  a series  every  term  of  which 
satisfies  the  boundary  conditions,  but  with  undetermined  coefficients. 
Substitute  the  series  into  equation  (15)  and  carry  out  the  integration. 
Then  minimize  the  resulting  expression  with  respect  to  the  undetermined 
coefficients.  This  yields  a system  of  homogeneous  equations  involving 
the  buckling  load  A non-trivial  solution  is  obtained  if  we  equate 

the  undetermined  coefficients  of  these  parameters  to  zero,  from  which 

N is  determined. 

ycr 
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The  boundary  conditions  vill  be  satisfied  if  the  deflection,  v» 
is  assumed  in  the  following  form: 

v-  sin—  sin  2*  F Z sin  — sin  222  (17) 

a b nrt  n-1  a b 


As  in  Part  II,  let  us  Introduce  a stress  function,  F,  in  the  form 

oo  oo 


F ■ 2 2b  cos  222 

p*0  q=0  p,q  a 


cos 


SEE 


(18) 


where  b 


p»q 


K-  E h2, 
p»q 


N 


and  d 


xO 


xO 

h 


The  first  step  in  this  analysis  is  to  determine  the  coefficients 
of  the  series  (18)  corresponding  to  the  data  found  experimentally  and 
described  in  the  previous  sections. 


Since  it  was  experimentally  impossible  to  determine  the  strains 
at  the  edges  of  the  plate  we  vill  make  use  of  the  fact  that  along  cross 
section 


which  results  from  the  fact  that  it  is  assumed  that  no  traction  forces 
exist  at  the  boundaries.  The  coefficients  of  the  series  for  the  stress 
function  will  be  derived  by  the  method  of  collocation,  which  is  best 
described  by  actually  deriving  the  coefficients  for  a specific  set  of 
test  data.  The  problem  used  for  illustration  purposes  will  be  the  same 
one  that  va»  solved  by  the  finite  difference  method  in  the  previous 
section. 


The  points  at  which  the  strains  were  taken  will  be  numbered  as  is 
shown  in  Fig.  8 just  as  was  done  in  the  previous  section.  The  values 
of  the  strain  in  micro-inches  per  inch  will  be  rewritten  for  clarity. 

Gage  Position  Strain 

0 60 

1 80 

2 10 

3 15 
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Two  series  may  be  written  one  at  the  section  x * a/2  (containing 
points  0 and  2)  the  other  at  the  section  x * a/4.  Since  e^  dy  * 0 

then  the  series  will  contain  only  two  terms  each  and  be  of  the  form 


at  x * a/2  “ a2  cos  cos 


at  x » a/4  * p2  008  + P4  003 


(20) 

(21) 


Using  the  formulas  listed  by  Scarborough  (7)  in  his  section  on 
Harmonic  Analysis  and  making  use  of  the  fact  that  by  virtue  of  (19) 
<Xq  and  must  be  identically  zero.  The  values  of  the  coefficients 

are  found  to  be 


d2  - - 70 


o4-  -10 


- 15 


h - - 95 

The  required  two  dimensional  series  will  be  of  the  form 


(22) 


2jrx“7 

CO0  ~~  / 

^ a / 


(23) 


Then  substituting  (21)  and  (22)  into  (23)  the  following  conditions  are 
arrived  at  for  the  coefficients 


at  x * a/2 
at  x » a/4  a^2 

Solving,  we  obtain  . 


a0,2-*2,2  “ “ 70 


- 95 


■0,2  " - 95 

a2,2 

•0,4*- 15 

a2,4 

a0,4~a2,4  * ~ 10 
a0,4  " - 15 

® - 25 
* - 5 


Assuming  that  CyQ  0,  Hooke's  law  yields 

„ _ E 


"x0~  x_y2  Cx0 


(24) 
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then  using  Eq.  (18)  to  relate  the  stress  function  (and  coefficients) 
to  the  membrane  stresses  we  find 

, . e * - S 2 q2  b cos  cos  92H  (25) 

tr2(l-v2)  30  p»q  a b 


o'  2 *" 

-5-5 5“  » - 2 S q K cos 

hV(l-v2)  10  *>»« 


ES  eos  SE 


then 


£ -EiS  x 10' 

.2  2/.  2x  2 X AU 

h it  (1-v  ) q 


(the  factor  of  10”°  appears  as  a conversion  from  micro-inches  to  inches) 
then 


v * •8835 

K2,2  “ *2325 

*0,4  * *°349 

K0  . « .0116 


b0,2  * *8835  E h 

b2  2 - .2325  E h5 

b0,4  * *0349  E h< 

b . - .0116  E h2 


Assuming  the  deflection  w to  be  suitably  represented  by  the  series 


w * sin  sin2  ^ 1 v,  - sin  — + w-  . sin  ^ 
a b l lfl  a 3,1  a 


The  strain  energy  is  calculated  according  to  Eq.  (13)  and  found  to  be 


128a* 


^ 227.92w1>;l2  ♦ 892.59w3fl2  - 258.08*.^  ^ (29) 
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Calculating  T according  to  Eq.  (16)  results  in 


128 


| 17.484  vl9 


x2+  11.656  w3>12  - 11.656  wx 


and  using  the  non-dimensional  buckling  parameter  defined  as 


K * K 


A 


V "-^7  f 17,484  Vl.l2*  11,656  "3.12  - 11,656  W1,1W3,1 } 


(30) 

Thus  the  governing  equation  for  the  non-pre stressed  plate  are 
8(M  T)  d(U-N  T) 

— 1 — ■»- 1 1 * 0 and  ■ 0 (31) 


aw. 


av- 


er 


and 


'1,1  3A 

(455.84  + 34.968K)  v^  - (258.08  + II.656K)  - 0 

(258.08  + 11.656K)  vlfl  - (1785.18  4 23.312K)  v3>1  - 0 


Setting  the  determinant  of  the  coefficients  of  these  two  equations 
equal  to  zero,  the  conditions  for  a non- trivial  solution  gives 


K2  4 98.680K  4 1099.9  - 0 
from  which  we  find 

K - - 12.81  or  - 85.87 


or 


12.81 


(32) 


(33) 
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which  agrees  with  the  solution  given  by  Timoshenko  on  p.  323  of  Ref*  (19)* 
From  Eq.  (16)  Wz  can  be  calculated  and  making  use  of  the  values  for 

bp,„  <*)  la 

( W7*30  “I,!2  * “3.12  " M64*28  wl,lW3,x} 


(31) 


then  minimizing  the  expression  U-W  the  following  two  equations  are  arrived 
at 


3(P-Vx-KyT) 

^1,1 

8(P-Wx-NyT) 


(1310*44  - 34.968K)  - (1522.36  - 11.656K)  v3>1  - 0 

(35) 

- (1522.36  - 11.656K)  w1>;l  ♦ (5691.96  - 23.312K)  w3>1  » 0 


upon  setting  the  determinant  of  the  coefficients  equal  to  zero  the 
following  quadratic  equation  is  obtained: 

K2  + 285.73  ♦ 7568.6  * 0 (36) 


yielding  the  solution 


or 


K « - 29.545,  - 256.19 


N 


ycr 


- 29.545 


ir^D 


(37) 


Thus  on  the  basis  of  the  Rayleigh-Ritz  analysis  we  may  conclude  that 
the  stresses  occurring  in  the  .020  inch  thick  plate  because  of  the  pre- 
stressing  procedure  will  increase  the  buckling  load  approximately  130$* 
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A similar  set  of  calculations  can  be  performed  for  the  second  *040 
inch  thick  plate,  the  data  for  which  is  listed  in  Table  10*  In  this  case 
it  is  found  that  the  prestressed  plate  will  buckle  symmetrically  in  three 
half-waves  even  though  the  non-pre stressed  plate  buckles  anti-symme trie ally 
in  two  half-wave s. 

Assuming  the  deflection  v to  be  suitably  represented  by  the  series 
w » sin  — sin  sin  ( w,  0sin  — + v,  - sin  \ (38) 

a b b l lt2  a 3»2  a J 


the  buckling  load  for  the  non-prestressed  plate  is  found  to  be 


(39) 


which  is  within  0.6$  of  the  accepted  value  listed  in  Reference  (12). 

For  the  prestressed  plate  with  the  same  assumed  deflection  the 
buckling  load  is 


N_ 


ycr 


11.63 


tt^D 


If  however,  for  this  same  prestressed  plate  the  deflection  is  assumed 
to  be  suitably  represented  by  the  series 

V « sin  SS  ata  EE,ta  20.  f u r sta  SS  * w »ln  23 

a b b t 1,J  a 3,3  a 

the  buckling  load  for  the  prestressed  plate  is  found  to  be 


Nycr  * - 10.04  (41) 

Thus  in  the  case  of  the  .040  inch  plate  whose  data  is  described  in 
Table  10  the  plate  will  buckle  symmetrically  in  three  half-waves  at  a load 
19.8$  above  that  of  the  non-prestressed  plate. 
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Similarly  using  the  data  for  the  *032  Inch  thick  plate  whose  data 
is  listed  in  Table  12  it  is  found  that  this  plate  will  also  buckle 
symmetrically  into  three  half-waves* 


For  the  deflection  w assumed  as  the  series 


(40)  the  buckling  load  is: 
(42) 


while  if  assumed  as  the  anti-symmetric  series  (38)  it  is 


N 


ycr 


11.12 


Thus  In  the  case  of  the  *032  inch  plate  a symmetric  buckling  mode 
of  three  half-waves  occurrs  at  a load  13%  above  that  of  the  nonpros  tressed 
plate. 


SECTION  VII 

BUCKLING  LOADS  FOR  PLATES  WITH  INITIAL  IMPERFECTIONS 

Before  proceeding  to  describe  the  experiments  for  the  determination 
of  buckling  loads  for  the  non-prestressed  flat  plates  as  well  as  the 
prestressed  plates,  it  seems  desirable  to  review  the  existing  methods 
for  experimentally  determining  the  buckling  loads  for  plates  with  Initial 
imperfec tions » 

The  three  most  common  methods  for  determining  the  buckling  load  from 
experimental  methods  are  the  top  of  the  knee  method,  the  strain  reversal 
method,  and  by  Southwell’s  plot. 

In  the  "top  of  the  knee  method",  the  buckling  stress  i®  taken  as  the 
stress  corresponding  to  the  top  of  the  knee  of  a curve  of  load  versus 
the  lateral  deflection.  If  the  lateral  deflection  cannot  be  readily 
measured,  any  other  quantity  that  increases  in  substantially  the  same 
manner  as  the  lateral  deflections  may  be  plotted  instead.  Ctoe  such 
quantity  is,  the  difference  of  strains,  in  the  direction  of 
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loading,  on  the  two  side a of  the  "buckle  crest".  The  method  is 
illustrated  In  Fig.  9 and  the  buckling  loads  obtained  by  the  "top 
of  the  knee  method"  are  Indicated  on  the  curves  by  small  circles  for 
three  different  amounts  of  Initial  imperfections.  It  should  be  noted 
that  the  larger  the  initial  Imperfections  the  more  difficult  it  is  to 
judge  the  "top  of  the  knee"  and  consequently  the  buckling  load. 

The  "strain  reversal  method"  is  best  explained  with  the  use  of 
Fig.  10  which  shows  the  variation  of  strain,  c,  at  opposite  faces  of 
the  plate.  As  in  Fig.  9 the  method  is  Illustrated  for  three  different 
amount  of  initial  deflections  and  the  critical  stresses  found  by  this 
method  are  indicated  by  small  circles.  The  innermost  curve  represents 
a perfect  plate.  The  strains  are  elongations  (or  contractions)  which 
could  be  found  by  the  application  of  electric  strain  gages  on  opposite 
faces  of  the  plate  at  one  of  the  expected  crests  of  the  buckling  mode. 
Before  buckling  the  plate  is  under  a uniform  compression  and  both  gages 
are  under  the  same  strain  depending  only  upon  the  load  and  the  Modulus 
of  Elasticity.  At  the  buckling  load  the  initially  perfect  plate  suddenly 
deforms  thus  adding  to  the  compressive  strain  in  the  gage  applied  to 
the  concave  face  while  imposing  a tensile  strain  (thus  decreasing  the 
compressive  strain)  in  the  gage  applied  to  the  convex  face  of  the  plate. 
The  critical  load  is  based  on  the  above  behavior  of  the  perfect  plate 
and  is  consequently  defined  as  the  load  at  which  the  extreme  fibre 
strain,  e^,  on  the  convex  side  of  the  buckle  crest  stops  increasing 
and  begins  to  decrease.  It  can  be  seen  from  Fig,  10  that  again  the 
larger  the  initial  imperfections  the  more  difficult  it  is  to  choose  the 
critical  load. 


Mr 


Top  of  Knee  Method 


Strain  Reversal  Method 


Fig.  9 


Fig.  10 
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Southwell* 8 method  has  found  universal  acceptance  in  determining 
the  buckling  load  of  colims  but  haa  been  unsatisfactory  in  the  case 
of  plates  whose  lateral  motion  is  restrained  at  all  four  edges.  In 
brief,  Southwell's  method  depends  on  the  linear  relationship  between 
the  deflection  of  the  structure  and  the  ratio  between  the  deflection 
and  load,  d/P*  The  slope  of  the  straight  line  plotted  with  these 
coordinates  in  the  theoretical  buckling  load.  Southwell  (8)  showed 
that  there  is  a Identifiable  straight  line  for  slender  columns  while 
Donnell  (9)  and  Hu(lO)  have  shown  that  the  linear  relationship  will 
not  exist  as  long  as  there  are  extensions!,  strains  present  during  the 
buckling  process.  Donnell  (9)  has  further  shown  that  these  extensions!, 
strains  do  exist  in  the  case  of  a simply  supported  plate  whenever  the 
sum  of  the  initial  imperfection  and  the  deflections  due  to  the  load  P 
are  of  the  order  of  the  plate  thickness,  as  is  usually  the  case* 

Because  of  the  relatively  large  initial  deformations  present  in 
the  pre stressed  plates  the  three  commonly  used  methods  for  experimentally 
determining  the  buckling  load  mentioned  above  are  found  to  be  very 
difficult  to  use.  For  this  reason  a method  developed  by  Toshiki  (ll) 
is  tried  in  the  determination  of  the  buckling  load  of  the  plates  in 
this  investigation* 

Toshiki* s method  is  based  upon  the  fact  that  after  buckling  the 
deflection  of  a theoretically  perfect  plate  is  proportional  to  the 
square  of  the  load  and  therefore  the  load  versus  deflections  curve  is 
a parabola.  This  parabola  is  shown  in  Fig*  9 for  a plate  with  no  Initial 
imperfections  and  it  can  be  seen  from  the  figure,  that  for  plates  with 
initial  imperfections,  these  curves  are  assymptotlc  at  large  deflections* 

When  the  load  is  plotted  against  the  square  of  the  deflection,  or 
any  other  quantity  varying  in  substantially  the  same  manner,  the  post 

buckling  curve  will  be  as  straight 
line  intersecting  the  vertical  axis 
of  the  buckling  load.  In  Fig.  11 
we  have  plotted  the  same  lines  as 
were  plotted  in  Fig.  9 but  now 
(cx  - c2) 2 was  used  as  the  abscissa 

instead  of  (e^  - e2).  A theoretical 

analysis  demonstrating  the  parabolic 
post  buckling  behavior  of  clamped 
plates  is  carried  out  in  Appendix  A* 


O A 8 (.2  16  20 

(Volue  proportional  to  bending  strain) 


Toshiki* s Method 
Fig*  11 
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SECTION  VIII 


TESTS  TO  DETERMINE  THE  BUCKLING  LOADS 

Having  shown  the  probability  of  raising  the  buckling  load  by 
prestressing,  and  also  that  the  prestressing  caused  no  objectionable 
Initial  deflections,  a testing  program  was  Initiated  for  the  purpose 
of  actually  measuring  the  buckling  load  of  prestressed  and  non-pre stressed 
flat  plates.  It  was  recognized  that  the  desired  buckling  frame  should 
have  the  following  properties: 

1,  The  deflections  at  the  edge  of  the  plate  should  be  zero, 

2,  The  slope  of  the  plate  at  the  edges  should  be  zero, 

3,  The  stringers  (side  supports)  should  be  flexible  in  order 
that  they  transmit  no  load, 

4,  The  load  should  be  applied  uniformly  to  the  ends  of  the 
plate. 

The  stringers  designed  for  the  testing  program  are  shown  in 
Figs,  12  and  13,  The  plate  is  clasped  between  the  cantilevered  "fingers" 
(shown  in  Fig,  13)  which  are  formed  at  the  end  of  the  t ' ' vide  roots. 
Since  the  width  to  depth  ratio  of  these  thin  roots  is  about*  ;?OsI  the 
rigidity  ratio  of  the  cantilever  beams  is  about  400:1,  As  a cone  : >-smee 
of  this  high  "rigidity  ratio"  the  stringers  are  very  flexible  in  the 
direction  of  the  applied  load  and  yet  possesses  a great  deal  of  lateral 
and  rotational  rigidity. 


General  View  of  Buckling  Jig  Detailed  View  of  Stringers 

Fig.  12  Fig.  13 
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I The  flexibility  of  the  cantilevered  "fingers"  reduces  to  a 
minimum  the  applied  load  being  transmitted  by  the  stringers  while  their 
lateral  rigidity  will  produce  almost  zero  displacement.  To  produce 
the  condition  that  the  edges  of  the  plate  will  have  nearly  zero  slopes, 
a four  bar  linkage  was  designed  which  Is  shown  in  Figs,  12  and  14* 

The  two  stringers  are  attached  to  the  rod  end  bearings  at  the  ends  of 
each  pair  of  rods  shown  in  Fig,  14  in  a manner  as  shown  in  Fig,  12, 

Thus  they  complete  two  similar  parallelograms  each  of  which  has  a 
cooresponding  side  along  a common  line.  By  thus  restraining  one  side 
to  a common  line  the  opposite  sides,  which  are  the  stringers,  are 
restrained  so  as  to  always  be  parallel  to  each  other  although  there  is 
no  other  restriction  on  their  relation  rotational  motion. 


In  order  to  ensure  that  the 
four  bar  linkage  did  not  restrain 
the  movement  of  the  plate  in  the 
direction  perpendicular  to  the 
applied  load  (the  x-direction  as 
defined  by  Fig,  1)  electric  strain 
gages  were  placed  on  the  plate 
to  measure  these  inplane  strains. 

The  test  showed  that  the  linkage 
offered  no  restraint  to  the 
expansion  (or  contraction)  of  the 
plate. 

The  actual  testing  of  pre stressed 
plates  was  carried  out  in  the 
following  manner: 

The  side  edges  were  clamped 
between  the  stringers  'and  the 
plate  flattened  by  moments  applied 
to  the  stringers  by  means  of  monkey 
wrenches.  The  four  bar  linkage 
was  then  attached  to  the  stringers 
thus  forcing  the  stringers  and 
plate  to  lie  in  a plane  after 
removing  the  wrenches.  The  ends 

of  the  plate  were  then  inserted  between  the  steel  angles,  attached  to 
the  end  plates,  and  firmly  bolted  into  position.  The  turnbuckle  connecting 
the  end  plates,  was  then  tightened  until  the  end  plates  were  parallel. 


Detailed  View  of  Four  Bar  Linkage 
Fig.  14 


The  lower  end  plate  rested  on  the  base  of  the  testing  machine  while 
the  upper  end  plate  was  loaded  through  a round  swivel  block  attached  to 
the  measuring  head  of  the  machine.  The  swivel  head  was  adjusted  so  as 
to  be  parallel  to  the  upper  base  plate  of  the  testing  jig. 
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It  was  found  that  If  the  steel  angles  (see  Fig.  12)  were  clamped 
very  tightly  to  the  ends  of  the  plate  the  friction  between  the  angles 
and  the  plate  would  transmit  the  load  uniformly  from  the  end  plates  to 
the  plate  under  test.  This  information  enabled  us  to  forego  any  more 
elaborate  method  of  transmitting  the  load  to  the  plate. 

The  plate  was  loaded  to  about  one  half  the  expected  buckling  load 
and  then  unloaded  to  twenty  pounds  to  insure  the  proper  seating  of  the 
plate  in  the  testing  jig.  The  strain  gage  indicator  was  then  zeroed 
and  the  test  begun.  The  two  sets  of  strain  gages  were  read  alternately 
at  fifty  pound  loading  increments. 

The  loads  measured  during  the  tests  were  the  total  loads  applied 
to  the  testing  jig,  that  is,  the  sum  of  the  loads  applied  to  the  free 
plate,  the  part  of  the  plate  supported  between  the  stringers,  and  the 
stringers  themselves.  This  latter  part  of  the  load  should  be  negligible 
as  the  stringers  were  designed  to  be  flexible  in  the  direction  of 
application  of  the  load.  The  load  going  down  the  supported  edge  strips 
were,  however,  not  negligible  and  were  found  both  experimentally  and 
analytically. 

The  experimental  investigation  of  the  load  carried  by  the  supported 
edge  strip  was  carried  out  by  testing  .040  in.  non-prestressed  flat 
plates  of  the  same  dimensions  as  the  prestressed  plates  but  with  four 
different  widths  of  supported  strips.  The  normal  testing  procedure 
had  a one  inch  edge  strip  supported  by  the  stringers  while,  in  addition, 
for  calibration  purposes,  tests  were  run  with  l/2  in.,  3/4  in.,  and 
1 l/4  in.  strips  supported  by  the  stringers. 

In  one  series  of  tests  the  plates  were  cut  to  size  so  as  to  have 
all/4  inch  strip  supported  by  each  stringer  and  then  repeatedly  trimmed 
down  in  width  so  as  to  have  different  widths  of  supported  edge  strips. 

The  results  of  these  tests  are  shown  in  Table  15, 

In  some  cases,  it  was  observed  that  the  buckling  lead  for  plates 
with  stringer  grip  of  1 1/4  in.  buckled  at  a lower  load  than  plates  with 
1 in.  grip.  The  cause  of  these  seemingly  peculiar  results  is  probably 
due  to  the  variation  of  boundary  conditions,  with  the  width  of  the 
stringer  grip,  and  with  each  test. 

In  another  series  of  tests  the  plates  were  not  retested  but  instead 
the  plates  were  originally  out  to  sizes  for  testing  the  four  different 
widths,  l/2,  3/4,  1,  and  1 1/4  inches.  The  results  of  these  tests  are 
shown  in  Tables  16,  17,  18,  and  19. 
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The  results  of  Tables  16,  17,  18,  and  19  are  shown  graphically  in 
Fig,  15,  Each  buckling  load  (found  by  Yoshiki’s  method)  is  shown  in 
addition  to  their  averages  shown  for  each  initial  radius  of  curvature. 

Two  curves  marked  experimental  are  shown.  They  are  straight  lines 
drawn  through  the  experimental  values  by  the  method  of  least  squares. 

For  one  line,  indicated  by  "Experimental  (l/2  - 1 l/4) " , a least  square 
"average"  is  drawn  through  all  the  points;  while  for  the  other,  indicated 
by  "Experimental  (3/4  - 1 l/4)",  the  points  for  1/2  inch  width  of  plate 
supported  between  each  of  the  stringers  are  omitted.  The  results  of 
the  l/2  inch  supported  width  tests  are  omitted  in  this  one  case  since 
it  is  felt  that  the  clamping  is  too  weak  to  approximate  the  clamped  edge 
condition.  It  can  be  seen  in  Fig.  15  that  the  straight  line  found 
using  the  l/2  inch  supported  width  gives  a "aero  width"  buckling  load 
about  115  lbs.  lower  than  that  found  omitting  the  l/2  inch  width  and  about 
70  lbs.  less  than  that  for  the  ideal  theoretical  plate. 


Width  of  Supported  Strip,  Inches 


Stringer  Calibration  Tests 
(.040  Aluminum  Plates) 

Fig.  15 


Since  it  can  be  seen  that  the 
line  representing  the  buckling 
loads  of  the  theoretical  plate 
falls  well  within  all  the  experimental 
points  (except  perhaps  the  ones 
tested  with  l/2  inch  supported 
edges)  this"theoretieal"  line  will 
be  used  for  calibration  purposes. 

The  line  labeled  theoretical 
is  obtained  by  assuming  that 
supported  edge  strips  undergo  the 
same  total  compression  as  an  ideal 
plate  would  at  the  inception  of 
buckling.  Then  the  load  per  unit 
length  carried  by  the  edge  strip 
will  be  equal  to  the  buckling  load 
per  unit  length  of  the  plate. 

S.  Levy  (12)  has  found  this  for 
aspect  ratio  3:2  to  be: 

Ny  « 8.33 

2 

where  h » .040,  a ■ 9,  v ■ .1, 
ana  E - 10.5  i 106  psl 


X -63.2lb*/ln- 

y 
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Thus  the  buckling  load  for  our  plate  with  a » 9 in.  is  569  lbs. 
and  for  strips  of  the  widths  tested  are  shown  in  Table  20.  Table  20 
is  calculated  by  multiplying  the  buckling  load  per  unit  length  by  the 
width  of  the  supported  strip. 

Now,  if  each  of  these  values  are  added  to  the  buckling  load  of  the 
theoretical  plate  (569  lbs.)  the  values  of  the  buckling  loads  for 
theoretical  plates  with  the  various  supported  widths  (stringer  grips) 
are  obtained.  These  are  the  values  that  determine  the  plot  labeled 
"theoretical". 

Very  close  agreement  between  the  "theoretical"  and  the  experimental 
curve  can  not  be  expected  since  the  calculations  of  the  "theoretical" 
curve  in  no  way  takes  into  account  the  variation  of  degree  of  clamping 
with  the  different  widths  of  stringer  grips ; and  in  an  actual  plate  the 
end  shortening  will  be  greater  than  in  the  theoretical  case.  Therefore 
a larger  load  will  be  carried  by  the  supported  strips  than  originally 
computed.  Moreover,  no  account  has  been  taken  of  load  carried  by  the 
stringers.  From  these  considerations,  the  experimentally  determined 
buckling  loads  will  be  larger  than  those  calculated  from  theory. 

Nevertheless,  in  spite  of  those  objections,  the  plot  of  Fig.  15 
still  gives  a reasonable  indication  of  the  degree  cf  clamping  that  is 
attained  by  the  buckling  jig. 

Having  shown  that,  after  accounting  for  the  load  transmitted  by 
the  strip  of  plate  clamped  between  the  stringers,  the  buckling  loads  are 
close  to  the  values  computed  from  the  accepted  theory,  the  program  of 
testing  prestressed  plates  was  begun  and  is  described  in  the  following 
paragraphs. 

Because  of  the  un symmetrical  buckling  of  prestressed  plate  two  sets 
of  strain  gages  were  used  in  order  to  increase  the  likliehood  of  the  gages 
being  at  the  crest  of  the  buckling  modes. 

The  results  of  testing  of  non-prestressed  flat  plates  were  very 
consistent  and  the  results  of  a representative  test  is  shown  in  Fig.  16. 
Both  (cj-Gg)  and  are  plotted  and  since  these  are  proportional 

to  deflection  and  the  square  of  deflection  respectively  the  methods 
previously  described  may  be  used  in  determining  the  buckling  load.  The 
plate  shown  is  numbered  F40-1,  5-H  where  F indicates  that  it  was  an 
non-prestressed  flat  plate  (prestressed  plates  are  indicated  by  the  letter 
C),  the  first  set  of  numbers  indicates  the  plate  thickness  to  be  .040 
inches,  the  second  set  of  numbers  indicates  that  the  aspect  ratio  was 


WADC  TR  54-6 


28 


Typical  Load-Deflection  Diagram 
for  Ordinary  Flat  Plate 

Fig.  16 


"Normal"  Type  Buckling  of 
Prestressed.  Plate 

Fig.  17 


1.5,  while  the  third  series  of 
numbers  indicates  that  this  was  the 
eleventh  plate  of  this  type  tested. 

The  plots  of  the  test  data 
of  prestressed  flat  plates  were 
of  three  general  types  shown  In 
Figs.  17,  18,  and  19.  The  curves 
of  Fig.  17  are  what  will  be  called 
the  "normal"  type,  that  is,  they 
are  of  the  type  expected  theoretically. 
The  curves  of  Fig.  18  will  be  called 
the  "double  parabola"  type  in  that 
the  post  buckling  behavior  is 
characterised  by  two  parabolas. 

Instead  of  one.  Thus  if  plates 
of  the  "double  parabola"  type  are 
analyzed  by  means  of  the 

plot  there  will  ire  suit  two  different 
buckling  loads  corresponding  to 
the  two  parabolic  post  buckling 
behaviors  of  the  plate.  The 
significance  of  these  two  parabolas 
is  not  obvious,  and  is  a puzzling 
characteristic  of  pre stressed  plates. 
The  curves  of  Fig.  19  will  be 
called  the  "mode  jump"  type,  that 
is,  they  occur  when  the  plate 
buckling  form  changes  suddenly  during 
the  test.  The  change  is  from  two 
half-waves  to  three  half-waves 
and  is  sometimes  so  sudden  as  to 
be  accompanied  by  a loud  popping 
sound. 

' Flat  plates  with  an  aspect 
(length  to  width)  ratio  of  1.5 
should,  according  to  theory,  buckle 
in  two  half  waves.  It  is  felt 
that  initial  imperfections  are  a 
major  cause  for  the  formation  of 
three  wave  form.  If  the  plates 
had  longitudinal  initial  deflections 
of  the  form  shown  in  Fig.  20a 
which  is  shown  to  be  the  case  for 
prestressed  plates  (see  Fig.  20), 
the. plate  could  be  forced  to  buckle 
into  three  instead  of  two  half 
waves. 
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"Double  Parabola"  Type  of  Buckling  "Mode  Jump"  Type  of  Buckling 
of  Prestressed  Plate  of  Prestressed  Plate 

Fig.  18  Fig.  19 

The  hypothesis  can  be  visualised  with  the  help  of  Fig.  20.  If  the 
initial  deflections  of  one  half  wave  shown  in  Fig.  20a  are  added  to  the 
usual  buckling  form  of  two  waves  shown  in  Fig.  20b,  an  unsymmetrical  two 
wave  form,  shown  in  Fig.  20c,  results.  This  wave  form  is  actually  a 
characteristic  of  the  prestressed  plates.  It  is  possible  that  the  longer 
half  wave  of  Fig.  20c  becomes  unstable  and  this  buckles  into  two  more 
waves  giving  a total  of  three  half  waves  shown  in  Fig.  20d. 


(o)  (b)  (c)  (d) 

Effect  of  Initial  Deflections  on  Buckling  Mode 
Fig.  20 
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The  above  hypothesis,  although  conjectural,  is  somewhat  substantiated 
by  experience  with  non-prestressed  flat  plates.  During  the  testing  of 
non-prestressed  flat  plates  it  was  observed  that  the  plates  "popping" 
into  three  waves  were  those  with  the  largest  initial  deformations  in  the 
form  of  one  wave.  After  noticing  this  the  testing  of  non-prestressed 
plates  were  restricted  to  those  that  appeared,  to  the  naked  eye,  to  be 
flat.  Thereafter,  no  trouble  was  encountered  with  this  "popping"  phenomenon. 
Unfortunately,  initial  deflections  are  inherent  in  the  prestressing  process, 
thus  necessitating  a larger  number  of  tests  to  compensate  for  the  tests 
discarded  because  of  "popping". 

Another  possibility  for  explaining  this  "popping"  is  the  formation 
of  inplane  stresses  that  would  be  of  such  a nature  as  to  cause  the  plate 
to  buckle  into  three  half-waves.  The  probability  of  this  hypothesis  being 
a major  factor  is  strengthened  considerably  by  the  results  of  the  Rayleigh- 
Ritz  anallses  shown  in  the  previous  section. 

It  was  shown  that  for  some  stress  distributions  occurring  during  the 
prestressing  process  (particularly  that  found  by  the  experiments  performed 
on  the  second  ,040  inch  thick  Aluminum  plate)  the  plate  will  buckle  into 
three  half-waves.  In  some  plates  (particularly  the  ,020  inch  thick  aluminum 
plate)  the  energy  level  at  which  the  plate  will  buckle  is  very  nearly  the 
same  for  the  symmetric  or  antisymmetric  modes  of  failure.  It  was  also  shown 
that  there  is  a wide  scattering  in  the  stress  distribution  found  experimentally. 
Since  there  will  also  be  some  variation  in  the  degree  of  clamping  at  the 
boundaries  and  some  non-homogeneity  of  the  plate  and  also  since  .the  prestress 
distribution  will  vary  somewhat  upon  deflecting  under  the  applied  buckling 
load  it  is  quite  likely  that  in  some  cases  a plate  will  buckle  in  two  half- 
waves while  in  other  cases  an  apparently  similar  plate  will  buckle  in  three 
half-waves. 

The  results  of  buckling  tests  performed  on  non-prestressed  flat  plates 
of  ,040,  .032,  and  ,051  inch  thickness  are  shown  in  Tables  18,  21,  and  22, 
respectivley.  The  results  of  buckling  tests  performed  on  prestressed  plates 
of  thicknesses  of  ,040,  ,032,  ,051  inches  are  shown  in  Tables  23,  24,  and 
25,  respectively, 

The  results  of  the  buckling  tests  performed  on  non-prestressed  flat 
plates  of  ,032  and  ,051  inch  thicknesses  which  were  listed  in  Tables  21  and 
22  are  shown  graphically  in  Figures  21  and  22  respectively.  These  figures 
are  similar  to  Fig.  15  which  showed  the  tests  of  .040  non-prestressed  flat 
plate  tests  for  calibration  purposes. 
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Since  the  calibration  tests  showed  that  very  good  results  were  obtained 
by  assuming  that  the  supported  edge  strips  undergo  the  same  total  compression 
as  an  ideal  plate  would,  the  calibration  tests  were  dispensed  with  in  the 
case  of  plates  of  *032  and  *051  inch  thicknesses* 
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Width  of  Supported  Strip,  Inches 


*032  Ordinary  Flat  Plate  Tests 
Fig.  21 


•051  Ordinary  Flat  Plate  Tests 
Fig.  22 


The  theoretical  buckling  loads  of  .032  and  .051  inch  thick  non-pre stressed 
flat  plates  with  different  widths  of  supported  edge  strips  are  shown  in 
Figs.  21  and  22  along  with  the  experimentally  determined  values  (by  Yoshiki's 
method).  The  theory  indicates  that  for  .032,  .040,  and  .051  inch  plates  the 
load  supported  by  the  1 inch  supported  strips  at  both  edges  will  be  63,  126, 
and  262  pounds  respectively. 

Making  use  of  these  values  and  subtracting  them  from  the  loads  listed 
in  Tables  23,  24*  and  25  ( and  shown  graphically  ih  Figs.  23,  24,  and  25) 
the  net  buckling  loads  were  determined.  The  average  experimental  "total" 
buckling  loads  are  listed  in  Tables  26,  27,  and  2 8 for  the  .040,  .032,  and 
•051  inch  thick  prestressed  plates  respectively  while  the  net  loads  are  listed 
in  Table  29. 
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Buckling  Loads  of  *040  Prestressed 
Plat  Plates 


Buckling  Loads  of  .032  Prestressed 
Flat  Plates 


Fig.  23 


Fig.  24 


Buckling  Loads  of  .051  Pre stressed 
Flat  Plates 

Fig.  25 


From  the  data  listed  in  the 
above  tables,  it  is  felt  that 
Yoshikl’s  method,  though  convenient 
to  use  for  non-Jrre stressed  plates, 
may  not  be  applicable  to  the  pre- 
stressed plates  as  indicated  by 
the  occurrence  of  double  parabola 
in  the  plate.  As  a result,  the 
top  of  the  knee  method  was  used. 

The  buckling  loads  thus  determined 
indicate  the  loads  at  which  the 
deflections  suddenly  become  large. 
The  net  loads  supported  by  the 
"free"  plate  are  listed  in  Table  29 
together  with  their  percentage 
increases  over  the  buckling  loads 
of  the  corresponding  non-prestressed 
plates,  both  measured  by  the  top 
of  knee  method.  For  purpose  of 
comparison,  the  values  determined 
by  Yoshiki’s  method  are  given  in 
the  parentheses  after  the  net 
buckling  loads  in  Table  29* 
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SECTION  IX 


Remarks  on  the  "Buckling"  Criteria 

The  current  methods  of  testing  make  use  of  the  form  of  the  Load-strain 
diagram  rather  than  the  actual  values  of  the  deflections  in  the  diagram* 

These  methods  all  have  as  a basis  the  load-strain  relations  existing  in  a 
theoretically  flat  plate  which  does  not  deflect  until  the  buckling  load  is 
reached  at  which  point  the  deflections  suddenly  occurs.  The  actual  plate 
with  initial  deflections,  does  not  esdiibit  this  sudden  increase  in  deflection 
and  in  fact  in  some  cases  the  evidences  of  the  behavior  of  the  ideal  plate 
is  indiscernable.  As  was  mentioned  in  the  Introduction,  a factor  of  major 
significance  is  the  ratio  of  plate  deflection  to  applied  load.  It  can  be 
seen  from  Fig.  26  that  a lower  value  of  this  ratio  does  not  necessarily 
occur  with  a higher  buckling  load  of  the  equivalent  ideal  plate.  In  this 
figure  two  Load  - deflection  diagrams  are  shown  along  with  their  extensions 
(shown  by  dashed  lines)  which  would  represent  the  diagram  of  the  equivalent 
ideal  plate.  The  intersection  of  this  dashed  line  with  the  ordinate  is 
the  buckling  load  which  all  of  the  measuring  methods  attempt  to  find.  It  can 
be  seen,  however,  that  in  the  case  shown  the  plate  with  the  highest 
theoretical  buckling  load  also  has  the  largest  deflection  to  load  ratio, 
a distinctly  undesirable  characteristic. 

Because  of  the  relatively  large  deflections  inherent  in  prestressed 
plates  the  importance  of  low  values  of  Deflection  to  Load  ratios  is 
especially  important.  It  is  felt  that  the  "buckling  load"  in  such  a case 
should  be  defined  as  the  load  at  which  the  maximum  deflection  in  the  plate 
reaches  a certain  prescribed  value.  Further  research  on  prestressed  plates 
with  this  new  "buckling"  criteria  would  seem  to  be  of  great  value. 


Illustration  of  Buckling  Criterias 
Fig.  26 
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FORMULATION  OF  THE  PROBLEM 

The  theoretical  treatment  of  the  problem  consists  of  two  parts. 

First,  it  is  necessary  to  find  the  inplane  stresses  induced  in  the  plate 
due  to  the  clamping  of  a curved  plate  into  a flat  frame.  Next,  with 
these  inplane  stresses  determined,  the  buckling  load  will  be  calculated 
for  a flat  plate  with  such  initial  stresses. 

The  first  problem  can  be  formulated  by  considering  the  process  of 
clamping  the  curved  plate  into  a flat  frame  to  occur  in  three  steps. 

First,  the  cylindrically  curved  plate  is  bent  into  an  identically  flat 
form  by  pure  moments  applied  at  the  edge.  Secondly,  the  edges  of  the 
plate  in  this  identically  flat  form,  will  be  clamped.  Note  that  the 
clamping  will  induce  no  additional  deflections  or  stresses. 

Thirdly,  the  pure  moments  applied  in  step  one  will  be  removed  by 
the  application  of  equal  and  opposite  moments  thus  leaving  the  clamping 
action  at  the  edges  as  the  only  restraining  force  on  the  plate. 

It  may  be  pointed  out  that  initial  bending  stresses  do  not  affect 
the  buckling  load  of  plates  but  inplane  (or  membrane)  stresses  do.  Since 
the  first  step  produces  only  pure  bending  stresses  and  the  second  step 
induces  no  stresses  in  the  plate,  our  analysis  will  be  concerned  exclusively 
in  analyzing  the  deflections  and  stresses  induced  by  the  loadings  of  the 
third  step.  This  analysis  will  yield  the  deflections  and  inplane  stresses 
ocurring  when  a cylindrically  curved  plate  is  clamped  into  a flat  frame. 

The  problem  of  determining  the  buckling  load  for  a prestressed  flat 
plate  problem  can  be  solved  by  the  use  of  the  Rayleigh-Ritz  energy  method. 

In  using  the  energy  method,  it  will  be  assumed  that  when  buckling  occurs 
these  inplane  stresses  will  remain  unchanged. 


SECTION  XI 

BENDING  OF  A CYLINDRICALLY  CURVED  PLATE  INTO  FLAT  FORM 

If  the  curved  plate  originally  is  in  the  form  of  a cylindrical  surface, 
it  is  possible  to  bend  it  into  an  identically  flat  form  by  inextensional 
bending.  This  is  because  both  the  cylindrical  and  flat  surfaces  are 
developable  surfaces.  The  problem  of  bending  one  developable  surface  into 
another  one  can  be  found  in  many  texts  on  Elasticity,  for  example,  ref.  13. 
The  edge  moments  necessary  to  bend  such  a plate  (Fig.  27)  are 


II 

w 

(-  ♦ V i-) 

R R 

(43) 

X 

y 

M = D 

(A-  .,1) 

(44) 

y 

•V 

Rx 

WADC  TR  54-8 


35 


are  the  changes  in  curvature . In  our  case  — * 0 and 


Mx  » vD/R  (45) 

Hy  - B/R  (46) 

In  order  to  apply  equal  and  opposite  moments  to  the  plate  when  the 
plate  has  been  clamped  at  the  edges,  it  is  necessary  to  represent  these 
moments  in  the  form  of  equal  and  opposite  pressures  applied  near  the  edges 
of  the  plate  (see  Fig.  28).  This  can  be  done  by  expressing  the  pressure 
distributions  in  terms  of  Dirac  delta  functions  (14).  According  to  usual 
notation  for  these  delta  function^  c f (x-e)  indicates  the  delta  function 
which  intersects  the  abscissa  at  x = e.  <?’  indicates  the  first  derivation 
of  <$, . 


where 


and 


R 


Thus 


R 


Coordinate  System 
Fig.  27 


Representation  of  Edge  Moments 
by  Pressure  Couples 
Fig.  28 
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The  pressure  distributions  cf(x)  and  <h  (y)  Bust  have  the  following 
properties: 


j S^fr-e)  dy  « 
0 

V - 

V <$2(y-b+c)  dy  = 0 (47) 

0 

a 

f <f*^(x-e)  dx  * 

(V*2(x- a+c)  dx  ■ 0 

/ 

0 

/ 

0 

V 

J (y-c)  /^(y-e)  dy  = - 
0 

v 

\ (y-b+c)  cf»2(y-b+c)  dy  « M^ 
0 

b*  r 

i (x-e)  o’^(x-c)  dx  “ - 

V f 

r (x-a+e)  tf»2(x-a+c)  dx  « ^ 

0 

/ 

0 

In  addition,  the  derivative  of  the  Dirac  delta  function  has  the 
following  useful  integration  properties: 

J f*(y)  /^(y-c)  dy  = Mx  tlistSLlS'lss^.  ~ ^ f»(c) 

0 


V 

J f(y)  </*2(y- b+c)  dy  » - Mx  tlfa-ZrP)  f»(b-t) 

0 2 

(49) 

a 

/ f(x)  {T*1(x-c)  dx  = My  ^My  f»(c) 

0 

f(x)  /*2(x-a+c)  dx*  tllS=S^L±  H1S^=Q1  * . ^ f '(a_e) 
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This  delta  function  pressure  distribution  can  be  written  in  terms 
of  a Fourier  series  and  the  Fourier  coefficients  easily  found  by  using 
the  above  integration  properties  of  Dirac  delta  functions*  Let 


«/*  (x)  » ^(x-c)  + </*2(x-a+e) 

<f*(y)  ■ /^(y-e)  + <f*2(y-b+e) 

d*  may  be  written  as  Fourier  series  by  noting  that 

cf'(x)  - ~ + £ a cos  — + 2 bm  sin  — 

2 m a n a 

<f*  (y)  *=  — + So  cos  222  + 2d.  sin  222 


(50) 


where 
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v " “ r>«  °os  ~ dx 
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sin  — dx 


) cos  222  dy 
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0 
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dn  ■ - jV*(y)  sin 


222  (Jy 

b 
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(51) 


Making  use  of  the  integration  properties  of  the  Fourier  coefficients 
(51),  may  be  evaluated. 
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Similarly, 


(52) 


c_  * 0 
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b«*f  p 


(x-e)  sin  dx  + — 
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(x-a+e)  sin  dx 
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a |_  X a a X a a J 


2mr  P_  *7 

~2  Mx  “ 003  m 7rJ 


b * 0 for  even  m 
m 


and  bm  * ^ for  odd  m 


similarly 


d_  * 0 for  even  n 
n 


d « M for  odd  n 
n b2  7 


(53) 
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Substituting  (50)  and  (53)  in  eqs.  (50)  we  find 


<f»(x)  * Z mainSS 

sr  7 m-1,3,5  a 


5(y)  • ^ M Z n sin  222 
b*  x n-1,3,5  b 

We  will  now  proceed  to  use  the  above  results  in  determining  the 
Fourier  series  representing  the  applied  pressure,  p , which  in  turn 
approximates  the  effect  of  a moment  equal  and  opposite  to  the  moments 
and  M required  to  open  the  plate  to  an  identically  flat  form* 


Sy  definition. 


Pe  ■ <Wx)  + £'(y) 


Substituting  eqs.  (45) » (46),  and  (54)  in  (55),  we  find 

p-M  (r  SL  sin  222  + £ 2*  ain  2227 

7 l mKL,3,5  a*  a n=l,3,5  b*  b ) 


Since  it  is  more  convenient  in  the  later  analysis  to  use  a double 
Fourier  series,  we  shall  expand  4frH^/a2  and  4n’Myv/b2  into  sine  series  as 


follows ! 


4tfM 

— ^ * Z e sin  522 
a s=l  b 


4rrM  v 


Z f sin 
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It  is  easy  to  verify  that 


- for  odd  s 
a 


even  s 


- for  odd  r 
r 

even  r (58) 


P-16M  2 2 (£  1 + £ v^)  8ln  2S  8in  222 

e y r»l,3,5  s*l,3,5  s z r 1>  a b 

(59) 

In  addition  to  the  above  presatire  equivalent  to  the  applied  forees 
(hereafter  referred  to  simply  as  the  equivalent  pressure),  there  will  also 
be  a clamping  pressure,  p , which  will  be  the  equivalent  of  the  moments 
oc curing  at  the  edges  of  the  damped  plate.  This  pressure  will  be  expanded 
in  a trigonometric  series  containing  arbitrary  coefficients  ifcich  will  be 
determined  by  the  use  of  the  boundary  conditions  due  to  the  clamping  of 
the  edge. 

The  clamping  effect  will  be  replaced  by  equivalent  edge  moments  ny 
and  ay  (at  i*  0,  a and  y * 0,  b respectively).  Using  Eq.  (56) 

.o  4^*  ___  «o  4^m_ 

p * 2 — r sin + 2 — 8 sin  ^2 

r-1,3,5  » a e-1,3,5  b*  b 

(60) 

m and  m may  be  expressed  by  a trigonometric  series  in  terms  of  arbitrary 

* y 

coefficients  K_  and  T_ 

m = vM  ^ 2 T sin 23 

7ir  r«l,3,5  a 


n - M ^ 2 T sin  521 

7 7 ir  s«l,3,5  8 b 


16M 


0 for 

16M  v 
2. 

b2 
0 for 


Hence, 
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combining  Eqs.  (60)  and  (6l)  gives  the  equivalent  clamping  pressure 

00  -T  K - 

p - 16M  2 2 (-§  r + v -§  a)  ain  — sin  sa 

° y r-1,3,5  s»l,3,5  & b*  a b 

(62) 

The  addition  of  the  equivalent  applied  pressure,  p , and  the 
equivalent  clamping  pressure,  pc,  will  be  called  the  total  applied 
pressure,  p^*  Thus, 


p-  - p + p « l6M  2 2 (T  + -)+  v^K  + -)  sin^^sin22 

T e c y r-1,3,5  s-1,3,5  * B s b*  r r a b 


(63) 

To  simplify  the  writing,  let 

p-  - 2 *2  p„  a sin  22  sin  522  (64) 

1 r-1,3,5  s-1,3,5  r,s  a b 


where 


It  should  be  pointed  out  that  (54)  are  divergent  series  and  conse- 
quently Eqs*  (59),  (60),  (62),  and  (65)  consist  of  divergent  series*  The 
paradox  of  obtaining  a convergent  solution  to  an  equilibrium  problem  under 
the  action  of  a load  represented  by  a divergent  series  has  already  been 
noted  by  S*  Levy  (15). 

The  explanation  of  this  paradox  would  require  further  mathematical 
research*  It  can  be  shown  that,  for  instance,  if  the  divergent  represen- 
tation (Eq.  65)  for  the  applied  pressure  is  summed  according  to  the  Borel 
summation  (one  of  many  types  of  summations  useful  in  summing  divergent 
series)  it  is  indeed  found  to  be  zero  at,  for  instance,  the  center  of  the 
plate,  as  it  should  he*  The  physical  meaning  behind  such  a means  for 
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summing  divergent  series  seems  rather  nebulous  and  no  concrete  conclusions 
can  be  drawn  from  such  a summation  yielding  the  desired  results* 


The  view  that  his  paradox  is  not  an  accidental  one  is  strengthened 
considerably  by  noting  its  occurrence  in  a rather  simple  linear  problem* 

The  pure  bending  of  a simply  supported  beam  of  length  "a"  by  equal 
moments.  My.,  applied  at  its  ends  can  be  solved  by  replacing  these  moments 

by  equivalent  pressures  exactly  in  the  same  manner  as  has  been  done  in  the 
previous  section  for  the  clamped  plate*  Solving  this  problem  by  the  method 
under  question  will  lead  to  the  exact  solution* 

Using  equation  (54)  the  moment  equivalent  pressure  may  be  written 
as 

^ *V  ^ 

eq  & 7 art, 3, 5 . . a 

The  equilibrium  equation  for  a beam  under  a pressure  load  is: 


For  our  problem  the  boundary  conditions  at  x « 0 and  x * a are 

2 

w * 0 and  —r  B 0 
dx* 

These  boundary  conditions  are  satisfied  termwise  by  assuming  the 
deflection  in  the  fomr  of  the  following  infinite  series: 


w * Z sin  s-=t 
m-1  m a 

Substituting  the  above  two  infinite  series  into  the  equilibrium 
equation,  we  obtain 


I Z v (SL)  8ln2HS-„4|M  ? 

art  a a a 7 m®l,3»5j 


m sin 
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Equating  coefficients  we  find 


v»'-^  i 


or 


v‘7 


2 i 1 3Lsln=S 

El  iplf3|5i  • • m a 


For  purposes  of  comparison  we  find 


-2 

d w 

. 2 
dx 


» M <*> 

St  _2  y 


ir  El  m»l,3,5f  • • m 


Making  use  of  the  following  formula  No.  808  listed  by  B.O.  Pierce  (17): 


V4 


1 .toHS 


mFl,3,5»  • • m 


we  find 


.2  M 
d w j 

dx2  El 


which  is  well  known  to  be  the  exact  solution  for  the  second  derivative 
of  a beam  bent  by  pure  moments. 


The  success  of  using  a divergent  representation  for  the  pressure, 
which  arises  from  the  use  of  the  Dirac  delta  function,  seems  to  be  hinged 
upon  the  order  of  the  differential  equation.  Note  that  if  the  above  equation 
were  a first  order  differential  equation,  e.g. 


dx 


the  solution  for  w would  then  be  expressed  by  a divergent  series  and  the 
solution  could  not  be  correct* 
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SECTION  m 


THE  DETERMINATION  OF  THE  INPLANE  STRESSES 

Next,  we  shall  calculate  the  inplane  stresses  in  a flat  plate  clamped 
at  four  edges  and  under  some  distribution  of  lateral  pressure  which  is 
equivalent  to  the  edge  moments  applied.  In  this  case,  the  von  Kar man's 
large  deflection  equations  for  thin  plates  (19)  will  be  used.  They  are 

1 ( ay  ax  ay  axay  axayj 

vV-E  ((&-)  4»)7  (67) 

L axay  ax  ay*"  J 

where  the  pressure  p^,  is  the  equivalent  pressure  derived  in  the  previous 

section  and  the  stress  function,  F,  is  defined  in  terms  of  the  inplane 
stresses  by  the  following  relations: 

a2?  a2?  a2? 

‘SO  ‘ ’Vi  ■ 8J2  > VO  ' • ^ 

The  boundary  conditions  on  the  deflection,  w,  are  as  follows: 


at  x ■ 0,a 


at  y * 0,b 


The  boundary  conditions  on  the  stress  function  F are  such  that 

1,  The  resultant  load  (at  any  cross  section  of  the  plate)  must 
be  zero  in  the  x-direction  and  in  the  y-direction. 

2,  The  boundaries  of  the  plate  must  remain  straight. 

(69) 
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The  deflection,  w,  and  the  stress  function,  F,  will  be  assumed  in 
the  form  of  the  following  trigonometric  series  with  arbitrary  coefficients! 


oo 

v * 2 

ao 

2 v 

. rsnrx 
sin 

sinSa 

ml 

n-1  m»n 

a 

b 

ao 

F - S 

OO 

cos  ^22 

cos 

p»0 

q=0  P,q 

a 

b 

(70) 


(71) 


Note  that  the  pressure  p^,  has  already  been  expanded  in  the  form  of  a 
Fourier  series,  namely. 


PT  « 2 2 pT 

r=l,3,5  s=l, 3,5  1 


sin£S  sln  SSL 


(72) 


To  satisfy  the  differential  equations,  the  coefficients  b and 

P»<1 

wm  are  to  be  determined  as  follows.  We  first  substitute  (70)  and  (71) 

into  equation  (67).  By  equating  the  coefficients  of  like  terms,  b can 

Pt^ 

be  expressed  in  terms  of  w in  the  following  form. 


£ 


9 

"P,q  4(p2  - + q2  S)2  i=l  B± 

a b 


(73) 


where  are  functions  of  w^  and  are  listed  in  Appendix  B.  Next, 

substitute  (70),  (71),  (72)  into  equation  (66)  and  equate  the  coefficients 
of  like  terms.  Since  b are  now  functions  of  w , we  obtain 

p»q 


, 2 7T2  2 v2,2  hr4 

^.=(r  7 ? 777 


i=l 


r,s 


(74) 
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where  A.  are  function  of  w _ and  are  listed  in  Appendix  C.  The  derivation 
1 m,n 

is  the  same  as  that  carried  out  by  Levy  (20).  Note  that  in  this  problem 
due  to  symmetry  the  coefficients  p„  _ and  w„  _ are  restricted  to  odd 

subscripts  only.  Because  of  this,  it  is  found  that  b are  restricted 
to  even  subscripts  only. 

Before  we  attempt  to  carry  out  the  solution  of  wm  let  us  examine 
the  boundary  conditions.  * 

It  is  obvious  that  the  assumed  series,  (70),  is  termwise  identically 
zero  at  x * 0,  x * a,  y * 0,  and  y * bj  hence  the  first  set  of  the  boundary 
conditions  on  w is  satisfied.  Also 


— - - Z Z m w cos  — sin  221 

0x  a m=l  n=l  m»n  a b 

and 

0w  _ tt  ~ to  my 

5y  b m=l  n=l  ,n  a b 

0W 

To  satisfy  the  conditions  — * 0 at  x * 0 and  at  x * s,  we  must  have 

0x 


or 


TT 

0O 

Z 

oQ 

z 

a 

m«l 

n=l 

TT 

oo 

OO 

z 

Z 

a 

m=l 

n»l 

90 

z 

m v 

at 

m»l 

O o 

m,n 

m w sin  » 0 
m,n  b 

m(— l)  mw  sin  — — ^ 

m.n  . 


0 


Z 

npl 


m(-l)m  w 
^ m,n 


0 


0 


(75) 
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The  conditions  of  zero  slope  at  the  edges  y * 0,  b lead  to  similar 
equations 


Z 

n»l 


n w 


m,n 


Z 

n=l 


n(-l) 


n w. 


m.n 


0 


(76) 


The  two  relations  in  (75)  and  also  those  in  (76)  become  the  same 
if  m and  n are  either  odd  or  even. 

In  our  problem,  m and  n must  be  odd,  the  boundary  conditions  of 
zero  slope  therefore  lead  to  the  following  requirements  on  the  deflection 
coefficients 


Z m K,  » 0 

®=1,3,5  * 

(77) 

OO 

Z n w = 0 

n-1,3,5  m’n 

(78) 

Now  let  us  examine  the  boundary  conditions,  (69). 

The  first  condition  may  be  shown  to  be  satisfied  in  the  following 
manner.  Defining  the  load  in  the  plane  of  the  plate  as  px  in  the  x- 

direction  and  p^  in  the  y-direction.  Then 

b. 


/ 


h 


where  h is  the  thickness  of  the  plate  and  is  a constant.  But 


„ a2? 

v 


Then 


b 

/ 


£2 

8y2 


dy 


(8£)  _ (M) 

8y  8y 

y»b  y*0 


] 
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Now, 


F-  s ? b cos  ES  cos  322 

p*0  q»0  p,q  a b 


therefore  Px  “ 0 


Proceeding  in  a similar  fashion,  we  find 

- 0 


py-h  / 


a o . b 

£§  dx-h^ 

ax  ax 

o 


(79) 

(80) 


The  second  boundary  condition  can  be  shown  to  be  satisfied  as  follows* 
Displacement  of  edges  in  x-direction. 


Displacement  of  edges  in  the  y-direction 


(81) 


(82) 


where  r and  e are  the  inplane 
x y 

in  terms  of  the  stress  function 


strain  components  which  may  be  written 
as  follows 


.1 

e ay  ax 


y E dx  Qj 


(83) 
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above  relations 


(64) 


It  is  evident  that  these  relations  are  independent  of  x and  y and 
consequently  the  conditions  of  straight  edges  are  satisfied.  Since  the 
detail  of  calculation  is  not  obvious  and  is  not  included  in  Levy's  reports, 
it  is  included  in  Appendix  D. 

Having  examined  the  boundary  conditions,  let  us  now  proceed  to  solve 
for  the  coefficients  vm  n*  Since  the  problem  was  set  up  in  a manner 

analogous  to  that  of  Levy  amd  Greenman  (21),  much  of  the  numerical  work 

done  by  them  may  be  used  in  this  investigation.  First,  let  us  express 

the  coefficients  b in  terms  of  w_  • Tb is  is  listed  in  Table  30. 
p,q  m,n 

Next,  let  us  substitute  b in  terms  of  w_  and  p as  given  by  (65) 

* p,q  m,n  *r,s 

into  (74)*  We  obtain 

Wi^A  * (.089783  + ,011057K1,  + .078726^)  a2/Hh  - 
Wi>3A  ■ (.0049583  + .0027683^  + .0065702T3)  a2/Rh  - Dg/h3 

w3>1A  - (.0056107  + .00025863K3  ♦ .0055245^)  a2/Rh  - Dy^h3 

w3>3A  * (.0011085  + .OOO4O949K3  + .0029158T3)  a2/Hh  - D^A3 

w1>5A  - (.0010103  + .000786371^  + .0011198^)  a2/fch  - D^h3 

w5>1A  * (.0012757  +.  000035632K5  + .0012685^)  a2/Rh  - D^A3 


Substituting  the  expressions  for  w and  F into  the 
and  then  combining  terms,  we  find  finally 


ir  _ _ 22 

— 2 2 m V 

8a  bfI  n«l  * 


2 2 n2w  2 

7 8b  bfI  n»l  *n 


if  * - 
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w?>1A  - (.00047166  + .0000094360^  + .00047031^)  a2/Rh  - lyfe3 

+ 9 o 

“9,1A  * < .00022325  + .00000347781^  + .00022287^)  a /Rh  - DgA 
wn>1A  * (.00012265  + .0000015640K11  + ,00Q12250r1)  a2/Rh  - 
W13,lA  " ( .°00074433  + ,00000080346K13  + .0000743721^)  a2/Rh  - I>10A3 
w15>1A  « (.000048506  + .00000045389K15  + .0000484751^)  a 2/fch  - D^A3 
v5>3A  * (.00034197  .000082290K5  + .00097655T3)  a 2/fch  - Dj^ii3 

v?>3A  - (.00013996  ♦ .000024629^  + .000409323T3)  a2/Hh  - D13A3 
w9>3A  - (.000069267  + .0000095788!^  + ,0002Q462T3)  a2/fch  - D-^A3 
w11j3A  » (.000038948  + ,0000044292Ku  + .00011564T3)  a2/fch  - D^A3 
w13,3A  * (.000023969  + .0000023123^  + .000071346T3)  a 2/Bh  - D^A3 
w3>5A  - (.00033873  + .00028518K3  + .0012183T5)  a 2/Rh  - I>17A3 
w5>5A  - (.00014365  + .000088453K5  + .00062981T5)  a2/hh  - DlgA3 
K^jA  (.000068202  + .000031922^  + .00031821T5)  a2/Rh  - D19A3 
w9>5A  = (.000036358  + .0000135951^  + .00d7424T5)  a2/Hh  - D2QA3 

51 
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W11  5^  * ( *000021305  + .0000066087Kn  + .00010352T5)  a /Rh  - D21A 
w13  5^  " (.OOOO^S  + ,0000035555K13  + .000065824^)  a2/Rh  - D22A3 
W;L  ?A  - (.00035647  + . 000311 24^x  + .0003l659T7)  a2/Rh  - D23A3 
v3  ?A  - (.00013114  + .00017047K3  + .00052020T?)  a2/Rh  - D^A3 
w$  ?A  “ (.000068377  + .000073798K5  * .00037 533T?)  a2Ah  - D^A3 
?A  * (.000037394  + .0000322361^  + .00022953T?)  a2/Rh  - D,^3 
w9  ?A  - (.000021690  + .0000152871^  + .00013994T?)  a2/Rh  - D^A3 
W11  7^  * (.000013382  + .0000079213^  + .000088632T?)  a2/Rh  - D2gA3 
wx  ^ » (.000164988  + .00015166  + .00011998T9)  a2/Rh  - D.^3 

v3  ^ * (.000061214  + .00010253K3  + .00024334^)  a2/Rh  - D3QA3 
w5  ^ - (.000035683  + .000055796K5  + .00022071T9)  a2/fch  - D3].A3 
w?  9A  = (.000021786  + .0000287371^  + .00015914^)  a2/Rh  - D32A3 
W9,9^-  (.000013684  + .0000151671^  + .00010799T9)  a2/Rh  - D33A3 
w13  7A  * (.0000087224  + .0000044368^  + .000058668T7)  a2/Rh 
wn  <^A  * (.0000089101  + .000008423aKu  + ,000073302T9)  a2/Rh 
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v13  ■ (.0000060254  + .000004940^  + .000050810T9)  a2/Rh 

wl,ll/^  = C • °00089546  + .0000845691^  + .000054742Tn)  a2/Rh 
w3>uA  * (.000032830  + .000064388K3  + .00012504TU)  a2/Rh 
w5,uA  ■ (.000020208  + ,000040889K5  + .00013233Ti;l)  a2/Rh 
W7,uA  = («0°«>13327  + .0000240221^  + ,00010885T1;l)  a2/Hh 
w9>nA  = (.0000089504  + .0000139701^  + ,000081382TU)  &Z/Bh 
Wll,llA  * (.0000061323  + .0000083070^  + .000059148TU)  &2/Bh 
w13,llA  = (.0000043047  + .00000511301^2  + .000043024TU)  &2/Bh 
v1>13A  - (.000053951  + .0000517701^  + .000028355^)  a2Ah 
w3>13A  = (.000019488  + .000042391K3  + .00006954T13)  a2/Rh 
W5,13A  * (.000012288  + .0000299221^  + .000081944^)  a2/Rh 
w7>i3A  = (.0000085232  + .0000194691^  + .000074644^)  a2/Rh 
w9,13A  = (.0000060332  + .000012306^  + .000060659T13>  a2/Rh 
wll,13^1  * (.0000043226  + .0000077975^  + . 000046979T^3 ) a2/Rh 
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w13>13A  » (.0000031436  + .0000050327K13  + .000035834*^)  a2/Rh 

(85) 

where  the  expressions  D2,  . . • , D33  are  given  in  Table  31,  with 
the  assumption  that  the  non-linear  terns  not  involving  v,  -,  w.  ,,  v0  -, 

XfX  XpJ  JpX 

W3  3,  w^  ^ and  ^ are  negligible  as  compared  to  terms  involving  only 

these  six  terms.  Table  31  should  be  read  in  such  a manner  that  * 1.451 

wl,l3  “ #972  wl,l  wl,3  " ,192°  wl,l  w3,l  ■*■•••• 

In  order  to  satisfy  the  boundary  conditions  (68)  on  the  slope  of  the 
edges,  equations  (85)  are  substituted  in  equations  (77)  and  (78).  The 
resulting  equations  may  be  written  in  matrix  notation  as  follows: 


a 

Rh 


HftH  M + h N 

where  the  matrix  JA^  is  given  in  Table  32  partitioned  so  that 


3 
5 

l1^ 

^9 

^l] 


A 


S 


T.115376 
-.012713 
.00420168 
-.00194778 
-.0010647 
-.00064353 
-.00041703 
-.120621 
-.0123367 
-.00395848 
-.0018093 
-.000979049 
-.000586332 
-.000376233 
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K. 


lll 


13 


(86) 


and  C,  (i  ■ 1 to  14)  are  non-linear  functions  of  v _ as  listed  in 
l m,n 

Table  33*  Table  33  should  be  read  in  such  a manner  that  * 1*370 
wl,l3  + *29B  wl,l^wl,3  “ #069  wl,l  W3»l  + * • • 

With  the  non-linear  terms  not  involving  w^  w^  y w^  y w^  y 

W1  5*  w5  1 neSlected  the  problem  can  be  solved  as  follows*  First,  we 

note  that  equations  (86)  are  linear  equations  in  terms  of  the  unknowns 
Kj . It  is  therefore  possible  to  solve  these  equations  for  Kj  in  terms 

of  the  non-linear  functions  of  w„  In  matrix  form,  we  have 

m,n  w 


M - — (V 

L il  Rh  L ^ 

jl  y h l 

(87) 


where  A 


-1 


ij 


denotes  the  inverse  of  the  matrix  A. 


-1 


Table  34  shown  partitioned  according  to  A * 


P Q 
R S 


ij 


and  is  given  in 
The  inversion  of 


of  the  matrix  A^j  was  carried  out  by  the  use  of  high  speed  electronic 

calculating  machines.  Since  we  have  neglected  all  non-linear  terms  of 
wm,n  those  living  v^,  v3(1,  v3(3,  v15,  w^  with  Kj 

found  as  functions  of  these  deflection  coefficients,  we  can  substitute 
these  Kj  into  the  first  six  equations  of  (85)  and  obtain  a system  of 

six  non-linear  equations  in  terms  of  these  deflection  coefficients.  In 
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matrix  notations,  the  first  six  equations  of  (85)  may  be  written  as 
follows: 


Wl.l 

wl,3 

w3,l 

w3,3 

Wl,5 

W5,l 


E. 


.089783 

.0049583 

.0056107 

.0011085 

.0010103 

.0012757 


and 
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o o o o o o 
o o o o o o 
o o o o o o 
O O O o o o 


o o 


o 


o 


«0 


o 


Oi  to 

O *A 

fc  £ 

O 'O  O (V  o o 

8 8 


$ 

£>  o 

& 


*A 

c$ 

IA 

*A  o 

8 


o 


o o o o o o 
o o o o o o 
o o o o o o 
o o o o o o 


o o o o o 


o 


CA  0s 

$ £ 

*A  o 

(V  -st 

° I § ° 

• • 


o 


• • 
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By  substituting  (87)  inte  (88),  we  obtain 


To  simplify  the  notation,  let 
2 


a 

Rh 


and 


b M MfJ  -bN-N 


then 


1 

h 


N ■;[■.]•  H 


(89) 


The  above  matrix  operations  have  been  carried  out  and  we  obtain  the 
following  results: 


" .0242193  ~ 

H1 

-.00335124 

H2 

.00274422 

H3 

.00016766 

and  Hj  * 

H4 

-.00120163 

H5 

.00120163 

_H6_ 

(90) 


where  the  are  listed  in  Table  35. 

The  six  non-linear  equations  (88)  can  be  solved  by  an  iteration 
method  as  follows*  Consider  at  first  the  first  equation.  Assume  that 
all  deflection  coefficients  except  w^  ^ are  zero*  Then  we  have  a cubic 

equation  in  Wj  , . This  equation  can  be  solved  by  trial  and  error  which 

give 8 a first  approximation  to  w,  Next,  consider  the  first  two 
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equations  and  assume  that  all  deflection  coefficients  except  w^  and 
v,  - are  zero*  Substituting  the  value  of  w,  , obtained  previously 

If?  ±9‘L 

into  the  second  equation*  a first  approximation  to  the  value  of  vi,3 

can  be  obtained*  Next  substitute  this  value  of  into  the  first 

equation  and  calculate  the  second  approximation  to  the  value  of  w^^* 

This  process  is  repeated  until  the  desired  accuracy  is  obtained*  It 
may  be  pointed  out  that  after  the  second  cycle,  it  is  not  necessary 
to  compute  the  w^  ^ by  using  the  value  of  w^  ^ computed  In  the  previous 

cycle  and  a value  which  appears  to  give  a more  rapid  convergence  can  be 
used*  In  this  way  the  rapidity  of  convergence  can  be  accelerated*  §jr 
taking  more  equations  and  solving  for  more  unknowns,  all  six  unknown 
deflection  coefficients  can  finally  be  solved* 


It  was  found  that  the  convergence  of  the  solution  is  fairly  rapid 
if  the  iteration  of  the  first  equation  is  carried  out  according  to  the 
scheme 


(n+l) 


where 
the 


2iL 

h 

(»>. 


V' 


— (.0242193 

2 

a 

% 

iw. ; M 1/3 

.620 

Rh 

h 

1 J 

(91) 

(a)  w (n) 

U 

(n) 

(n)  (n) 

w (n) 

L,1  »W1,3 

*3,1 

,w3,3  ,w1>5  , 

*w5,l 

) with 


(n)«  term  deleted;  and  the  remaining  five  equations  are  iterated 
^vl*l  j according  to  the  scheme: 


w (n+1>  2 
wi  _ a . 


where  H 


<n) 


H. 


h 

(n) 


Rh 


M 


+ H 


(n) 


(92) 


(v. 


(n)  « (»)  „ (»)  „ (») 


1,1  *1,3  *3,1  *1,5  *5,1 


(n) 


). 


In  the  above  equations,  the  superscript  (n+l)  indicates  the  (n+l)-th 
approximation  and  n indicates  the  value  assumed  before  the  (n+l)-th 
cycle  is  computed.  Since  these  equations  are  cubic  equations,  it  is 
possible  that  these  equations  may  have  more  than  one  set  of  solutions* 

To  show  the  uniqueness  of  the  solution  is  a very  difficult  task  although 
it  is  intuitively  clear  that  there  should  be  only  one  equilibrium 
configuration  in  this  problem.  This  is  partially  bora  out  ty  solving 
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the  first  equation  when  only  w^  ^ is  assumed  to  be  different  from  zero* 

It  was  found  that,  in  this  case,  we  hare  one  real  root  and  two  complex 
ones* 

The  iteration  of  the  six  equations  according  to  the  scheme  of 
Eqs.  (91)  and  Eq*  (92)  are  shown  graphically  in  Fig*  29*  In  this  figure 
the  solid  lines  lead  from  the  assumed  vilues,  w^(n)  to  the  iterated 

values,  w1^n+1^  while  the  broken  lines  show  the  next  guess  made  from  the 

results  of  the  previous  iterated  value  which  is  used  in  the  next  iteration* 

The  iteration  converges  to  the  values  for  w^  with  g 

— - 337*5  (a  - 9% 

Rh 

R * 6",  and  h ■ *040")  as  follows: 

■la?!-*,  1*824,  "1*2  m 0*116,  -2*1  ■ *527,  "2*2  ■ .0194,  -1*2  > -,071,-2*1b  .2025 
, h h h h > h h 


from 


These  values  for  Wj  are  used  in  computing  Ci  from  Table  33  and 
Table  31.  We  obtain  therefore. 


'8 


u10 

°u 

°X2 

°13 

cu 

and 


15.081T 
.4988h3 
*039917h3 
-.0002530b3 
*00001019h3 
.00001045b3 
-.0000002114b3 
20.33h3 
-1.07 5h3  ^ 
-*07069h3 
*009659h3 
*000113 5h3 
-*00003037h3 

m 

-#00000006266h; 


1.840 

1.830 

1.820 

-.120 

-.110 

-.100 


.530 

.520 

.210 

.200 

.030 

.020 

.010 


L 

v- 

V 
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A 
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Iteration  of  the  Six  Non-Linear 
Equations 


Fig*  29 
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9 Pm  9 I"m  9 Pm  9 Pm  9 I' 


Dx  * I6.9219h3 
D2  = ,90l8h3 

D3  = ,1300b3 

D4  - ,001054h3 

Dc  - -,00000934h3 
D6  ■ -,00000687h3 

D?  = -,00000024h3 
Dg  « -,00000002h3 

D9  = -,5360h3 
D1q  * -.129Ch3 
D11  = -•01286h3 
D12  * -,0004967h3 
D13  «=  -.00002564b3 
D14  * -. 000003 39h3 


D1C  » .OOOOOOllh3 
JL5 

* -,04668h3 
D1?  • -.006l35h3 
D18  “ -•00°57689h3 
D19  - -,00008891b3 
D20  ■ ,00001768b3 

Da  » ,00000798h3 

D22  « ,00000001b3 

« ,00008837h3 

D24  “ • 001891b3 

» ,0006802b3 

Dgg  « ,00007449h3 

- . 0000020 8h3 

D28  * -*00000251h3 


Using  these  values  for  and  D^,  the  K^'s  can  be  computed 
according  to  Eq,  (87)  by  which  we  obtain 


K1  » - 3032 


— K,,  * - 27664 
Rh 


Rh 


T?  » 2411 


K3  » -7765 


Rh 


ILL3  - - 32664 


3183 


Rh 


K5  - - 12703 

Rh 

- 278.3 

af 

Rh 

Tu  « 3817 

Krj  - - 16553 

A2. 

Rh  * 

» 924.4 

af 

Rh 

T13  » 4520 

^ - 22669 

— f5 

Rh  ? 

« 1694 
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t: 


Substituting  these  results  into  equations  (85),  the  values  of 

w„  _ are  as  follows: 
m,n 

"i.i^ 

« 1.77 

w9>1A  - .059 

w13,3A 

» -.0015 

vi,3/h 

» -.110 

wll,lA  “ *°32 

w3,5A 

« -.030 

V3,1A 

- .521 

w13,1A“  •08° 

w5,5A 

- -.0076 

- .018 

a -.014 

w7,5A 

» .034 

W1,5A 

= -.099 

j/h  ■ .018 

u9,5A 

a -.00076 

w5,1A 

= .201 

w9  ^A  » -.0046 

wU,5A 

» -.00028 

v7,1A 

* .133 

wn>3A  * -.0025 

w13»5A 

« -.000097 

w1,7A 

» -.060 

w7>9A  * *038 

W11,UA 

» -.0020 

w3,7A 

» -.027 

w^^  * .0046 

w13.11A 

= -.0014 

v5,7A 

= -.010 

w11,9A  ’ •°033 

W1,13A 

» -.011 

w,(7A 

■ .032 

h13>7A  “ .0024 

w3,13A 

a -.0078 

w9>7A 

• -.0018 

U1,11A  * -*017 

W5,13A 

- -.0056 

w11,7A 

« -.00094 

w3,11A  * ~'012 

w7,13A 

» .018 

v13,7A 

« -.00054 

V5tu/h  * “•0075 

u9,13A 

« -.0028 

- -.022 

w7,1XA  - •°22 

V11,13A 

» -.0019 

',3,9A 

T 

- -.00085 

w9>11A  - -.0030 

v13,13A 

» -.0014 

w5(9A  - .0059 
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Comparing  the  Tallies  of  these  deflection  coefficients,  it  is 
observed  that  it  is  reasonably  accurate  to  neglect  all  the  non-linear 
terms  in  these  equations  except  the  first  six*  The  derivation  of  the 
values  of  the  first  six  deflection  coefficients  so  computed  from  the 
results  obtained  earlier  is  probably  due  to  the  accumulated  errors  in 
the  substitution,  especially  since  the  matrix  inversion  . -1  was 
found  to  no  more  than  four  significant  figures* 


The  stress  function  coefficients  bp  can  be  computed  by  substi- 
tuting the  values  of  these  deflection  coefficients  into  Table  30  from 


which  ve  obtain 

V * •569Eh2 

b4,0 

■ 

*0048Eh2 

b0,4  * *0072Eh2 

b4,2 

6 

*0041Eh2 

b0,6  * -*0013Eh2 

b4»4 

m 

-*0017Eh2 

bQ  g - -*0006SEh2 

b4*6 

m 

-*0025Eh2 

bo,io  " •00027Eh2 

b4,8 

m 

*00034Eh2 

b0  12  " -»00006Eh2 

b6,0 

B 

*00l6Eh2 

b2,0  * *020Eh2 

V 

V 

-*00084Eh2 

b2,2  * *079Eh2 

b6,4 

m 

-*00098Eh2 

bn  . - *00Q4Eh2 

b6,6 

B 

*00022Eh2 

b2,6  * -*0019Eh2 

b8,0 

m 

.OOOlSEh2 

b2,8  * -•0012Eh2 

b8,2 

m 

-*00002Eh2 

b2  10  “ *00005Eh2 

b8,4 

B 

-. 000009Jh2 

b2,12  = -O0^2 

^0,0 

B 

*00002Eh2 

With  the  values 

of  w and  b determined,  we  thus  obtain  tht 

m»n  Pt9 

solution  to  our  problem  of  clanping  a curved  plate  in  a flat  frame* 
The  actual  distribution  of  the  deflection  and  stresses  can  be  found 
from  the  relations  previously  listed  and  shown  below: 


v - 2 2 v sin  — sin  222 

* a b 


(70)  bis 
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Unit  M&mbrane  Normal  Forces  tbit  Msmbrane  Normal  Forces 

(at  x « a/4  and  y * b/4)  (at  x ■ a/2  and  y = b/2) 

Fig.  34  Fig*  35 


Tbit  Membrane  Shear  Forces 
(at  x « a/4) 

Fig.  36 
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and 


a 

.6 

xO 

d 

H 

V 

dx 

» 92*1 

*xyO 

0x9y 

where 

F * £ L b cos  222  cos  222  (71)  bis 

P*^  A Tl 


The  variation  of  the  non-dimensional  unit  membrane  forces  and 

deflections  are  computed  using  the  finite  number  of  terms  obtained  above. 

They  are  plotted  along  four  sections  of  the  plate,  that  is  along 

x * a/4  and  a/2  and  along  y **  b/2  and  b/4.  The  values  are  plotted  for 

only  one-quarter  of  the  plate  (0-x—  a/2,  0-y— b/2)  in  Figs,  30-36 

because  of  the  symmetry  of  a , a and  w.  The  computed  non-dimensionalized 

x y 

forces  are  represented  by  the  solid  lines  of  Figs,  34-36.  In  Figs,  30- 
33  the  deflection  curves  are  shown  magnified  for  the  portions  of  the 
plate  near  the  boundaries.  It  can  be  seen  from  these  magnified  views 
that  the  boundary  conditions  on  the  deflections  which  theoretically 
require  zero  slope  and  deflection,  seem  to  be  well  approximated  by  taking 
the  finite  number  of  terms. 


SECTION  XIII 

DETERMINATION  OF  THE  BUCKLING  LOAD  FOR  A PRESTRESSED  PLATE 

Having  determined  the  inplane  forces  in  the  prestressed  plate,  the 
next  problem  is  to  find  the  buckling  load.  Here  the  use  of  the  term 
buckling  load  can  be  subjected  to  objection,  because  the  prestressed 
plate  will  have  initial  deflections  and  there  will  not  be  a definite 
load  at  which  the  "flat”  form  of  plate  becomes  unstable.  However,  if 
we  assume  that  these  initial  deflections  do  not  exist,  then  we  may 
determine  the  "buckling  load"  by  the  Rayleigh-Ritz  method. 

Let  w be  the  lateral  deflection  of  the  plate  due  to  buckling* 

Then  the  bending  strain  energy  in  the  plate  with  all  edges  clamped  when 
buckling  occurs  is 


WADC  TR  54-8 


64c 


IT 


2 

+ ^M()  dx  dy 

07 


(93) 


where  D is  the  flexural  rigidity  of  the  plate.  The  work  done  by  the 
inplane  forces  is 


W » 


1 

2 


+ 2N 


xy 


9w  6w  | 
dx  dyj 


dx  dy 


(94) 


where  Nx  = Nx0,  Ny  = Nycr  + NyQ,  the  subscript  0 refers 


to  the  original  inplane  forces  and  the  subscript  cr  refers  to  the 
buckling  value.  At  buckling,  we  have 


U - W « 0 

from  which  we  obtain 


(95) 


m „ SLz-g 

Nycr  T (96) 

where  ^ 

T * — jj  (— ) dx  dy 
2 yjk  0y 


and 


W « W - N 


ycr 


T 


It  must  be  emphasized  that  the  above  relations  are  correct  only 
if  there  were  no  initial  deflections.  With  the  presence  of  the  initial 
deflections,  the  values  of  the  inplane  forces  originally  due  to  the 
prestressing  will  be  changed  when  the  load  Nyer  is  applied  and  therefore 

the  above  method  of  calculation  may  induce  an  error  which  is  not  small. 
However,  in  the  absence  of  a better  method,  it  was  decided  to  try  to 
use  the  Rayleigh-Ritz  method. 

The  Rayleigh-Ritz  method  can  be  carried  out  as  follows.  First, 
the  deflection  w is  assumed  in  the  form  of  a series,  every  term  of  which 
satisfies  the  boundary  conditions  but  with  undetermined  parameters. 
Substitute  the  series  into  equation  (95)  and  carry  out  the  integration. 
Then  minimize  the  resulting  expression  with  respect  to  the  undetermined 
parameters.  This  gives  us  a system  of  homogeneous  equations  involving 
the  buckling  load  A non-trivial  solution  is  obtained  if  we  equate 
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the  determinant  of  the  coefficients  of  these  parameters  to  zero,  from  which 

N is  calculated* 

ycr 

The  boundary  conditions  for  a clamped  plate  are  that  the  deflection 
and  the  slope  normal  to  the  edge  are  zero  at  the  edge,  i*e* 


at  x 

* 0,a 

„ _ dw 

v * 

dx 

at  y 

- 0,b 

w.  2s 
5y 

To  satisfy  these  conditions,  we  may  assume  v in  the  following  form. 


w ■ sin  -=  sin 


25  a4n 


Z S w sin 


225  SEE 


b jbfI  n*l 


where  w _ are  undetermined  coefficients* 

nifiQ 

Substituting  (97)  in  (93),  we  find 

b,  a , 


0 0 


Us  2^  s i«E  pwl)2  C0B 


(n-1)2  cos  1 ,JH  8to 

b J a , 


a oinSOE 
b b 


£(n+l)^  cos  te&teX  - (n-1)^  cos  ^ x ^ 
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Performing  the  indicated  integration,  we  obtain 


D7T4  f 


2 f ^m,!2  (»Wn) 

V.  HFl  7 a 


128a' 


+ 2 2v,  2 (-)  (n4+6n2+l) 

n=l  i,n  b 


- £ Wl*n-2Wl»n  2 [(n+l)2*8»] 

- J wm-2,lwm,l  2 ^l)2-8n] 


+ 2 2 4v 

df!  n=*l 


sf/bx 

'] 


) (m4 +6m2 +!)+(—)  (n4+6n2+l) 


+ 2 (-)(m2+l)(n2+l) 


- 2 2 2w  _w 

m=l  n=3  m*n“2  ffi»n 


(“)  (n4-6n2+8n-3) 
b 


+ (— ) (m4+6m2+l)+2(-)  (m2+l)  (n-1)' 
a b 


2 2 2w  _ w 

np 3 n=l  m“2»n  m»n 


(-)  (m4-6m2+8m-3) 
a 


+ (S)  (n4+6n2+l)+2(-) (m-l)2(n2+l) 
b b 
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+ 2 2 

mF3  n«3 


wm-2,n-2  wm,n  j^b 


) (n^-6n  +8n-3) 


+ (i)(m4-6m2+8m-3)+2(-)(m-l)2(n-l)^ 


- 2 2 2w 


m»l  n*=l 


%n+2 


v F (- 

m,n  [ ^ 


(-)  (n4-6n2-8n-3) 


+ (-) (m4+6m2+l)+2(-) (m2 +1) (n+l)2] 
a b J 


2 2 2w 

m=l  n=l 


r,b 

r W | (— 

m+2,n  m,n  I N& 


')  (n»4-6m2  -8m-3) 


+ (-)  3 (n4+6n2+l ) +2  (-)  ( m+l)  2 (n2+l )/ 
b b -J 

* ”,  = “*-2,0*2  “m,n  J&  (■‘•AfWj) 

m“3  n=l  * La 


+ (S)3(n4-6n2-8n-3)+2(-)(m-l)2(n+l)2l 


_J 


0O  oo 

+ 22 

m=l  n=3 


Ha 

Wm+2,n-2  Wm,n  / ^ 


(— ) (n4-6n2+8n-3) 


+ (-)  (m4-6m2 -8m-3)+2(-)  (m+l)2(n-l)2  I 

a b J 
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+ 2 

m3 


2 wm4.9  n+o  wm  n / 1 (*)  (n^-6n2-8n-3) 
nc3  m+2,n+2  m,n  ^ b 


+ (-)  ( m^-6m2  -8m-3 ) +2  (-)  ( m+1 ) 2 (n+1 ) ' 
a b 


(99) 


It  may  be  mentioned  that  we  obtain  the  above  general  form  for  U if  all 
terms  in  the  assumed  series  for  w are  taken.  If  only  a finite  number 
of  particular  terms  are  taken,  the  coefficients  of  those  terms  omitted 
naturally  should  be  deleted  from  the  above  expressions. 


Substituting  the  assumed  form  for  w , into  the  expression  for  T 
and  integrating,  we  obtain 


T ■ - 


£a_ 

128b 


[2 


2(l+n)(2w„  „ -v  0 

' ' m,n  m-2,n 


w0  _)+(n  -1) (-2w  0 -2w_  _ 

2—uif  n 2+m,n  id^ 


2wm,2-n  + wm-2,n-2  + Wm-2,n+2 


- w 


m-2,2-n 


~ w2-m,n— 2*^2— m,n+2+w2— m,2-n+w2+m,n— 2+w2+m,2+n”^2+m,2-n 


3 


(100) 

The  work  expression  due  to  the  prestresses  existing  in  the  plate 
will  now  be  computed  assuming  that  these  stresses  remain  the  same  as 
before  the  buckling  load  is  applied  and  remain  constant  during  buckling. 
Thus 


8 

0x8y 


(101) 
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and 


OO  0° 


F » Z Zb  cos  222  cos  322 
p*0  q«=0  p,q  a b 


To  simplify  the  writing,  let 


1 1 


NxQ  (&)  dx  dy 


ax 


(— ) dx  dy 
ay 


V = + 
xy 


„ 3w  9w  , 

N — *-  dx 


xyO 


dx  dy 


dy 


(102) 


Then,  after  substitution  and  integration,  we  find 

hrA 


V. 


x 512a2  m=l 


Z Z Z Z v v,  . - £(n-t)2b  . -2(n+t)2b  . 

=1  n-1  k-1  t=l  m»n  k>t  b L s*n-t  s»n+t 


',ww)W(,,w)\^  [2-  £ %] 

f(2mk+2)  +(2mk+2)  +(2ri^2) 

IS  s s 
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+ (-mk+l-2m)  +(-mk+l-2m)  6^’Z 

+(-jnk+l+2m)  £m"lc+2  +(_nk+l+2m) 

s s 


■K-»k-l+2m)  ^-“-k  +(^nk-l+2n)  i“tk-2+(-mk-l^m)  tf***2 


where  £ j ■ 1 when  i = j and  & ^ = 0 when  i / j 

t__4  oo  OO  k'*  oo  p* 

2 2 2 2 2 v w.  (S)  j2(m-k)  d . 

2 m=l  n=l  k**!  t=l  m»n  *»*  b «-  **.' 


(103) 


V 


512a 


'(’^'■2)2bB-k-2,B'(2^Mk)2b2-1»k,s-(®-k+2)2bm-i*2,! 

*(*4Al-2,,J(",k)2lrt,,*(“W)2‘„k(2)t] 

I2-?  %]  f(2nt+2)  ^_t+(2nt*2)  ^~n+(2nt-2)  4”+t 


+(-nt+l-2n)  ^g”n~t+ ( -nt+l-2n ) ^”2+n~t+(-nt+l+2n)  &2+n-t 


+ ( -nt+l+2n ) ^t“2-n+(-nt-l+2n)  &2-n-t 

B S 

+(-nt-l+2n)  ^J2+n+t+(-nt-l-2n)  52+n+^j 


(104) 
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h / 

256a2 


00  00  00*0 
2 2 2 2 
df*1  n»l  k»l  t=l 


wm,nVk,t 


(S) 

b 


|[(2^)b2^>e*(2^k)b2_wk>s-(2*»k)b2+wt(11 

-(”h*-2)bm^-2,0’(”+lt-2)b1»k-2,S]  [2n(t'n)  ^ 

-5n(n-t)  ^+2n(n+t)  ^+^+(n+l)(24u-t)  ^ 

S 8 3 

♦ 6.-n) (2-n+t)  i2’n+t-(l^)(n+t+2)  $£+t+2+(n-l)  (n-t^)^^"2 

' 8.8  S 

-(l+n) (t-n-2)  ig”n“2+^.-n)(n+t-2)  ^+t"2] 

+ r^(^)bm^>s+2(^k)w,8+(2+^k)b2+BHk,3 

-(2^k)b2^k>s-(2+m+k)b2+^k>s+(m-k^)b^2>8 

-(n^k^b^^JJ  pl-t)(2-m-t)  ^^md+t)  (2+t-n)S2+t_n 
-m(l+t)  (n+t+2)  f*+t+2-m(l+t)  (n-t-2)  S1~t_2+m(t~l)  (t-n-2)  ^tHa”2 

8 S S 

+»(l-t)(„+t-2)?U-2]}  (2-  $ S”) 

J (105) 
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To  obtain  numerical  answers,  let  us  consider  a plate  with  the 
following  dimension  ratios: 


b/a  - 1.5,  — - 337.5 
Rh 


(106) 


These  ratios  are  those  corresponding  to  the  plate  tested  having  a 
thickness  of  .040  in.  and  an  initial  radius  of  curvature  of  6 in. 


Using  eqs.  (99),  (100),  (103),  (104),  and  (105),  and  taking  only 
one  term  in  assumed  form  of  deflection,  namely. 


we  find 


w * w,  « sin2  — sin  ^ sin 

1,2  a b b 


TT  , o*-,  iA>  2 

U - 1.361  — j vl  2 

a 9 


(107) 


T - -.156  ir2  w,  . 


- -.28;  A 


a2 

wl,2 

TT^D 

v 2 

2 

a 

Wl,2 

iAd 

2 

2 

a 

Wl,2 

(108) 


Substituting  the  expressions  (107)  into  the  buckling  criteria  yields 
the  following  equations: 


2(1.563  + .156  7 r^)  2 « 0 

a 9 


2 

or  N_  « 10.02 
7 bl 


(109) 
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The  buckling  load  for  a non-prestressed  plate  of  the  same  dimensions 

can  be  found  in  an  analogous  manner  by  equating  to  zero  the  workjlone  by 

in-plane  forces  due  to  prestressing,  i.e.,  omitting  y f y and  W • 

x y xy 

The  buckling  criteria  for  the  corresponding  non-prestressed  flat 
plate  yields 

.156  7T^Ny)  v1>2  * 0 


(HO) 


2(1.361 


or  N. 


■8.72  & 


From  the  above  calculations,  we  find  that,-  by  assuming  the  deflection 
in  the  form  of  (107)  and  assuming  that  the  in-plane  forces  remain  constant 
when  the  buckling  load  is  applied,  the  Rayleigh-Ritz  analysis  shows  that 
the  buckling  load  is  increased  14.9$  by  prestressing.  It  may  be  mentioned 
that  a more  accurate  value  for  the  buckling  load  of  a non-prestressed 
clamped  rectangular  plate  with  b/a  “ 3/2  is  (12) 


N 

ycr 


a 


which  is  4.7$  lower  than  the  value  just  computed. 


(in) 


To  stttfy  the  effect  of  taking  more  terms  in  the  assumed  deflection 
series  (97),  we  consider  next  the  following  deflection  form 


sin^S 


(112) 


Substituting  this  deflection  form  into  (99),  (100),  (103),  (104)  and  (105), 
we  obtain 

0 * I®7*11  m1,22-137*63u1,2W3,2+565'63w3,22J> 


T " ' ^ f1°“l.22-6-67“l . 2“3 ,2*6‘67W3 ,tj 
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Wx  ■ “ ^2  |^8#20w1,2  “31  * 3 5wl , 2w3 , 2+85 *40w3 , 2 ^ 

V--f^2  f3-13vl,22+7*88wl,2V3,2-2-70w3,22J 

W3cy  “ " ^~2  ^~2#27w1,2  +2*97w1,2w3,2"*11w3,2  ^ 


(U3) 

First,  let  us  calculate  the  buckling  load  for  a non-pre stressed 
plate.  Let 


K « N a2/ir2D 
ycr 

then,  we  have 


(114) 


U-N  T 
ycr 


64a2 


^(87.11+10K)w1>22-(137.63+6.67K)w1>2w3^2+(565,63+6,67K)w3>22J 


(115) 

Minimizing  the  above  expression  with  respect  to  the  undetermined 
coefficients  v3  2 and  2 we  obtain  the  following  two  equations  at  buckling, 

/(174*22+20K)w1  0-(137.63+6.67K)w_  .1-0 
dwi,2  64a"1  l 1,2  3»2J 


^ f -(137.63+6,67K)w1  2-(1131.26+13.33K)w3  27  - 0 
9w3  2 64a  n f 


(116) 
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Setting  the  determinant  of  the  coefficients  of  these  two  equations 
equal  to  zero,  the  conditions  for  a non-trivial  solution  give 


K2  + 104.06K  + 802.10  - 0 
from  which  we  find 


(117) 


or 


K - - 8.38,  or  - 95.7 


H - 8.38 

yer  2 

a 


(118) 


Although  this  answer  was  arrived  at  using  only  two  terms  of  the 
series,  it  is  only. 0656  higher  than  the  more  accurate  value  found  by 
Levy  (12)  and  listed  as  Eq.  (ill).  For  the  prestressed  plate,  we  have 


7-T-To  -To  -To  -TS 

x y xy  64a 


|(99.9+10K)w1>22-(158.1+6.67K)w1>2w3>2+(648.21+ 


6.67K)w. 


4 


(119) 


Minimizing  this  expression  with  respect  to  w.  „ and  w_  we  obtain  the 

J-t*  i,z 

following  two  equations 

(199.8*2CK)w1>2  - (158.12+6.67K)w3>2  » 0 


-(158.12+6.67K)w1>2  + (1296.4+13.34K)w3>2  * 0 (120) 
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) 


The  condition  for  a non-trivial  solution  is 

K2  + 119. IK  + 1052.7  = 0 
from  which  we  find 


(121) 


or 


K ■ - 9.61,  or  - 109.52 


N * 9.61  2-2 

ycr  , 


(122) 


Although  the  actual  values  of  w^  ^ or  2 canno^  b®  ^oun^  there 
ratio  can  be  calculated  using  either  of  Eqs.  (119). 

Hie  first  of  Eqs.  (119)  may  be  rewritten  as 

wi  0 

(199.8+20K)  - 158.12  + 6.67K  = 0 

W3,2 


or 


Substituting /K  * - 9.61  the  equation  becomes 


7.6 


1*2 


*3,2 


"1.2 


3,2 


- 222.22  = 0 


29.2 


(123) 


(124) 


for  the  prestresssed  plate. 

The  buckling  loads  and  ratio  of  deflection  coefficients  were  computed 
for  both  the  non-prestre3sed  and  prestressed  plates  using  other  coefficients 
and  combinations  of  coefficients  in  the  assumed  deflection  series  (109). 

The  results  are  tabulated  in  Tables  36  and  37  and  shown  graphically  in 
Fig.  37  for  various  percentages  of  the  calculated  amounts  of  pre stressing. 
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Note  that  due  to  orthogonality  relations  between  the  assumed 
deflection  form  and  the  stress  function  series  talcing  a limited  number 
of  terms  of  the  deflection  form  (wm  ) will  automatically  limit  the  number 

of  terms  describing  the  stress  function  (bp  g).  The  stresses  computed 

with  these  limited  number  of  terms  (due  to  considering  only  one  term, 
v1t2  '^n  deflection  series)  are  compared  with  those  computed  with 

a larger  number  of  terms  in  Fig.  34t  35,  and  36.  In  these  figures  the 
stresses  corresponding  to  the  limited  number  of  terms  are  shown  with 
broken  lines  while  the  more  exact  representations  are  shown  with  solid 
lines. 


Results  of  Rayleigh  Ritz  Analysis 
Fig.  37 
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SECTION  XIV 


. 5 

Q ^ 


ri,r?<T'r?Ciir?r?°° 


v ntn'rn 

iaiaxaiaiaiaiaiaia 


? 7 


IA  IA 

<?  <?  V 


rjl  O 


IA  IA 


V T 


ri,ti,7'f'?^cr'»‘V 

t^nn^stsfsf^sf 


sitvm?? 


ri,ciif?Yl?v?crt,?cr 

(VCV(VtV(V(V(VCV(V 


4A  lA  IA 

1A  IA  • • • 

.*  .*  1 1 7 1 7 


'j  n t ■?  i 'i  n 

HHHHHHHHH 
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These  results  were  duplicated  within  - .001  inches  for  similar  tests 


TABIE  1 - lb 


(Flatness  tests  on 

•032  inch  Aluminum  Plate) 

Position 

Deflection 

-3 

x 10  inches 

Deflection 

Position  x 10^  Inches 

Position 

Deflection 

-3 

x 10  inches 

1-1 

0 

2-1 

-18 

3-1 

-12 

1-2 

-15 

2-2 

-38 

3-2 

-27 

1-3 

-23 

2-3 

- 

3-3 

-37 

1-4 

-27 

2*4- 

-49 

3-4 

-30 

1-5 

-25 

2-5 

—46 

3-5 

-28 

1-6 

-18 

2-6 

-35 

3-6 

-20 

(minus  indicates  downward  deflection) 


TABLE  1 - lc 


(Flatness  tests  on 

.040  inch  Aluminum  Plate) 

Deflection 

Deflection 

Deflection 

Position 

-3 

x 10  inches 

-3 

Position  x 10  inches 

Position 

-3 

x 10  Inches 

1-1 

0 

2-1 

0 

3-1 

+3 

1-2 

-6 

2-2 

-9 

3-2 

-1 

1-3 

-14 

2-3 

-19 

3-3 

-10 

1-4 

-12 

2-4 

-17 

3-4 

-13 

1-5 

—8 

2-5 

-11 

3-5 

-13 

1-6 

-5 

2-6 

-5 

3-6 

-6 
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Note:  Gage  Numbers  refer  to  those  indicated  in  Fig, 


TABLE  3 


Results  of  Tests  on  24  S-T  Al.  Alloy  Plate 
Radius  of  Curvature  - 9 l/4  inches 


e (10^)  In. /in. 
61 


Gage 

Test  1 

Test  2 

Test  3 

Test  3 

(using 

corrections) 

1,2 

-43.5 

-26 

-24 

(1.10) 

+42 

3,4 

+61.5 

+71.5 

+71.0 

+71 

5,6 

+75.0 

+81.5 

+84.5 

+84.5 

7,8 

-140.5 

-98 

-in.  5 

(8,9) 

+76.5 

7-2, 8-2 

X 

X 

X 

+13.5 

Note:  + = tensile 

- = compressive 
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tabie  a 

Test  Data  for  Determining  Initial  Stresses 
in  ,040  Aluminum  Plate 
Initial  Radius  « 8,70  inches 

®a  = membrane  strain  (micro-inches  per  inch)  see  Fig,  4 for  placement  of  gages 
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-76  -85  -47  1/2  -26  -24  +9 

-75  1/2  -84  -49  -29  -34  1/2  -11  l/2 
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-50  1/2  -30  1/2  -14  1/2  -7  1/2 

-62  1/2  -43  -25  l/2  -17 


TABLE  7 (continued 
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TABLE  10  (continued) 
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TABLE  11 


Membrane  strains  for  second  *040  Aluminum  Plate 

Gage  Position  Membrane  strain 

see  Fig»  5 micro-inches  per  inch 


A-l 

-10 

A-2 

+65 

A-3 

-10 

k-U 

-15 

B-l 

+10 

B-2 

+10 

B-3 

0 

B-4 

-60 
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TABLE  12 


Membrane  strains  for  first  ,040  Aluminum  Plate 


Gage  Positions 

Membrane  strain 

Assumed 

(see  Fig,  5) 

(see  Fig,  4) 

micro-inches  per  inch 

Value 

C-l 

A-l 

-16 

-26 

-21 

C-2 

A-2 

-10 

-10 

-9 

C-3 

A-3 

+14 

+12 

. * 

+10 

A-n4 

+21 

+19  1/2 

+20 

B-7 

B-l 

-16 

-18 

-21 

B-6 

B-2 

-21 

-23 

-26 

B-5 

B-3 

+43  1/2 
+38  1/2 

+41 

B-4b 

B-4 

-18  1/2 
-19 

-19 
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TABLE  13 

Membrane  strains  for  *032  Aluminum  Plate 

Membrane  strain 


Gage  Positions  micro-inches  per  inch  Assumed 


(see  Fig.  5) 

(see  Fig.  4) 

Test  1 

KHI 

Test  3 

C-l 

A-l 

-130 

-131 

-52  1/2 
-54  1/2 

-40 

-26  1/2 

-70 

C-2 

A-2 

+10  1/2 
+18  1/2 

+25  1/2 
+25 

+35  1/2 
+46 

+30 

C-3 

A-3 

+37  1/2 
+43 

+15 

+21  1/2 

+30  1/2 
+33  1/2 

+30 

C-4 

A-4 

+10 

+10  1/2 

-22  1/2 
-14 

-3 

-3  1/2 

-5 

B-7 

B-l 

-189 

-176  1/2 

-140 
-135  1/2 

-36 

-32  1/2 

-120 

B-6 

B— 2 

+38  1/2 
+35  1/2 

+51 

+49  1/2 

+53  1/2 
+54 

+45 

B-5 

B-3 

-45 

-46 

-42  1/2 
-43 

-28  1/2 
-29 

—40 

B-4b 

B— 4 

+24  1/2 
+24  1/2 

+15  1/2 
+14 

+27  1/2 
+28 

+20 
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TABLE  14 


Membrane  strains  for  .051  Aluminum  Plate 


Gage  Positions 
see  Fig.  5) (see  Pig.  4) 

Membrane 
micro-inches 
Test  1 

strain 
per  inch 
Test  2 

Assumed 

Value 

C-l 

A-l 

-42  1/2 

-76 

-60 

-37 

-75  1/2 

C-2 

A-2 

-44 

-85 

-60 

-35  1/2 

-84 

C-3 

A-3 

-32  1/2 

-47  1/2 

-35 

-21 

-49 

C-4 

-25  1/2 

-26 

-25 

-18  1/2 

-29 

B-7 

B-l 

+43  1/2 

+56 

+50 

+44 

+51 

B-6 

B-2 

+68 

+76  1/2 

+70 

+69  1/2 

+71 

B-5 

B-3 

-14  1/2 

+ 1/2 

-5 

-2  V2 

+2  1/2 

B-4b 

B-4 

-54  1/2 

-41  1/2 

-45 

-42 

-44  1/2 
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TABLE  15 


Calibration  Tests 
Performed  on  Flat  Plates 


Plate  Ho.* 
F, 40, 1.5, 30 


F, 40, 1.5, 31 


F,40,1.5,32 


F,40, 1.5,33 


Width  of  supported  edge  Buckling  load,  lbs. 


1 1/4 

715 

1 

755 

1 1/4  (a) 

670 

1 1/4  (b) 

650 

1 

675 

3/4 

730 

1/2 

600 

1 1/4 

760 

1 

740 

3/4 

715 

1 1/4 

815 

1 

625 

3/4 

770 

1/2 

640 

* The  designation  of  the  plate  number  is  as  follows:  F indicates 

non-preatressed  flat  plates,  the  first  set  of  numbers  indicates  the 
plate  thickness  is  .040  in.,  the  second  set  of  numbers  Indicates  that 
the  aspect  ratio  is  1.5,  while  the  last  series  of  numbers  indicates 
that  this  was  the  30th  of  the  type  tested. 
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TABLE  16 


Buckling  Tests  of  Flat  Plates 
With  l/2  Inch  Supported  Edges 

no. 

Buckline  load.  lbs. 

F, 40,1.5,22 

565 

23 

3 waves 

24 

560 

25 

610 

29D 

560 

31D 

600 

33D 

640 

Average  value  ■ 589  lbs. 

TABLE  17 

Buckling  Tests  of  Flat  Plates 
With  3/4  Inch  Supported  Edges 


Plate  No. 

Buckling  load,  lbs, 

F, 40,1.  5, 26 

510 

27 

580 

28 

670 

29A 

crimped  Plate 

31C 

730 

32C 

715 

33C 

770 

Average  value  = 663  lbs. 
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TABLE  18 


Plate  No. 


Buckling  Tests  of  Flat  Plates 
With  1 Inch  Supported  Edges 


Buckling  Load  lbs.  (Toshikl's  Method) 


.1 

750* 

2 

770* 

3 

760* 

4 

- 

5 

745 

6 

800* 

7 

- 

8 

765 

9 

745 

10 

705 

11 

755 

12 

755 

13 

815 

14 

745 

15 

600 

16 

3 waves 

17 

650 

18 

730 

19 

680 

20 

3 waves 

21 

690 

30B 

755 

31B 

675 

32B 

740 

33B 

625 

Average  value  * 726 

* Gage  poorly  located  - determined  buckling  load  by  load  vs.  bending 
graph,  not  shown  in  Fig.  15* 
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TABLE  19 


Buckling  Tests  of  Ordinary  Flat  Plates 
With  1 l/A  Inch  Surroorted  Edges 

Plate  No. 

Buckling  Load  lbs. 

F,40,1.5,30A 

715 

31A 

660 

32k 

760 

33k 

815 

Average  Value  « 738  lbs. 


TABLE  20 

Theoretically  Calculated  Loads  Carried  by  Supported  Widths 
Width  of  supported  strip  Load  carried  by  strip 


l/2  inch 

63  lbs 

3/4  inch 

95  lbs 

1 inch 

126  lbs 

1 l/4  inch 

158  lbs, 
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TABLE  21 

Buckling  Tests  of  Non-tire  stressed  Flat  Plates.  .032  Inch  Thlcknesi 
With  1 Inch  Supported  Edges 


Plate  No. 


F,  32, 1.5,1 


2 

3 

4 

5 

6 

7 

8 


Buckling  Load  (Toshlkl’s  Method) 
320 
300 
400 
365 
350 
255 
270 
245 


Average  buckling  load  313  lbs. 
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TABLE  22 


Buckling  Tests  of  Non-Prestressed  Flat  Plates.  .051  Inch  Thickness 


F, 51,1.5,1  980 

2 1490 

3 1150 

4 1100 

5 1200 


Average  buckling  load  1175 


WADC  TR  54-8 


105 


TABLE  23 


Buckling  of  .040  Prestressed  Plates 
Initial  Radius  Buckling  Load  Buckling  Behavior 


Plate  No. 

Inch 

Ton  of  knee 

Toshiki’s  Method 

0,40,1.5,27 

5 

800 

? 

Mode  jump 

28 

? 

tt  ft 

29 

900 

810 

Normal 

30 

700 

740 

»i 

31 

700 

740 

Mode  jump 

32 

Approx.  700 

7 

tt  it 

33 

Approx.  1900 

680 

34 

900 

980 

Normal 

35 

Approx.  1200 

7 

Mode  jump 

36 

5 

900 

860 

Normal 

0,40,1.5,2 

6 

? 

1480 

it 

3 

1800 

1530 

n 

4 

900 

1070 

n 

13 

900 

770 

Mode  jump 

24 

500 

7 

Normal 

25 

750 

790 

n 

37 

650 

580  or  2180 

Double  Parabola 

38 

800 

710  or  1460 

n n 

39 

650 

810  or  920 

n n 

40 

700 

640  or  1840 

it  ii 

a 

550 

7 

Normal 

42 

6 

700 

650  or  1580 

0,40,1.5,5 

7 

850 

7 

Normal 

6 

1150 

7 

Mode  jump 

7 

900 

880 

Normal 

8 

850 

850 

n 

9 

700 

760 

tt 

10 

600 

620 

it 

11 

7 

650 

620 

it 

12 

7 

650 

620 

n 
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TABLE  23  (continued) 


Plate  No. 

Initial  Radius 
Inch 

Buckling  Load 

Top  of  knee  Yoshiki's  Method 

Buckling  Behavior 
Type 

0,40,1.5,16 

8 

900 

890 

Normal 

17 

700 

820 

n 

18 

1100 

910 

u 

19 

750 

710 

n 

20 

1100 

? 

Mode  jump 

21 

1050 

980 

Normal 

22 

800 

820 

ti 

23 

8 

600 

610 

11 

0 ,40,1.5,43 

10 

800 

480 

n 

44 

800 

540  and  1900 

Double  Parabola 

45 

800 

580 

Normal 

46 

800 

640 

ir 

47 

800 

? 

n 

48 

950 

700 

n 

49 

800 

450 

n 

50 

10 

1000 

680 

n 
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Plate 

C,32, 

c,32,: 


C,32, 

c,32,: 

* Not 


TABLE  24 

Buckling  of  .032  Inch  Prestressed  Plates 


Initial  Radius  Buckling  Load  Buckling  Behavior 


No.  of  Curvature  Ton  of  knee 


1.5,1 

5 

500 

2 

5 

450 

3 

5 

1700 

4 

5 

300 

5 

5 

550 

6 

5 

200 

7 

5 

- ■ 

8 

5 

- 

9 

5 

250 

1.5,10 

6 

450 

11 

6 

250 

12 

6 

(400+)  * 

13 

6 

400 

14 

6 

400 

15 

6 

250 

32 

6 

350 

33 

6 

200 

34 

6 

200 

35 

6 

400 

36 

6 

375 

37 

6 

- 

L.5,16 

7 

375 

17 

7 

400 

18 

7 

400 

L.5,19 

7 

- 

20 

7 

400 

21 

7 

400 

22 

7 

500 

23 

7 

200 

used  in  finding  averages 

• 

Toshikl 1 s Method  Tyne 


510 

Normal 

430  or  640 

Double  Parabola 

1650 

Normal 

300  or  870 

Double  Parabola 

400 

Mode  jump 

215 

Normal 

250 

it 

220  or  710 

Double  Parabola 

225  or  610 

W If 

440 

Normal 

230 

tt 

340 

Mode  jump 

430  or  1060 

Double  Parabola 

380  or  1190 

Normal 

- 

if 

170 

Mode  jump 

220 

Normal 

200 

tt 

320 

tt 

270 

Mode  jump 

- 

it  n 

- 

Normal 

350 

tt 

460 

tt 

400 

tt 

- 

Mode  jump 

390  or  760 

Double  Parabola 

330 

Normal 

250 

Normal  (3  half  wave 
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TABLE  24  (Continued) 


Initial  Radius  Buckling  Load  Buckling  Behavior 


Plate  No, 

of  Curvature 

Top  of  knee 

Yoshiki’s  Method 

Type 

0,32,1.5,24 

8 

350 

300 

Normal 

25 

8 

375 

320  or  490 

Double  Parabola 

26 

8 

550 

385 

Normal 

27 

8 

- 

300 

Mode  jump 

28 

8 

525 

465 

Normal  * 

29 

8 

400 

360  or  530 

Double  Parabola 

30 

8 

- 

330 

Mode  jump 

31 

8 

350 

290 

Normal 

* Not  used  in  finding  averages. 
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TABLE  25 


Buckling  of  .051  Inch  Prestressed  Plates 


Initial  Radius  Buckling  Load  Buckling  Behavior 


Plate  No. 

of  Curvature 

Yoshiki’s  Method 

0,51,1.5,1 

6 

2200 

? 

Mode  jump 

2 

6 

2000 

? 

n it 

3 

6 

1400 

? 

n n 

4 

6 

2000 

? 

it  it 

5 

6 

2800 

? 

n it 

6 

6 

1400 

1200 

Normal 

7 

6 

1900 

1200  and  2000 

Double  Parabola 

8 

6 

1900 

1500 

Normal 

0,51,1.5,9 

7 

1400 

1420 

it 

10 

7 

1700 

1400  and  2700 

Double  Parabola 

11 

7 

1400 

1220 

Normal 

12 

7 

1600 

1400 

n 

13 

7 

• 1400 

1280 

it 

H 

7 

2000 

? 

Mode  jump 

15 

7 

1400 

1790 

Normal 

16 

7 

1400 

1200  and  2800 

Double  Parabola 

0,51,1.5,17 

7 3/4 

1200 

900 

Normal 

18 

7 3/4 

1800 

? 

Mode  jump 

19 

7 3/4 

1500 

1000 

Normal 

20 

7 3/4 

1500 

1080 

n 

21 

7 3/4 

1400 

1250 

n 

22 

7 3/4 

1600 

1200 

it 

23 

7 3/4 

1600 

? 

it 

24 

7 3/4 

1800 

? 

n 

0,51,1.5,25 

9 

1300 

1150 

ft 

26 

9 

1500 

1180 

n 

27 

9 

1400 

1210 

rt 

28 

9 

1450 

1160 

ft 

29 

9 

1600 

? 

ft 

30 

9 

1500 

1220 

tt 

31 

9 

1500 

1220 

n 

32 

9 

1200 

1060 

rt 
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TABLE  28 


Avera ge  Buckling  Loads  of  .051  Inch  Prestressed  Plates 
Initial  Top  of  Knae  Method 

Radius.  In.  Average  Load  No.  of  Tests  Averaged  Average  Load  No.  of  Teats  Av. 

6 1950  8 1300  or  1567  3 

7 1538  8 1387  or  1801  7 

7 3/4  1550  8 1086  5 

9 1469  8 1163  6 
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TABLE  29 

Peroent  Ihorease  of  Buckling  Load 

Initial  Radius  Net  Buckling  * Increase  of  Prestressed  Over 

of  Curvature  Thickness  Load  Ordinary  Plate  (Nat  Buckllgg  Loadp) 


5 

.032 

499  (405  or  589) 

101*  (63*  or  137*) 

6 

.032 

263  (235  or  379) 

67*  (-5*  or  53*) 

7 

.032 

317  (298  or  360) 

28*  (20*  or  45*) 

8 

.032 

360  (279  or  321) 

45*  (12*  or  29*) 

.032 

(248)* 

- 

5 

.040 

8a  (676) 

40*  (13*) 

6 

.040 

683  (777  or  1236) 

47*  (29*  or  106*) 

7 

.040 

687  (619) 

15*  (3*) 

8 

.040 

749  (694) 

25*  (16*) 

.040 

(600) 

25*  (16*) 

10 

.040 

718  (455  or  652) 

19.66*  (-24*  or  8.66*) 

.040 

(600)* 

- 

6 

.051 

1426  (1038  or  1305) 

56.2*  (13.7*  or  42.9*) 

7 

.051 

1205  (1125  or  1539) 

39.64*  (23.22*  or  68.56*) 

7 3/4 

• 

o 

ut 

H 

1288  (824) 

a*  (-9.74*) 

9 

.051 

1207  (901) 

32.2*  (-1.3*) 

.051 

(913)* 

mm 

* Used  for  purposes  of  comparison 
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TABLE  30 

Coefficients 


0,2  " 2g#44  ( ^1  ,l2+lgw3, 1+50w5,  1“4w1  ,1W1 , 5“36w3 , 1W3 , 3 


^w1,5W1,7“4w1,7w1,9”4w1,9W1,11W1,13"36w3,3W3,5 


+98“7,1+162"9,12+242“U,12+338U13,X2'36"3,5U3,7) 


£ 

bn  / “ nwi  o-16w,  ,vn  _-16w-  ,v,  _-l6v,  rv,  rt-l6w. 


0,4 


455,1 


1,1  1, 3 1,1  1,5  1 , 3 1,71 , 5 1,91 


"l6wl,9Wl,13“l6wl,llWl,15+l44w3 ,1W3 ,3“144w3 ,1W3, 5 


-144w3,3W3,7  . . . ) 


E 


0,6  * — -36“l,3W1>9 


-36“l.5Ul.ll-36“X,7Ul,13-36“l>9“l,15*324“3.1U3.5 


'324u3,1u3,7+162u3,3  -324w3,3w3.9  ' * • > 


rt  * 


E 


0,8  ?282 


(64w1,1W1,7“64w1,1w1,9+64w1,3W1,5“64w1,3W1,1I 


“64w1,5W1,13”64w1,7w1,15+576w3,1w3,7‘576w3,1W3,9  * * * \ 
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,7W1,11 


WADC  TR  54-8 


£ 


O.io  * ~ (lOOtf-  _w.  Q-100w-  «w_  - _ +100w.  _w_  _-»100w-  _w 
* v 17778  X»1  1?9  1»1  1»11  1,3  1,7  1,3  1,13 


*5°“i.5  -100»1i5>'1>15^OOu3>1v3i^OO»3j1w3(U 


4900i,3,3w3,7  • • • > 


M2  ‘ ^ (K4W1.1W1,11-I«’«l.l’'lli3^1,3‘'if9-U4w1>3v1(15 


*1^1>5'<1.7*1296“3.1''3,11-1S96>'3,1U3.13+1296“3,3V3,9 

*1296w3,5v3,7  • • • ) 

(arl.l2-4vl,lw3,l+18“l,32-36wl,3w3,3 

t50“1,52-100“1,5V3.5t98“1,72-1’S,7“3,7 


+162“l,9  <Z42w1,U2*338uX>132-^“3(1''5i1  . . . ) 


E 


'2>2  " 300.4  (16“l.lV1.3+16“l>lW3tl^4“1,3u1.5-64w1,3U3.1 


■1S.3U3,5',:U4ul,5Ul,7'144“l,5U3,3^4wl,5W3,7 


-576,,l,9V3,3-784''l,nw3,5  * • * > 
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WADC  TR  54-3 


TABUS  30  (Contlmcd) 


V mi 


+100wl,3Wl,7+100wl,  3W3,1^W1,3W3  ,7+196w1,  5W1,9 

"196wl,  5V3,l“36wl,  5W3,9*324wl,7Wl,n“324wl,7W3, 3 

+484wl , 9W1 , 13^84w1 , 9W3 , 5+676w1  , UW1 , 15+196w3 , 1W5 , 3 * * # * 


”2,6  • zt  (-16'l,lWl,5+64vl,l’'l,7+64vl,lV3,5-1S.lW3,7 


3600 


+144wl,  3W1,9+U4W1, 3W3 , 3+256w1,  5Wl,n+256wl,  5W3 ,1 
“1S,5W3,n+400wl,7wlA3^°^l,7w3,l^6w1.9Wl,15 


"576wl,9W3,3"784vl,nW3, 5 * 


. ) 
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TABUS  30  (Continued) 


9474 


^^1,1W1»7+100w1,1W1,9+10CV1,1W3,7‘36w1,1W3,9 


^vl,  3W1 , 5+196w1 , 3V1 , ll+196wl , 3W3 , 5"^W1 ,3*3,11 

+324v1,5W1,13+324w1,5W3,3^4w1,5V3,13+4S4w1,7W3,1 

+4^1,7W1A5^1,9W3A'J9^1,UV3>3  * * • > 
b2,10  " 2U20  (^W1AV1,9+^WI,1W1,11+U4w1,1W3,9"64w1,1W3,11 

”16w1,3W1,7+256w1,3V1>13+256w1,3W3,7“16w1,3V3,13 

+40^l,5^,15+4^l,5W3t5+^,7w3,34784wl,9W3,l 


_1(^wl  ,nW3,1‘"l296w1  >13W3, 3“l600w1 ,1.5W3 , 5 * * * * 


b2,12  * a6l6  (*’100wl,lWl,ll+196wl ,1V1,13+196w1  ,1*3  ,ll“100lfl,  1*3 ,13 


“36vl,3Wl,9+324wl, 3W1,15+324v1,3W3,9^4w1, 5W1,7 


+^vl>5W3,7+67^lf7w3>5^1,9W3>3+n5SfllW3,l 


“J44wl,13W3,l“176wl,15W3> 3 * * * \ 
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TABLE  30  (Continued) 


E 


b4,°  “ 23^  *l6wl ,1W3 ,l"l6wl,  1W5 ,l+144l,l, 3W3 , 3~l44wl 


U(XVl,5V3,5+mwl,7w3.7  * • * > 

b4,2  * ^ (“4wl ,1W3 , , +36wL ,1W3, 3+36wI,1W5 , 1‘ 


,-43 


+100,'lf3w3,l+19fe'l,3w3,5_196wl,3w5fl't324wl, 

t<S4u1.5V3.7-484w1.5v5,3*676“l,7u3,5 


+1156w^  ij|  * • » ) 


“4,4 


E_ 

4807 


(64ul,lW3,5*64“l>lU5,3-64”l,3U3,l+256“l 


+256k1.3W5,1*256"1,5w3,1+576"i,5M3,9 


■576wl,5w5,l',576wli7v3,3'fl024wl,9w3,5 


'4,6  ' Zt  ('4kL.1W3.5+100'i,1w3,7+100“i,1W5,5-36w1, 


9216 


,3W5,3 

5,3 

5W3,3 

,3W3,7 

) 

>3W3,3 
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TABLE  30  (Contlntwd) 


+324,'l,3v3,9*324wl,3’'5,3-19S,5w3,l*476v1.5w3,U 


t676“l,5w5,l*484wl,7"3.1t900w1.9w3,3*u4”l,llv3.5  * • • > 


'4,8  * * (‘16wl,lW3,7+144wl,lH3,9'l6“l.,3v3,5+400wl,3''3,U 


17778 


+784*1 ,9W3 , l+l296wi , nw3 , 3 • • • ) 


E 


6,0  ' (36ul,l“5,l+324wl,3',5>3tl8w3,l  +162w3(3  . . . ) 


6.2  ‘ dt  (-1S,l''5.1^4u1>lK5.3t256w1.3u5,l"784w1.5W5.3 


12844 


+L44w3,1W3,3+576w3,3W3, 5 * * * > 


6.4  ‘ WW7  ('4wl.lW5,3+100“1>lw5.5-196“l<3'‘5,lt47S.5W5,l 


'36*3,lw3,3'f324''3,lM3,5'f900w3,3w3,7  • • • > 
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V - ~ (-U*'l.3v5.3-576^.5w5,l-:UAM3,lV3,5*576*3,l''3.7 


♦I296»3j3»?>9  • • • ) 


36864 


^64w3,1W5,1  * * * ^ 


b8,2  " 38940  (’'4v3,1V5,1+196w3,1V5,3+324w3,3V5,1  * * * 5 


45511 


("16w3,1W5,3”:U4w3,3W5,1  * * * * 


hlA  A * ~ (50W.  . • • • ) 

omm  5»1 
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TABLE  31 


D1 

Coefficients 

d2  d3 

°4 

D5 

D6 

w 3 
1,1 

1.451 

-.02705 

-.001497 

0 

0 

0 

vl,lvl,3 

-.972 

.448 

• 01149 

-.00805 

-.01488 

0 

V191V3,1 

-.1920 

.0410 

.264 

-.0364 

0 

-.000711 

w1,12w3,3 

0 

-.0542 

-.0689 

.1298 

.00957 

.000843 

vl,lvl,5 

0 

-.0875 

0 

.00832 

.1492 

0 

vl,ly5fl 

0 

0 

-.00517 

.00325 

0 

.0945 

W1,1W1,32 

5.38 

0 

-.0642 

0 

.0417 

0 

V1,1W3,12 

11.30 

—.414 

0 

0 

0 

.00695 

w w 2 
1,13,3 

10.50 

0 

0 

0 

.0751 

.0213 

w w 2 
Wl,l  1,5 

10.50 

0 

-.1382 

0 

0 

0 

w w 2 
Wl,l  5,1 

29.35 

-.811 

0 

0 

0 

0 

W1,1W1,3W3,1 

.983 

-.459 

-.232 

.205 

.0604 

.00918 

W1,1W1,3W3,3 

-1.301 

0 

.387 

0 

-.0681 

-.0247 

W1,1W1,3W1,5 

-2.095 

.490 

.0992 

-.0593 

0 

0 
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TABUS  31  (Continued) 


Di 

D2 

D3 

D4 

D5 

D6 

W1,1W1,3W5*1 

0 

0 

.0667 

-.0948 

0 

—.0626 

W1,1W3,1W3,3 

-5*89 

1.382 

0 

0 

-.2015 

-.01457 

W1,1W3,1V1»5 

0 

.354 

0 

-.1753 

-.1680 

0 

W1,1W3,1W5,1 

-.443 

.238 

.1009 

-.0559 

0 

0 

W1,1W3,3W1,5 

1.342 

-.400 

-.332 

.1307 

0 

.0271 

W1,1W3,3W5,1 

.525 

-.642 

-.1058 

.1637 

.1200 

0 

W1,1V1,5W5,1 

0 

0 

0 

.1041 

0 

0 

wl,33 

0 

1.632 

0 

-.0640 

0 

0 

W1,3^W3,1 

-2.755 

0 

.740 

0 

-.1118 

-.0209 

w1,32w3,3 

0 

-1.297 

0 

.668 

0 

0 

w1,32w1,5 

2.935 

0 

-.1837 

0 

.848 

0 

w1,32w5,1 

0 

0 

-.1520 

0 

0 

.1998 

W1,3W3,12 

-4.96 

2.64 

0 

0 

-.2141 

-.0253 

W1,3W3,32 

0 

4.52 

0 

0 

0 

0 
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TABLE  31  (Contlntted) 

D1  D2  D3  D4  D5 


W1,3W1,5 

0 

4.99 

w v 2 
Wl,3  5,1 

-9.72 

5.17 

W1,3W3,1W3,3 

16.59 

0 

W1,3W3,1W1,5 

4.25 

-1.310 

W1,3W3,1W5,1 

2.855 

-1.087 

W1,3W3,3V1,5 

-4.78 

0 

W1,3W3,3W5,1 

-7.68 

0 

W1,3W1,5W5,1 

0 

0 

w 3 
3,1 

0 

0 

W3,12v3,3 

0 

0 

w3,12w1,5 

0 

-1.257 

w3,12w5,1 

2.155 

-.655 

W3,1W3,32 

0 

0 

w w 2 

w3,l  1,5 

-5.92 

0 

v w 2 
V3,1W5,1 

0 

0 

WADC  TR  54-8 


0 

-.397 

0 

0 

0 

-.431 

0 

0 

.517 

-.704 

.450 

0 

-.367 

.465 

•286 

.850 

0 

-.915 

to 

& 

• 

0 

-.456 

.469 

-.396 

0 

2.02 

-.324 

0 

-1.839 

2.165 

0 

0 

0 

1.047 

0 

0 

0 

4.09 

0 

0 

1.719 

0 

0 

5.76 

-1.120 

0 

123 


D6 

0 

0 

♦1209 

.0647 

0 

-.1031 

0 

-.1952 

0 

0 

0 

.793 

0 

-.0886 

0 


v-  -W,  ,w_  _ >14*18  3.040 

3fl  3,3  lf5 


w3,1,'3.3V5,1  "4*53  3,135 


-.733  -.585 


w3,1W1»5W5,1 


0 1.679 


-.637  -.782 


w3,32w5,l 


6.63 


1.451 


W3,32w1,5 


1.850  .1281 


W3,3W1,5 


0 -2.685 


1.608 


W3,3W5,1 


-2.12 


3.50 


3f3  1,5  5»1 


8*44  -2*68  -1*205 


*1,5%, 1 


-.642 


1.920 


W1,5W5,1 


0 -2.53 


2.51 


2.12 
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TABLE  31 
®7  ” D12* 

D7  DS 

( continued) 

°9 

D10 

Du 

D12 

w 3 
1,1 

0 

0 

0 

0 

p 

0 

w1,12w1,3 

0 

0 

0 

0 

0 

0 

v1,1m3,1 

0 

0 

0 

0 

0 

.00000802 

Mlflv3,3 

0 

0 

0 

0 

0 

-.001499 

wl,l  wl,5 

0 

0 

0 

0 

0 

0 

v 2w 
Vl,l  W5,l 

-♦000203 

0 

0 

0 

0 

-.02025 

W1,1W1,32 

0 

0 

0 

0 

0 

0 

v w 2 
wl,l  3,1 

-.000233 

0 

0 

0 

0 

-.001217 

v w 2 
wl,l  3,3 

-.00250 

0 

0 

0 

0 

0 

v v 2 
1,1  1,5 

0 

0 

0 

0 

0 

0 

v w 2 
Wl,l  5,1 

0 

.000509 

-.0000531 

0 

0 

0 

W1,1W1,3W3,1 

0 

0 

0 

0 

0 

-.00650 

W1,1W1,3W3,3 

0 

0 

0 

0 

0 

0 

W1,1W1,3W1,5 

0 

0 

0 

0 

0 

0 
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TABLE  31  (Continued} 
R,  Dg  D, 

D10 

D11 

D12 

wl,lwl, 3*5,1 

.00462 

0 

0 

0 

0 

.0931 

wl,l*3,lw3,3 

.001197 

0 

0 

0 

0 

.01295 

W1,1W3,1W1,5 

0 

0 

0 

0 

0 

•01490 

W1,1W3,1W5,1 

.00471  - 

-.0002050 

0 

0 

0 

0 

W1,1W3,3W1,5 

0 

0 

0 

0 

0 

-.01303 

/3,3W5,1 

-.00216 

.000730 

0 

0 

0 

0 

W3A 

0 

0 

0 

0 

o 

-.0616 

w 3 
1,3 

0 

0 

0 

0 

0 

0 

vl»3^w3,l 

0 

0 

0 

0 

0 

0 

*1,3%,  3 

0 

0 

0 

0 

0 

-.0404 

V1932v195 

0 

0 

0 

0 

0 

0 

v 2w 
"1,3  5,1 

-.0C718 

0 

0 

0 

0 

0 

W1,3W3,12 

-.001727 

0 

0 

0 

0 

.0315 

W1,3W3,32 

0 

0 

0 

0 

0 

.0714 

W1,3W1,52 

0 

0 

0 

0 

0 

0 
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TABLE  31  (Continual 


D? 

D8 

D9 

D10 

hi 

°12 

w v 2 
1, 3*5,1 

0 

-.002445 

.000416 

0 

0 

0 

W1,3W3,1W3,3 

-.01307 

0 

0 

0 

0 

0 

w1,3W3,1W1,5 

0 

0 

0 

0 

0 

-.0391 

w1,3w3,1W5,1 

-.01845 

.001781 

0 

0 

0 

0 

W1,3W3,3W1,5 

0 

0 

0 

0 

0 

0 

*1,3*3, 3*5,1 

.0296 

-.00582 

0 

0 

0 

0 

*1,3*1,5*5,1 

.0313 

0 

0 

0 

0 

.1491 

w 3 

*3,1 

0 

-.0000201 

0 

0 

0 

0 

¥3,12w3,3 

0 

.000292 

0 

0 

0 

0 

w3,12w1,5 

0 

0 

0 

0 

0 

-.0493 

*3,1%, 1 

0 

0 

-.0000288 

0 

0 

-.1859 

W3,1W3,32 

0 

-.001701 

0 

0 

0 

0 

w w 2 
*3,1*1, 5 

0 

0 

0 

0 

0 

0 

W3,1W5,12 

.01522 

0 

0 

-.00001529 

0 

0 

W3,1V3,3W1,5 

.02185 

0 

0 

0 

0 

•0822 
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TABLE  31  (Continued) 


D_ 

7 

D8 

D9 

D10 

D11 

D12 

W3,1W3»3W5»1 

0 

0 

.000374 

0 

0 

.823 

w3,1W1,5W5,1 

0 

0 

0 

0 

0 

0 

3 

w3,3 

0 

0 

0 

0 

0 

0 

w3,32w1,5 

-.01810 

0 

0 

0 

0 

0 

*3,3%, 1 

0 

0 

-.001024 

0 

0 

0 

W3,3V 1,52 

0 

0 

0 

0 

0 

-.1239 

W3,3W5,12 

-.0339 

0 

0 

.0001128 

0 

0 

W3,3W1,5W5,1 

-.0697 

.01172 

0 

0 

0 

0 

3 

Wl,5 

0 

0 

0 

0 

0 

0 

w1,52w5,1 

-.0393 

0 

0 

0 

0 

0 

W1,5W5,1 

0 

0 

0 

0 

0 

0 

t 3 

0 

0 

0 

0 

-.00000263 

-.502 

t 
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TABUS  31  (continued) 

<*13  - V 


D13 

Du 

D15 

D16 

D17 

D18 

wl,l3 

0 

0 

0 

0 

0 

0 

W1»1^W1,3 

0 

0 

0 

0 

.0000427 

0 

wl,i%,l 

0 

0 

0 

0 

0 

0 

wl*l^w3f3 

0 

0 

0 

0 

-.01541 

0 

w1,12w1,5 

0 

0 

0 

0 

-.00948 

0 

wl,l  v5,l 

.00000851  0 

0 

0 

0 

0 

W1,1W1»32 

0 

0 

0 

0 

-.00253 

0 

W1,1W3,12 

.00000241  0 

0 

0 

0 

0 s 

v v 2 
1,1  3,3 

0 

0 

0 

0 

0 

0 

v w 2 
Wl,l  1,5 

0 

0 

0 

0 

0 

0 

W1,1W5,12 

0 

-.0000757 

.000001522 

0 

0 

0 

W1,1W1,3W3,1 

0 

0 

0 

0 

-.0371 

.000461 

W1,1W1,3W3,3 

0 

0 

0 

0 

.0451 

-.001522 

WL,1W1,3V1,5 

0 

0 

0 

0 

0 

0 
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TABLE  31  (Continued) 


Du 

D15 

Dl6 

D17 

^8 

W1,1W1,3W5,1 

-.00368 

0 

0 

0 

.0227 

.000510 

W1,1W3,1W3,3 

-.001042 

0 

0 

0 

0 

-.00318 

w1,1w3,1W1»5 

0 

0 

0 

0 

.1190 

0 

W1,1W3,1W5,1 

-.001486 

.00000902 

0 

0 

0 

0 

W1,1W3,3W1,5 

0 

0 

0 

0 

0 

0 

W1,1W3,3W5,1 

.00793  - 

-.000534 

0 

0 

-.0522 

0 

W1,1W1*5W5,1 

.00897 

0 

0 

0 

-.0972 

.0725 

w 3 
wl,3 

0 

0 

0 

0 

0 

0 

w1,32w3,1 

0 

0 

0 

0 

.0456 

-.00323 

w1,32w3,3 

0 

0 

0 

0 

0 

0 

w1,32w1,5 

0 

0 

0 

0 

-.0847 

0 

w1,32w5,1 

0 

0 

0 

0 

-.0401 

.0266 

W1,3W3,12 

-.001942 

0 

0 

0 

0 

-.00828 

2 

W1,3W3,3 

-.01305 

0 

0 

0 

0 

0 
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TABLE  31  (Confirmed) 


D13 

DU 

D 

15 

°16 

D,_ 

17 

DX8 

W1,3W1,52 

0 

0 

0 

0 

0 

0 

W1,3W5,1^ 

0 

*00568 

-.000708 

0 

0 

0 

w1,3W3,1W3,3 

0 

0 

0 

0 

0 

.01880 

wlt3W3,lWl,5 

0 

0 

0 

0 

0 

0 

W1,3W3,1W5,1 

.0312 

-*002445 

0 

0 

-.1507 

0 

W1,3W3,3W1,5 

0 

0 

0 

0 

.467 

-.0599 

W1,3W3,3W5»1 

0 

0 

0 

0 

.1499 

0 

W1,3W1,5W5»1 

“■o0203 

0 

0 

0 

0 

0 

3 

w3,l 

0 

0 

0 

0 

0 

0 

w3,l  w3*3 

0 

-,0002820 

0 

0 

— ©408 

0 

w3,12w1,5 

2 

*00582 

0 

0 

0 

0 

.0449 

W3,l  w5,l 

0 

0 

.0000000316  0 

0 

0 

W3,1W3,3 

0 

0 

0 

0 

.426 

0 

w3,1W1,52 

0 

0 

0 

0 

0 

0 

w w ^ 

3,15,1 

-0 00583 

0 

0 

.0000000165  0 

0 

/ 


/ 
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Jontinuec 


v-  ,w_  awn  _ -.01328  0 

3,1  3f3  If 5 


v3,1W3,3W5,1 


-.000378  0 


W3,1W1,5W5,1 


-.0496  .00730 


W3»3^W1»5 

w3f3^w5»l 


W3,3W1,5 


•.001499 


W3,3W5,1 


.0738 


-.0001408  -.803 


W3,3W1,5W5,1 


,0465  -.00822 


W1,5^W5»1 


W1,5W5,1 


0 -.01147  .00209 
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D19 

TABLE  31 
<D19  - V 

D20  D21 

D22 

it 

D23 

B24 

V 3 
1,1 

0 

0 

0 

0 

0 

0 

wl,l  Wl,3 

0 

0 

0 

0 

0 

0 

wl,l  w3,l 

0 

0 

0 

0 

0 

0 

wl,l  W3,3 

0 

0 

0 

0 

0 

0 

wl,l  wl,5 

0 

0 

0 

0 

-.00456 

.0000900 

wl,l  W5,l 

0 

0 

0 

0 

0 

0 

wl,lwl,3 

0 

0 

0 

0 

-.00474 

.00001822 

w w ^ 

Wl,l  3,1 

0 

0 

0 

0 

0 

0 

wl,lw3,32 

-.0001518 

0 

0 

0 

-.0360 

0 

W1,1W1,5 

0 

0 

0 

0 

0 

0 

W1,1W5,12 

0 

0 

0 

0 

0 

0 

W1,1W1,3W3,1 

0 

0 

0 

0 

0 

0 

W1,1W1,3W3,3 

0 

0 

0 

0 

.01611 

-.01428 

W1,1W1,3W1,5 

0 

0 

0 

0 

.0332 

-.00454 
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TABLE  31  (Continued) 


D19 

D20 

D21 

D22 

D23 

D24 

W1,1W1»3W5,1 

.000510 

0 

0 

0 

0 

0 

W1,1W3,1W3,3 

.00001970  0 

0 

0 

0 

0 

W1,1W3,1W1,5 

0 

0 

0 

0 

.03245 

-.02a5 

W1,1W3,1W5,1 

0 

0 

0 

0 

0 

0 

W1,1W3,3W1,5 

0 

0 

0 

0 

-.0272 

.0251 

w1,1W3,3W5,1 

-.00379 

.0000325 

0 

0 

0 

0 

W1,1W1,5W5,1 

-.00804 

0 

0 

0 

0 

.0228 

wl,33 

0 

0 

0 

0 

0 

0 

w1,32w3,1 

0 

0 

0 

0 

.0208 

-.01120 

W1,32w3,3 

0 

0 

0 

0 

0 

0 

w1,32w1,5 

0 

0 

0 

0 

0 

0 

w1,32w5,1 

-.001700 

0 

0 

0 

0 

.01730 

W1,3W3,12 

.000284 

0 

0 

0 

0 

0 

W1,3W3,32 

0 

0 

0 

0 

0 

0 

W1,3W1,52 

0 

0 

0 

0 

.0750 

-.01729 
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TABLE  31  (Continued) 


D, 


D, 


W1,3W5,1 


w1,3W3,1W3,3 


W1,3W3,1W1»5 


19  20  "21  "22 
0 -.001925  .0001578  0 


-.002005  0 


0 


0 


v,  ,w  -.01086  .000548  0 

J-fJ  Jp-L 


W1,3W3,3W1,5 


W1,3W3,3W5  1 •0l628  -.001613  0 


W1,3W1,5W5,1 


w3,l 

w3,l  w3,3 
w3,12w1,5 
W3,l  W5,l 


W3,1W3,3 


W3,1W1,5 


w3,lw5,l 


0 0 0 


.00001907  0 


-.00548  0 


0 


0 


0 


0 -.0002285  0 


0 


0 

0 

0 

0 

0 

0 


0 


23 

0 

-.1502 

-.1124 


0 

•1490 

0 


0 

-.1013 

0 

0 

0 

0 


D24 

0 

0 

.0698 

0 

0 

-.0667 

-.0819 

0 


0 

-.1757 

0 

0 
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TABLE  31  (Continued) 

D19 

D20 

D21 

D22 

D23 

D24 

W3,1W3,3W1,5 

0 

0 

0 

0 

.273 

0 

w3flW3,3W5,l 

0 

0 

.0000348 

0 

0 

0 

W3,1W1,5W5,1 

•0594 

-.00831 

0 

0 

.2105 

-.1268 

m3,33 

0 

0 

0 

0 

0 

0 

w3,32m1,5 

-.01320 

0 

0 

0 

0 

0 

w3,3  W5,l 

0 

0 

-.0001857 

0 

-.219 

0 

W3,3W1,52 

0 

0 

0 

0 

.1600 

.1152 

W3,3W5, l2 

-.0314 

0 

0 

00001683 

0 

0 

W3,3W1,5W5,1 

0 

0 

0 

0 

•.344 

•1475 

3 

Wl,5 

0 

0 

0 

0 

0 

0 

wl,5  W5,l 

0 

0 

0 

0 

0 

0 

W1,5W5,12 

0 

.01672 

-.00285 

0 

.2995 

0 

2 

0 

0 

0 

0 

0 

0 

9 
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D25 

TABLE  31 
(D25  ~ D30 

D26 

) 

D27 

D28 

D29 

D30 

V 3 

wl»l 

0 

0 

0 

0 

0 

0 

wl,l  wl,3 

0 

0 

0 

0 

0 

0 

wl,l  W3,l 

0 

0 

0 

0 

0 

0 

wl,l  W3,3 

0 

0 

0 

0 

0 

0 

wl,lS.,3 

0 

0 

0 

0 

0 

0 

w1,12w5,1 

0 

0 

0 

0 

0 

0 

W1»1W1,3 

0 

0 

0 

0 

0 

0 

W1»XW3,1 

0 

0 

0 

0 

0 

0 

wl,lw3,32 

-.000641 

0 

0 

0 

0 

0 

wl,lwl,5 

0 

0 

0 

0 

.00706 

-.000692 

W1,1W5,12 

0 

0 

0 

0 

0 

0 

W1,1W1,3W3,1 

0 

0 

0 

0 

0 

0 

W1,1W1,3W3,3 

.0000797 

0 

0 

0 

0 

0 

W1,1W1,3W1,5 

0 

0 

0 

0 

-.00418 

.0000860 
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TABLE  31  (Continued) 


D25 

D26 

D 27 

D28 

D29 

D30 

Wl,lwl,3w5,l 

0 

0 

0 

0 

0 

0 

W1,1W3,1W3,3 

0 

0 

0 

0 

0 

0 

W1,1W3,1W1,5 

.000920 

0 

0 

0 

0 

0 

W1,1W3,1W5,1 

0 

0 

0 

0 

0 

0 

W1,1W3,3W1,5 

-.00352 

0 

0 

0 

.01081 

-.00813 

W1,1W3,3W5,1 

0 

0 

0 

0 

0 

0 

W1»1W1,5W5,1 

-.0204 

.001021 

0 

0 

0 

0 

3 

wl,3 

0 

0 

0 

0 

-.000494 

0 

wl,3  w3»l 

.000468 

0 

0 

0 

0 

0 

*1,3%,  3 

0 

0 

0 

0 

.00505 

-.00376 

w1,32w1,5 

0 

0 

0 

0 

0 

0 

w1,32w5,1  • 

-.01268 

.000548 

0 

0 

0 

0 
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TABIE  31  (Continued) 

D25 

D26 

D-„ 

27 

D28 

D 29 

D30 

W1,3W3,1W3,3 

-•00796 

.000274 

0 

0 

0 

0 

w1,3W3,1W1,5 

-.00827 

0 

0 

0 

.0247 

-.01657 

W1,3W3,1W5,1 

0 

0 

0 

0 

0 

0 

W1,3W3,3W1,5 

0 

0 

0 

0 

0 

0 

V1,3W3,3W5,1 

0 

-.00808 

-.000491 

0 

0 

0 

W1,3W1,5W5,1 

.0538 

-.00580 

0 

0 

0 

.0328 

w 3 
w3,l 

0 

0 

0 

0 

0 

0 

w3,12w3,3 

0 

0 

0 

0 

0 

0 

v391m1,5 

-.00999 

.000997 

0 

0 

0 

0 

w3,lS,l 

0 

0 

0 

0 

0 

0 

w3,lw3,32 

0 

0 

.00000880 

0 

0 

0 

W3,1W1.52 

0 

0 

0 

0 

-.0303 

.0250 

w3flw5,l2 

0 

0 

0 

0 

0 

0 
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TABUS  31  (Continued) 


D25 

D26 

D27 

D28 

D29 

D30 

W3,1W3,3W1,5 

.0224 

-.00362 

0 

0 

-.0875 

0 

w3,1W3,1W5,1 

0 

0 

0 

0 

0 

0 

W3,1W1,5W5,1 

0 

-.01682 

.00204 

0 

0 

0 

w 3 
w3,3 

0 

0 

0 

0 

0 

-.0270 

w3,32wl,5 

0 

0 

0 

0 

0 

0 

w3,32w5,1 

-.0774 

0 

0 

.0000268 

0 

0 

W3,3W1,52 

-.02025 

0 

0 

0 

0 

0 

2 

W3,3W5,1 

0 

0 

0 

0 

0 

0 

W3,3W1,5W5,1 

0 

.01793 

-.00331 

0 

.1218 

-.0603 

w 3 
wl,5 

0 

0 

0 

0 

0 

0 

w1,52w5,1 

0 

0 

0 

0 

0 

-.0409 

W1,5W5,12 

0 

0 

-.00563 

.000862 

0 

0 

w 3 
5,1 

0 

0 

0 

0 

0 

0 

WADC  TR  54-6  HO 
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TABLE  33 
Coefficients 


(°i  - C?) 


C1 

C 

2 

°3 

C4 

°5 

C6 

'M3 

1.370 

-.001497 

0 

0 

0 

0 

0 

wl,l2wl,3 

.298 

-.01245 

0 

0 

0 

0 

0 

w1,12w3,1 

-.0690 

.155 

-.000687 

0 

0 

0 

0 

w1,12w3,3 

-.1148 

.2434 

-.00365 

0 

0 

0 

0 

W1,12w1,5 

.452 

-.02181 

0 

0 

0 

0 

0 

w1,12w5,1 

0 

.00458 

.0338 

-.000177 

0 

6 

0 

W1,1W1,32 

5.55 

-.0767 

0 

0 

0 

0 

0 

W1,1W3,12 

10.06 

0 

.00330 

-.000226 

0 

0 

0 

W1,1W3,32 

10.59 

0 

.01681 

-.00296 

0 

0 

0 

wl,lwl,52 

10.55 

—.1442 

0 

0 

0 

0 

0 

W1,1W5,1 

26.92 

0 

0 

0 

.000282 

-.0000485 

0 

W1,1W1,3W3,1 

-.092 

.1975 

-.00802 

0 

0 

0 

0 

W1,1W1,3W3,3 

-1.528 

.512 

-.0318 

0 

0 

0 . 

0 

W1,1W1,3W1,5 

-.430 

-.1096 

0 

0 

0 

0 

0 

W1,1W1,3W5,1 

0 

-.1042 

.0742 

-.00387 

0 

0 

0 
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TABLE  33  (Continued) 


C2 

°3 

°4 

°5 

C6 

C7 

W1,1W3,1W3,3 

-2.75 

0 

.00838 

-.001830 

0 

0 

0 

wl,lw3,lwl,5 

•449 

-.0860 

-.00736 

0 

0 

0 

0 

W1,1W3,1W5,1 

.271 

-.0668 

0 

.000252 

-.0001779 

0 

0 

wl,lw3,lwl,5 

.049 

.163 

-.0344 

0 

0 

0 

0 

W1,1W3,3W5,1 

-.801 

.1243 

0 

.00268 

-.000710 

0 

0 

W1,1W1»5W5,1 

0 

-.0141 

.0349 

-.00614 

0 

0 

0 

w 3 
1,3 

4.89 

-.192 

0 

0 

0 

0 

0 

2 

wl,3  w3,l 

-3.168 

.890 

-.0338 

0 

0 

0 

0 

wl,3^w3,3 

-3.85 

1.970 

-.1208 

0 

0 

0 

0 

wl,3^wl,5 

7.167 

-.6072 

0 

0 

0 

0 

0 

wl,3^w5,l 

0 

-.2314 

.2440 

-.01184 

0 

0 

0 

W1,3W3,12 

1.89 

0 

.0278 

-.00268 

0 

0 

0 

W1,3W3,32 

13.31 

0 

.1975 

-.03894 

0 

0 

0 

W1,3W1,52 

15.49 

-1.312 

0 

0 

0 

0 

0 

W1,3W5,12 

3.64 

0 

0 

0 

.00497 

-.000919 

0 

W1,3W3,1W3,3 

18.12 

0 

.1592 

-.02118 

0 

0 

0 

H4 
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TABLE  33  (Continued) 


C1 

C2 

C3 

°4 

C5 

C6 

°7 

W1,3W3,1W1,5 

-.24 

.985 

-.0978 

0 

0 

0 

0 

W1,3W3,1W5,1 

1.02 

.274 

0 

.0208 

-.002814 

0 

0 

W1,3W3,3W1,5 

-9.28 

3.129 

-.400 

0 

0 

0 

0 

W1,3W3,3W5,1 

-8.92 

1.161 

0 

.0544 

-.01045 

0 

0 

W1,3W1,5W5,1 

0 

-.997 

.4212 

-.05233 

0 

0 

0 

w 3 
3,1 

0 

1.048 

0 

0 

-.0000201 

0 

0 

W3,l2w3,3 

0 

2.616 

0 

0 

-.000499 

0 

0 

2 

w3,l  wl,5 

.75 

0 

.0060 

-.00296 

0 

0 

0 

w3,12w5,1 

.190 

0 

.235 

0 

0 

-.0000287 

0 

W3,1W3,32 

0 

4.99 

0 

0 

-.002782 

0 

0 

W3,1W1,52 

-6.11 

1.875 

-.1193 

0 

0 

0 

0 

W3,1W5,12 

0 

2.40 

0 

-.00227 

0 

0 

-.0000015 

W3,1W3,3W1,5 

-3.95 

0 

.1526 

-.0352 

0 

0 

0 

W3,1W3,3W5,1 

1.21 

0 

.582 

0 

0 

-.000586 

0 

W3,1W1,5W5,1 

2.600 

-.374 

0 

.03046 

-.00537 

0 

0 
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C1 

C2 

TABLE  33  (Continued) 
°3  °4 

C5 

C6 

C7 

W 3 
w3,3 

0 

4.11 

0 

0 

-.00450 

0 

0 

w3,3^wl,5 

14.32 

0 

.361 

-.0826 

0 

0 

0 

W3,3  W5,l 

7.93 

0 

.965 

0 

0 

-.001765 

0 

W3,3W1 t5 

-9.14 

5.60 

-.5116 

0 

0 

0 

0 

W3,3W5,1 

0 

4.37 

0 

.0305 

0 

0 

-.000225 

W3,3W1,5W5,1 

-.91 

2.237 

0 

.1084 

-.02459 

0 

0 

„ ^ 
Wl,5 

9.60 

-1.03 

0 

0 

0 

0 

0 

wl»5^w5  »1 

0 

-.844 

.708 

-.0631 

0 

0 

0 

W1,5W5,12 

2.86 

0 

0 

0 

.00978 

-.00426 

0 

w 3 
W5,l 

0 

0 

.614 

0 

0 

0 

0 

WADC  TR  54-8 
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TABLE  33  (Continued) 


w 3 
wl,l 

C8 

1.446 

C9 

-.02705 

°10 

0 

C11 

0 

°12 

0 

°13 

0 

cu 

0 

w1,12v1,3 

-.938 

.424 

-.01475 

0 

0 

0 

0 

wl,l  W3 A 

.596 

-.0682 

0 

0 

0 

0 

0 

wl,l  w3,3 

-.2025 

.328 

-.0367 

0 

0 

0 

0 

w1,12w1,5 

0 

-.0625 

.1208 

-.00429 

0 

0 

0 

w1,12w5,1 

.456 

-.0914 

0 

0 

0 

0 

0 

W1,1W1,3 

5.19 

0 

.0341 

-.00469 

0 

0 

0 

W1,1W3,1 

11.33 

-.420 

0 

0 

0 

0 

0 

W1,1W3,32 

10.59 

0 

.0740 

-.0392 

0 

0 

0 

V1,1K1,52 

10.08 

0 

0 

0 

.00498 

-.001064 

0 

wl,lw5#l2 

29.35 

-.812 

0 

0 

0 

0 

0 

W1,1V1,3W3,1 

.33 

.1235 

-.0486 

0 

0 

0 

0 

W1,1W1,3W3,3 

-.264 

0 

.0596 

-.02633 

0 

0 

0 

W1,1W1,3W1,5 

-1.797 

.3121 

0 

.0196 

-.00392 

0 

0 

W1,1W1,3W5,1 

-.0806 

.1553 

-.0710 

0 

0 

0 

0 

U7 
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C8 

C9 

o 

H 

O 

C11 

°12 

°13 

°I4 

W1,1W3,1W3,3 

-5.95 

1.440 

-.2173 

0 

0 

0 

0 

wl,lw3,lwl,5 

0 

-.0974 

.1305 

-.0294 

0 

0 

0 

W1,1W3,1W5,1 

-.109 

.0600 

0 

0 

0 

0 

0 

W1,1W3,3W1*5 

.482 

-.073 

0 

.0305 

-.01234 

0 

0 

W1,1W3,3W5,1 

.199 

-.100 

-.0628 

0 

0 

0 

0 

W1,1W1,5W5,1 

0 

.0671 

.0146 

-.0264 

0 

0 

0 

3 

wl,3 

0 

1.440 

0 

0 

-.000494 

0 

0 

wl,3  w3,l 

-•640 

0 

.00885 

-.0105 

0 

0 

0 

wl,3  w3,3 

0 

.510 

0 

0 

-.00600 

0 

0 

wl,3%.,5 

2.384 

0 

.594 

0 

0 

-.000705 

0 

w1,32w5,1 

.493 

0 

.0008 

-.0077 

0 

0 

0 

W1,3W3,12 

-5.07 

2.78 

-.2535 

0 

0 

0 

0 

W1,3W3,32 

0 

4.79 

0 

0 

-.0368 

0 

0 

W1,3W1,52 

0 ' 

3.79 

0 

.0231 

0 

0 

-.000376 

W1,3W5,12 

-9.74 

<* 

5.21 

-.447 

0 

0 

0 

0 
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C8 

°9 

TABLE  33  (Continued) 
°10  °11 

°12 

°13 

°u 

W1,3W3,1W3,3  17,10 

0 

.600 

-.1880 

0 

0 

0 

W.  « W_  2§46 

,2*X  Xj5 

-.156 

0 

.0556 

-.0180 

0 

0 

W1,3W3,1W5,1  1,641 

.504 

-.2372 

0 

0 

0 

0 

W1,3W3,3W1,5  ~2,75 

0 

.1865 

0 

0 

-.00744 

0 

W1,3W3,3W5,1  ^,89 

0 

.0931 

-.104 

0 

0 

0 

W1,3W1,5W5,1  ,65° 

-.585 

0 

-.018 

-.0029 

0 

0 

w3,l3  6.06 

-.972 

0 

0 

0 

0 

0 

“3,12“3,3  -5'51 

6.49 

-1.223 

0 

0 

0 

0 

^.lS.S  ° 

-1.463 

1.233 

—.1442 

0 

0 

0 

u3.lS.l  6J2 

-1.584 

0 

0 

0 

0 

0 

w3,lw3,3  12,25 

0 

1.276 

-.527 

0 

0 

0 

2 

w3  ^ “1*206 

0 

0 

0 

.0227 

-.00750 

0 

w3,lu5,l2  17*» 

-3.401 

0 

0 

0 

0 

0 

v3,1W3,3k1,5  -14*86 

3.36 

0 

.360  - 

■.1282 

0 

0 

W3,1W3,3W5,1  ”7,45 

7.25 

-2.035 

0 

0 

0 

0 
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W3,1W1,5W5,1 

w ^ 
w3,3 

w3,32w 1,5 
w3,32w5,1 
W3,3W1,52 
W3,3W5,12 

w3f3Wl,5W5,l 

w 3 
1,5 

V1,52v5,1 

W1,5W5,12 


C8 

0 

0 

5.78 

11.17 

0 

-6.59 


4.44 


0 

.639 

0 

10.60 


C9 

TABLE  33  (Continued) 

C10  cu 

°12 

°13 

°u 

-.514 

1.014 

-.270 

0 

0 

0 

4.34 

0 

0 

-.0810 

0 

0 

0 

2.015 

0 

0 

-.02834 

0 

0 

1.180 

-.606 

0 

0 

0 

1.52 

0 

.0844 

0 

- 0 

-.002889 

11.01 

-2.629 

0 

0 

0 

0 

.52 

0 

.195 

-.1153 

0 

0 

0 

1.300 

0 

0 

0 

0 

0 

0 

0 

-.0070 

-.00089 

0 

-2.61 

2.63 

-.3294 

0 

0 

0 

-2.51 

0 

0 

0 

0 

0 
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3473*5510  7595*1523  12235.519  16933.003  22159*176  26453.076  31252.328 
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Non  Linear  Terms 


K2- 

I!3- 

H4 

S' 

n 

.09894 

.032192 

-.011851 

.0273624 

.00704138 

-.306 

-.04684 

.04367 

.02588 

-.00761590 

-.0660 

-.237 

.0402 

.00881520 

.005427 

.1052 

♦073 

-.0838 

-.02150 

-.0018C7 

•1162 

-.00704569 

-.02217 

-.1244 

-.00151685 

-.0340709 

.02047 

-.01082 

-.00391952 

-.0906 

.338491 

.1619 

-.043514 

.0784 

.0231184 

.875 

.249522 

-.10613 

.196889 

.05005 

.619752 

.212407 

-.0732423 

I 

.1467 

.0278 

.633899 

.3341 

-.080868 

.219829 

.0462826 

.634673 

.532422 

.145061 

.457 

.252 

-.180 

-.0744 

-.00633 

-.1491^3 

-.357 

.040213 

.0242 

.0266 

-.406 

-.1576 

.0829 

.00165457 

-.0121891 

.0339956 

-.0727 

.1054 

! -.0C780271 

.0634 

-1.156 

-.162592 

.15480 

.1385 

-.0a26 

-.330 

-.0095535 

.1551 

.1947 

-.00161405 

-.196 

-.1106 

.0554 

.00916206 

-.00157925 
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-.683 

-.i:  08925 

1.16854 

1.640 

-1.25Q44 

-.161 

.150962 

3.54 

3.03659 
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.979 
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•591r95 

,U60 
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.418 
-.703 
-.886 
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{ Continued  ) 
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w*  2V?  o 
„-.;aa i 
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-.6-1812 
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h3 
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.153433 
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0 
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■ 

P.09 
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1.45 
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_#7rQpQ  7 
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B 
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—2.16 

0.365473 

-.0366604 

^.3u1.5w5.1 
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h3  ’ 
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TABUS  36 

Results  of  Rayleigh-Ritz  analysis  for  buckling 
parameter,  K,  using  selected  terms  of  series  (97) 

K - 


K 

Coefficients 

(Non-Prestressed 

Retained 

Wl,3 

9.91 

wlf2 

8.72 

K Percent  Increase  of  pre- 

(Pre -Stressed  stressed  plate  over  non- 
Plate)  urestressed  plate 

10.99  6.9$ 

10.02  14.9* 


*1,2  "*  "1,4 

8.68 

9.97 

14.9* 

"1,2  "3,2 

8.38 

9.61 

14.7* 

Accepted  value 

of 

8.33 

- 

mm 

Ref.  22 

TABLE  37 

Ratio  of  deflection  coefficients  used  in 

Rayleigh-Ritz  analysis 

Ratio  for  Ratio  for 

Coefficients  Hon-Prestresaed  Plate  Pre-Streased  Plate 


"l.A.4 

-19.1 

37.7 

"l,2^W3,2 

12.3 

29.2 
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SBCTIOtf  XV 

Summary  and  Conclusions 

The  proposed  method  of  prestressing  has  been  investigated  and 
shown  to  be  capable  of  raising  the  buckling  load.  It  has  been  shown  that 
the  buckling  load  is  quite  sensitive  to  changes  in  the  magnituie  and 
distribution  of  the  stresses  produced  by  the  prestressing  process* 

The  experimental  program  was  divided  into  two  parts.  In  the 
first  the  deflections  and  stresses  induced  in  the  prestressing  process 
were  determined,  and  the  effect  these  stresses  would  have  upon  the  buckling 
load  was  calculated.  The  calculations  were  performed  by  the  Finite 
Difference  Method  and  by  the  Rayleigh-Ritz  Method.  For  one  plate,  which 
was  the  only  plate  for  which  the  Finite  Difference  Method  was  used,  this 
method  predicted  that  for  the  stresses  found  the  Buckling  Load  would  be 
raised  over  100$  while  the  Rayleigh-Ritz  method  predicted  a raise  of 
approximately. 130$.  Because  it  was  felt  that  the  Finite  Difference  Method 
was  unreliable  for  the  small  number  of  net  points  used  and  also  because 
of  the  lengthy  computations  it  requires,  the  Rayleigh-Ritz  Method  was  decided 
upon  for  determining  the  buckling  loads  of  the  plates.  The  Rayleigh-Ritz 
method  showed  an  increase  of  about  20$  and  13$  using  the  data  of  two  of 
the  other  tests.  It  is  interesting  to  note  that  for  the  latter  two  pre- 
stressed plates  a symmetric  three  half-wave  buckling  mode  was  predicted 
and  for  the  non-pre stressed  plates  of  this  (3:2)  aspect  ratio  an  anti- 
symmetric two  half-wave  buckling  mode  is  predicted. 

Thus  it  was  shown  that  relatively  small  in-plane  stresses  (less 
than  1000  psi)  would  have  a considerable  effect  on  both  the  buckling  mode 
and  load. 


The  second  (and  main)  phase  of  the  experimental  program  was 
that  of  actually  testing  prestressed  plates  to  determine  their  buckling 
loads.  All  the  plates  tested  were  of  the  same  length  and  width  but  were 
of  various  thicknesses  and  had  various  initial  radii  of  curvature.  A 
number  of  non-pre stressed  plates  were  also  tested  for  purposes  of  comparison 
and  calibration  of  the  buckling  jig. 

The  average  increase  of  the  buckling  load  by  prestressing  was 
38$  measured  by  the  "Top  of  the  Knee"  method,  although  there  was  considerable 
scatter.  Some  plates  exhibited  raises  in  their  buckling  load  of  over 
100$  while  other  plates  had  their  buckling  loads  slightly  lowered. 

A need  for  a more  realistic  evaluation  of  the  strength  of  plates 
has  been  pointed  out  since  the  present  concept  of  the  buckling  load  is 
based  upon  the  behavior  of  the  ideal  plate.  In  actuality  there  is  no 
sudden  decrease  in  the  strength  but  only  a gradual  decrease  necessitating 
some  criteria  probably  based  on  the  magnitude  of  the  deflection  as  a ratio 
to  the  applied  load. 
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The  theoretical  program  also  was  divided  into  two  parts.  The 
first  part  consisted  of  finding  the  stresses  and  deflections  induced  by 
the  prestressing  procedure,  while  the  second  part  consisted  of  finding 
their  effect  upon  the  Buckling  Load  by  the  Rayleigh  Ritz  Method. 

Von  Kerman’ s non-linear  equations  were  used  in  finding  the 
deflections  and  stresses.  These  equations  normally  are  used  only  for 
deflections  of  the  order  of  the  plate  thickness,  however,  since  in  this 
case  the  major  deflection  is  from  one  developable  (cylindrical)  surface 
to  another  (flat)  developable  as  to  take  advantage  of  this.  The  methods 
used  in  solving  the  non-linear  equations  is  basically  one  used  by  Levy 
in  which  the  terms  are  expanded  in  Fourier  series  which  transforms  the 
differential  equations  into  algebraic  equations  in  terms  of  the  trigonometric 
coefficients.  The  resulting  non-linear  algebraic  equations  are  then  solved 
by  an  iteration  method. 

The  buckling  load  of  one  particular  is  then  determined  by  the 
Rayleigh  Ritz  method  and  it  is  found  that  in  this  case  the  prestresses  found 
would  raise  the  buckling  load  approximately  15$.  This  conclusion,  however, 
pan  be  subjected  to  questioning,  because  in  the  use  of  Rayleigh-Ritz  method 
some  assumptions  are  made  which  may  not  be  actually  valid. 
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SEC  TIC'!  XVII 
APFEilDIXSS 


APPENDIX  A 

VALIDITY  OF  YOSHIKI'S  METHOD  FOR  DETERMINING  THE  BUCKLING  LOAD  OF  PLATES 

In  the  section  reviewing  pertinent  literature  the  method  developed 
by  M*  Yoshiki  (11)  was  explained  with  reference  to  Figs*  1 and  3*  The 
method  is  based  on  the  parabolic  post  buckling  behavior  of  plates  and 
which  can  be  proven  for  clamped  plates  in  the  following  manner: 

The  total  potential  energy  in  a plate  is  composed  of  the  strain 
energy  due  to  bending,  membrane  stresses  and  the  movement  during  buckling 
of  the  applied  load  according  to 

W " V1  + V2  “ T (Al) 

For  a rectangular  plate  with  clamped  edges  the  strain  energy  due 
to  bending  is 

^ Jf(\72 w -V2  wQ)2  dx  dy  (A 2) 

while  that  due  to  axial  compression  is 


V 


but  since  c » f(x)  only,  and  c*  = E e 


where  Wq  is  the  initial  deflection* 

Denoting  e as  the  total  strain  me as  tired  from  the  unstrained  perfect 
plate,  i*e.  be  is  the  end  shortening  of  the  plate  then 


0 O 


^ * “ 


b a 

//■ 

0 0 


<5y  cy  dx  dy 


(A3) 
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or 


V - • - “ IT-)*  - C-2)  ] *y 

7 2b  y [ 0y  dy  J 


then  eliminating  from  the  expression  for  V2  the  following  expression 
is  arrived  at 


Ebh 


J0\  2bJ0  L87  ^ J J 


dx 


and 


T * ~ p b e 
2 


where  p = - a N 


Assuming  a deflection  during  buckling  of  the  form 


w = w sin  — gin  EZ 

m’n  a b 


sm  MS  sinSIS 


and  initial  imperfections  of  the  form 


. irx  . vy  V* 

w0  “ wmO*nO  Sin  - sin  

aba 


sin 


dw 


(M) 


(A5) 


(A6) 


(A7) 


which  satisfy  the  boundary  conditions  w = 0 and  — “ 0 

0n 

The  energy  expressions  become 

\ « vffi>n2^)3  (2+i*)  2(n2+l)2+8n2+(i)(2+fj)2(m2+l)2 

+8m2+8(-) (m2 +1) (n2+l)^ 

b J 

+wmO,nO  (2+i“°)2(no  +l)2+8n^  +(-)  (2+/°)2(m2+l)2 
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+Sm^+8(S)  (mg+D 

V>VnC»f  ['^1)2+(n-1)2]  £o+<ntl>2£(« 

-(n-1)2  ^-„o+('"1)2^«o3(f[v1>2t(V1)2J Cl 

-("o-D  < V1 ) 2,to-2*  'V1  > 2^moJ 

“m, nV.nO ['f"*1)2<V1)2*(”-1)2(V1)j  4 

-M)2(vi)^-(rt)2(vi)2^(,-i)2(^)2^il3j 

^ ft  I. 

2 n0_<i2+n0“  nO-2  *2-n0} 

'2(^)  Wo  |j(n* l)2(n0+l)2Hn-l)2(n0-l)^g 

-M%0*i)l0-(»*i)V>2^M)2(»o-i)!L 

^ fin  fin  fin  fin  ")  ^ 

2V-a2+mO-4MO-2+J2-moj  "m.nV.nO 
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+ (n0+l)2^  7 

0 24n0j  (A 8) 

which  nay  be  written  in  the  form 

V1  * °lWm,n  +a2W mO,nO+03Wm,nw mO,nO  (A9) 


where  a^,  c^*  and  are  constants  which  are  functions  of  m,n,mO,nO,a, 
and  b only  and  “ 0 if  m * n+l+-p  or  mO  ■ nO+l+-p  where  p is  any  even 
integer* 

For  the  special  case  where  v and  wQ  are  of  the  sane  form,  l*e* , 
m-mQ^r,  n«nQ*p  then 


T1  * °4  (vB,n  - wmO,nO)' 


(Aid) 


where 


Drr4 

4*64a2 


f (S)3(2f  (p4»6p2.l)+(S)  (2*^)  (r4+6r2-l) 


+4(~)  (r2+l)  (p2+l)  (1+^  i{)j 


Also 
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V ■ “P4  -‘-...‘f «’  H)H 

+4n4-(n-l)  2( 3n2+2n+3) ^ 

4+8n2+4n4-(n0-l) 2(  3n2+2nQ+3)  Jg] 

|jl2+8  ^mO+2 

+ ^0-2^^-ro0+  ^2-moj  j4+^n^“^  ( n_1 ) 2^+4n2+4ii2n2 

-2n2(n-l)2^+(n0-l)4^^0-2(l+n2)(n0-l)2^°J 

-»o,»o2  (m) 

which  may  be  written  in  the  fora 

V2  * fiW.»)nA+p2''«D>r,04tP3wm,n2w»0,n02^4o2tp58U».,n2+P6*w-D,n02 

(A12) 
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Again  for  the  special  case  m ■ *q  * r and  n * »q  ■ p 
T2  * ^.*2*  !&(**,»2'**0,nO2)2+  ^*»,i>2"  v»0,n0 


(113) 


where 


y2  - fsb  f^-4jf(P2+D2(p-i)2  f<p-i)2v^} 


The  total  potential  energy  defined  according  to  Eq.  Cl  can  be  found 
for  the  special  case  where  m ■ Bq  ■ r and  n ■ Hq  ■ p to  be 

W ■ a4('W,»0)2*  V*  ^2(w»,n2j,'B0,n02)2+  ^3(V»2'*»0.n02),-^>be 


Nov  minimizing  with  respect  to  e and  v 


(A14) 


m,n 


J ■ 2 * J*; ~ * P b * 0 (H5) 


* “ 2a4(vm,n“*»0,n0>^  ^Wv/^nO2'*2  *3wm,n*  * 0 


m,n 


(A16) 
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Solving  Bqs*  (C15)  and  C16)  fop  e there  results 


• * 7lP  + 2^wm,n^“wmO,nO  * 


(il7) 


and 


al^wm.n*JWm0.n0^  - , 2 2X 

• * *—  * a2<V»  “V.nO  > 


m,n 


Equating  Eqs.  (A17)  and  (A18)  there  results  for  P 


(418) 


- w, 


P * Mwm,n  “^mC^nO  ^ + h 


nO.nO 


a,n 


-h 


now  setting  wm  * & there  results 


2w_ 


2-«„r,  J • ftl*  -22l?2> 


m,n  mO,nO 


(419) 


(420) 


and 


p- L £V ^f22 ) ♦ p, ^2*=2  - $ 


• w 


% 


(A2L) 


However  in  the  post-buckling  region  the  deflection  w becomes  large. 


relative  to  the  initial  deflection  Vq,  that  is 

M0,n0  ^ jtf|SO 

wm,n 


(A22) 


and 


Pl  f2  - P2 


(423) 
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Since  the  strain  gage  readings  taken  during  the  buckling  tests 
were  assumed  proportional  to  the  deflection  caused  by  the  applied  load, 
Eq,  (A23)  demonstrates  the  parabolic  post-buckling  behavior  of  clamped 
plates  assuming  the  plate  to  buckle  in  a mode  similar  to  the  initial 
deflections* 

For  the  prestressed  plate  of  aspect  ratio  3/2  which  is  under 
consideration  the  deflection  forms  which  will  predominate  are 


w ■ w,  9 sin2  — sin  22  sin  ^ 

1,2  a b b 


(A24) 


w^w,  , sin2  — sin2  ^ 
0 10,10  a b 


(A25) 


For  these  forms  of  deflections  the  bending  and  membrane  strain 
energies  can  be  found  to  be 


“St  V1  ' W1.22f246(b) 


+ W1  1 2 fi 

10*10  t 


&(-)  + 48(-) 

b a 


+ 32  (*)J 


-216  w,  t w1  (-) 

10,10  L’z  b 


fcA7* 


8192b4  2 
/ 


(A26) 


2b4  2 4 4 2 2 

j e +87 5w,  a +140W-,  i —68  Ow^  -»  w«  . 

TT4  1,2  10,i0  10»10  1,2 


ev1  22  + 2048  2 

u4  X’2  ^ Wl0,l0 


(A27) 
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The  conditions  that  ~ ■ 0 and  — 


0 will  yield 


75.,  A**. 

1,2  10*10  Eh  w1>2  Eah  Eh  a4  a 

letting  wif2“wl0tl0  “ cT  and  using  Eq*  (^2<^ 

75  ^ 2(1+  ^S^0^5$v  2.  621^  Vjo  _ l2flOFb4  + lgg 

V1*)  Eh  w,  „ Eah  Eh 


JjL23  ♦ 16  ^ ♦ HO  J-  0 


(A28) 


hot  according  to  relations  (A23)  resulting  frc«  the  fact  that  v » 
wrtD  nO  ney  he  written  as  m,B 


75<f2.12SS^  +12S  [^.164  +4o4]-o 

Eah  Eh  *-  a4  a2-* 


(A29) 


or  Just  as  in  Eq*  (i23) 


pmK<f2  ♦ ^2 


(A30) 


Thus  it  is  demonstrated  the  parabolic  nature  of  the  post-buckling 
behavior  of  flat  plates  with  initial  deflections*  Since  the  above  did 
not  depend  on  the  actual  values  of  m,mQ,nQ  the  above  derivation  will 

apply  for  assuned  deflections  of  the  general  foms  of  Eqs*  (a6)  and 
(A7). 
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DERIVATION  OF  EQUATIONS  GOVERNING  STRESS  FUNCTION  COEFFICIENTS 


Upon  substituting  (71)  into  the  expressions^ 

V^-ir4  s L b 4 + 42° Ob23  <=0,32*  (a) 

p»0  q»0  & b a b 

also  by  substitution  of  (70) 

«»2  2 _4  ao  #*  n 

(0J£  ) - iHiH.  tl—  2 2 2 2 m(s-n)  (r-*)n 

5x0y  0x4  0y*  LeTb  m*l  n“l  r»2  s*»2  L 


-m(s-n) 


ij  wm,nWr-m, 


___  __  _A  oo  o®  00  ao 

cos  — cos  222  + 2-s-*  2 2 2 2 

S'"B  a b 4a'c:b  m«l  n=l  r°2  u«  - «» 


m(n-u)  (r-m)n+m(n-u)J 


rtrx  „_a  tgry 
cos  — cos  — 

a b 


7T2  ”,  s,  j*  \ [m)w(m)J  VnVt.i 

4a  b hfI  n=l  t«  s*2  ’ * 


cos  — cos  522  + 2^-5  2 2 2 2 

a b 4a  d hfI  n»l  t*=  -«®t f -op 


® Off  Otf 

2 2 2 2 m(n-u) 


tJTX  wry 

COS  — — COS  


Since  (Bl)  has  indices  running  from  zero  tooo(B2)  vd.ll  be  transformed 
into  a similar  form  making  use  of  the  identity  cos(-a)  * cos(+a)  and 
for  purposes  of  simplification  the  four  summations  of  (B2)  will  be  labelled 
2 A,  2 B,  2 C,  2D,  in  the  same  order  as  they  appear  above. 
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SA“  22  2 ^ 2 »(s-n)  |(r-«)n-ai(s-n)lw  w cos^^cos2^2 

4a  b n*l  n»l  r«=2  s«2  ^ J m»n  r“®*8H1  a h 


2 B * 5 SEES  m(n-u)  f(r-m)n+in(n-u)I 

4a V m=l  n=l  r«2  w=l  L J 


v v 


m,n  r-m,n-u 


cos  ^™*C08  2-2  + 


oo  oo 


■spj  SEES 
b 4a  o hfI  n*l  r=2  tpO 


m(n+u)  I (r-m)n+m(n+u)|  w w cos  — cos  212 

L J m*n  r-a^n+u  , 


<*><*>  oo  oo  r-  — , 

2 (j  *s  - p 2 2 2 2 m(s-n)  1 (ja-t)n+m(s— n)  / w w 

4a  d m=l  n=l  t=l  8=2  L J m»n 


^3  OO  OO  OO 


s-n 


tmc  gOTrjr 


COS'"  CO 


? ? ? , p —j 

+ 9 2 E E 2 E m(s-n)  I (m+t)n+m(s^i)/w  w , 

4a  d m=l  n=l  t*0  s*2  L -J  m,n  DH-t, 


s-n 


£22Lftt|  SEE 


COS”  1 cos 


y4  o*  00  <*,  M p 

ED*  o 2 E E 2 E m(n-u)  I (m-t)n-m(n-ti)  / 

4a  b hfI  n=l  t=l  u*0  L — / 


(B3) 


Vm-M-u  cos  — 


^ oo  oo  oo  OO  p ~7 

+ “ 2 2 SEES  m(n+u)  I (m+t)n-m(n+u)  I 
4a  b ro=l  n=l  t*0  u=l  L J 


Wm,nVt,n+u  C0S  ~ cos 

a b 
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S s 2 *2  m(n+u)  f(m-t)n-4n(n+u)7 
4ab  n=l  t*=l  w*l  *-  *“* 


trrx  tgry 
VnVt,n+ucos—  C0S-7 


^ 2 2 2 2 m(n-u)  £(ia*t)n-m(n-u)j 


+ 4aV 


ibpI  n»l  t=0  u“0 


tux  tpry 

W V . cos  cos  — i 

m,n  m+t,n-u  h 


therefor© 


a 2 2 *2  *2 

(2-H-)  -^~  = SA+ZB  + ZC  + ZD  (B4) 

0x9y  9x  0y 


Equating  coefficients  of  (EL)  and  (B4)  according  to  eq*  (67)  the 
following  results: 


2 _ 2 _ 

2 B.  + 2 0* 
k=l  k~l 


— p— i q-i  r ~f 

vhere  ^1*  Jx  nfx  m(,'n)  VW 


if  p 0,  1 and  q ^ 0,  1 


JL“  0,  if  p ■ 0 or  1,  or  q ■ 0 or  1 
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E.  ■ 
1 


p-1  ' <«  p "1 

2 2 m(n-q)  (p-m)n+m(n-q)  / w„  w 

m=l  n=a+I  L J m»n  3 


p-m,n-q 


if  p ^ 0,  1 and  q / 0 


Bl^  * 0,  if  p * 0 or  1,  or  q ■ 0 


p-1  00  r -j 

B0  * 2 2 m(n+q)  I (p-m)n+m(n+q)  j w v 

* ufI  n=l  L J m,n  ] 


p-mfn+q 


if  p / 0,  1 


B,  ■ 0,  if  p » 0 or  1 


r 


C,  ■ 2 2 m(q-n)  I (m-p)n+m(q-n) 

m=p+l  n=l  L 


) I v w 
J m,n  m-p,q-n 


if  p / 0 and  q / 0,  1 
C1  = 0,  if  p • 0,  or  q * 0 or  1 

00  q-i  r -i 

C - 2 2 m(q-n)  (m+p)n+m(q-n)  w w 

* ml  n=l  L J m»n  ®*1>»q-® 

if  q / 0,  1 
“ 0,  if  q ■ 0 or  1 
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D 


1“ 


oo  0O 

Z Z m(n-q) 

iBFl+p  n=l+q 


jj(nH.p)n-m(n-q)J  wm>n' 


m-p,n-q 


if  p / 0 


D-  « 0,  if  p * 0 


oo  mo  r —j 

Z Z m(n+q)  I (m+p)n-m(n+q'  / w v 

15=1  n**!  L J * ‘ 


m+p,n+q 


if  q^O 


D0  = 0,  if  q - 0 


n 


■ mm  OO  OO  p 

D,  * E Z m(n+q)  I (m-p)n-m(n+q)  / v w 

J m=l+p  n-1  L J n,m 


m-p,n+q 


if  p / 0 and  q + 0 


D,  ■ 0,  if  p *>  0 or  q » 0 


>.  * Z Z m(n-q)  (m+p)n-m(n-q)  I w_  w , 
4 m=l  n=l+q  L J m»n 


ifp/0mrq/0 


D,  * 0 if  p ■ 0 and  q = 0 


Making  the  indicated  transformations  (B5)  can  be  put  in  the  required 
form  as  shown  below: 
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E 


piQ  ./2  b 2 &\2 

F*4  4(p  - + q -) 


2 Bn 
- E 

n=l 


(B6) 


letting  m * k and  n ■ t in 

p-l  q-i 


\ • * J jkt(p-k) (q-t)-k2(q-t)2] ''k>t''J>_k>q-t 


if  p / 0,  1 and  q * 0,  1 


= 0,  if  p = 0 or  1,  or  q * 0 or  1 
letting  m « k and  n * t in  Cg 


B 


2 V [kt(k+p)(q-t)*k2(q-t)2jw]t>1.v(k+p)>(q_t) 


2 k=l  t=l 


if  q jf  0,  1 


B?  = 0,  if  q » 0 or  1 


letting  m ■ k+p  and  n * t in 


q-1 


B- 


W W [(k*p)(t)(k)<q'tH(ktp)2('1't)2JV(k+p).*,(q-t) 


if  q / 0,  1 and  p / 0 


Bj  * 0*  if  q * 0 or  1,  or  p * 0 
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letting  m = k and  n * t in  Bj 
P-1 


B * E 2 [kt(p-k)(t+q)+k^(t+q)^J  w. 

4 k**l  t«=l  L K* 


t (p-k),(t+q) 


if  p ft  0,  1 


* 0,  if  p * 0 or  1 


letting  m * k and  n « q+t  in  B^ 

P-1 
2 

k»l  t»l 


V 5,  I t(^)(p-)^2t2jwt((t+q)V(p.k)>) 


if  q ^ 0 and  p / 0,  1 


Bj  “ 0,  if  q ■ 0,  or  p * 0 or  1 


letting  a * k and  n B t in 


B6  * j*  J,  [«<**P)<W>-kW]  wk)tv(t+p)> 


(t+q) 


if  q / 0 


B,  ■ 0,  if  q « 0 
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letting  m * k and  n ■ t+q  in  Dg 

*i  ■ Jj.  r^)(-p)-v]vkj(wv(k+p)>t 

lfq^Oorp^O 
Ej-O,  ifp-q-0 
letting  m = k+p  and  n * t in 

BS  * Ji  £ [0«P)t*(t+q)-(k*p)2(t*,)2Jw(ktp)(tvki(ttij) 


if  q / 0 and  p / 0 


Bg  * 0,  if  q * 0 or  p » 0 

letting  m * k+p  and  n * t+q  in  D. 

4 

ifp/O 
= 0,  if  p * 0 
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APParon  c 


DERIVATION  OF  EQUATIONS  GOVERNING  RELATION  OF  PRESSURE  AND  DEFLECTION 

COEFFICIENTS 

Upon  substituting  (70)  into  the  expression 


V^v-ir4  5 Z v (SL  + 42s1bSS  sin  S22  (Cl) 

npl  n*l  m,n  b*  a b 


also  by  substitution  of  (70)  and  (71) 

a^  a2w  a2?  a2w  a2F  a2w 

— ? — ? + ““s  “ - 2 * 

a y 8x  ax  ay  axay  3x8y 


2 2 2?  /(inq-np)2b  w sin  sin  LsialHl. 

-1  nxl  p»0  q*0  l m’n  a b 


+(mq+np)2b  w sin  sin  ■CSrfl)lD[ 

P»*  m»n  a b 


♦Ul+nprH)  w sin-^S^S  sin 
p,q  m,n 
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Substituting  the  above  transformations  and  making  use  of  the 
identity  -1  oo  in  (C2)  the  following  equation  results; 

2 f(u)  - 2 f (-u) 


— oo 


a2F  a2w  a2?  a2w  * d2!  aS 

— 5 — 2 + — 2 -2 — 

ay  ax  ax  ay  axay  axa y 


— 5“o  2 2 2 2 |7s-n)m-(r-m)nj  v sin  — sin 

4a  d2  m*l  n»l  r=l  s«l  L J r"*»s“n  m'n  a 


srry 

b 


1 i 1 1 Fn^(M):3 


rrrx 


sin  SEE 
b 


- 2 2 2 2 |(n+u)m+(r-ai)n|  \ „ w sin  22  sin  22 

m=l  n»l  r«l  u=l  “*  r-m,n+u  m»  a v 


2 2 2 2 |(s-n)BH-(m-t)n}2b  . w sin  “2  sin 

nt=l  n*l  t=l  s=l  L J ®~t»8“11  m»n  - v 


••2  2 2 2 J(s-n)m+(iftft)n|  2b _ . w sin  ~2  sin 

m=l  n=l  t=l  s=l  L J m»n  a b 

+ 2 2 2 2 [(n-u)m-(m-t)n)  2b  v _ sin  ~2  sin  22 

ml  n-1  t*l  ir=l  L J m-t,n-u  a b 


oo  oo  oo  oo  r*  ~i  o 

2 2 2 2 (n+u)m-(m+t)nj  o w sin  — sin 

m=l  n®l  t*l  u=l  L J Eft»n“u  e 


HX 

b 
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2 ^ .S.  S_  w sin  sin  HI 

n*l  t»l  w*l  * ' - ^ 


OO  OO  oo  oo  |-  np 

2 2 2 2 I (n+u)m-(B-t)nrh  . 

»=1  n«l  t=l  tFl  L J nwc» 


w_  _ sin  — sin 


iSS  -4«  HLZ 


n+u  m,n 


7 £_ 

~6"b 
4aT>  n*l 


z \ 


(C3) 


Equating  coefficients  of  (Cl)  and  (C3)  according  to  equation  (66) 
and  noting  that 


pT  » 2 2 p sin  sin 

T r»l  s=l  r*s  a to 


odd  integers  only 


■2  22  to4 


Pr,.  ■ “"“S,.*  V 


s *» 


n=l 


r s 


^ » 2 2 [(s-n)n-(r-«)n]  \ 


mfI  n*l 


r-m,s-n  n,n 


r oo 


Ag  * 2 2 J(n-s)HH-(rHa)n|^b 


bf>1  n*s 


r-m,n-e  m,n 


- r oo  p 

■ -2  2 £(n+s)m+(r-n)r 


m*l  n*»-s 


r-m,n+s  ra,n 
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14-  Z Z [(.-.Wk)^^, 

Wtt  1^1 

V -Z  Z [~(a-n)«+(m+r)n] 
iff*— r n*i 

mm  ^ OO  p ^ 

*6  - JL  JL  [(,,-)"-(m-r)”J2b*-r,n-.w»,B 

■FT  n»a 

*8  * £,1 
if,*»  n*s 

S * [(^M**)^2^  W 

**  ”M  (04) 

Equation  (C4)  is  put  into  the  required  fora  by  making  the  indicated 
transformations 

letting  m ■ k aid  n * t in 

*i-  £ £ 

WADC  TR  54-8  l81 


JlLW 


s ,z  k*s?  vm  a &s*c  in.  Iks, 

tj.  : 


<d 


a ® k+5^ 


A* 


.;•  IV  fc-frJP* 


0»  na«»  l~  “|  . 

S 2 ■ (t+s)(k+r)-*t|  Tjl 

*0  fail  L -* 


" k*0  t*l  L 

letting  a » k and  a ■ t+s  in  . 

00  00  r- 

A,  * 2 Z It] 

**"  k=l  1j*0 

letting  in  • k and  n ■ t in  A^ 

w.  Ji  HHVa 


w 


V 


letting  a » k and  n * t+s  in  Ag 


t+S 


letting  n « k and  n « t in 

r od 

*»  ■ Ji  i, 


letting  a * k+r  and  n ■ t in  A. 


**’  So  Si  ^‘X^HVtV.i 


letting  a * k and  n ■ t in  A, 


*9*  SlJi 
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APPE1DXX  D 


SATISPACflOW  OF  STRAIGHT  EDGE  BOUND ART  CONDITIO! 

It  can  be  shown  that  the  assumed  form  of  the  strews  fractions,  F 
and  the  deflection  form  will  automatically  satisfy  the  boundary 
condition  that  the  edges  of  the  plate  must  remain  straight.  This  is 
don©  by  showing  that  the  end  shortening  of  the  plate  is  independent 
of  x and  y and  must  therefore  be  constant  along  each  edge. 

Substituting  (70)  and  (71)  into  (83)  and  then  into  (81)  the 
following  relationship  is  arrived  at 


/ 


)i? 


Z b 


E p*0  q*0 


p»q 


(v^-X  - %j)  COS  ^22  cos 

nab 


2 oo  oo  oo  <so  __ 

- 2-*  Z Z Z Z m m*  v wt  t cos  — 

2aZ  bfI  n»l  m»«l  n«-l  %n  m fn  a 


alszl 


cos  ^22  sin  SZ2:  sin  t cbc 

a b b J 


(Dl) 


It  can  be  seen  from  Eq,  (B6)  that  b^  * 0 and  also  that 

2 2 


m C° ° o 00  * . , . , _ 

bn.  Z ir  (q-t)q  w.  .vv  Z Z kqv.  .v.  .1 

0,q  4q4  a2  (k»l  t»l  k,t  k,q“t  k«l  t=l  k*q+tJ 


(D2) 


Integrating  (81)  and  setting  bQ  q » 0 


<f--“ 


E b 


2 b nq2eosSEl  - 2? 
q*l  0,q  b 4a 


Z Z Z 

d=1  n«l  n'» 


2 

m w 


m,n  m,n* 


sin— — *^sin 
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or 


<r  • - 


— % 2 b q2  cos  SZ2  - — 2 2 2 a2w  w , 

E b2  q»0  0,q  b 8a  dfI  n»l  n*»l  “*n  %n 


[cos  (n-n1)  ^ - cos  (n+n*)  22  j 

L b b J 


0>3) 


The  first  summation  of  (D3)  will  be  examined  separately.  Using 
(D2)  it  m^r  be  written 


oo 


f5 


E b q=l 


asz 


~ ^ „ 


4a  k*l  t»l  q«2 


k,t 


“k  a t cos  ^ ~ 2 2 2 ^k  twk  n+t  008  322 
k*q-t  b 4a  k-1  t»l  q=l  k,t  k,q+t  b 


or 


CO 


-=2  S b / COB  2“ 
E b q=l  °»q  b 


2 00  q— 1 oo  0 

— 2 2 2 k2w,  .w  .cos  322 

4a  k»l  t=l  q*2  k>t  ’q~%  b 


^ co  <1-1  oo 

+ — 2 2 2 x “ Wi.  +v,  . cos 

4a  k*l  t=l  q*2  q * iC»q“t  b 


2EZ 


2 ~ ex?  ^ 

+ — 2 2 2 k v.  tv,  cos  ^ 

4a  k=l  t*l  q*l  k,t  k»q+t  b 


which  may  be  written 
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oo 


TT  a “ v 2 au 

775  E,  b0,qq  cos  “ 

E b q=l  b 


— 2 V 2 k\  tw  tcos  322 

4a  ML  trt  q=2  k,t  k,q“t  b 


JL  ~ q-1  00  2 t 

2-  2 2 2k-  v,.  .w. 


4a  k=l  t=l  q=2  q 


k,t  k,q-t 


cos 


222 


whew  t q-t 


2 «>  q-1  «■»  2 t 

~ 2 2 2 k*  w,  . w.  cos 

4a  k»l  t=l  q=2  q k,t  k,q~t 


222 

b 


where  q-t  t 


£ ~ ” k?  -2 

4a  kcl  t=l  2 


w.  . cos 
k,t 


2tny 


where  q-t  “ t,  i*e«,  q * 2t 


2 op  co  00  0 

+ —222  kw  w.  cos  322 
4a  H M ,4  k.tk„n  b 


Here  note  that  if  we  restrict  the  first  summation  of  the  above  equations 
to  t q-t  then 


JJ  » q-1  » , 2 «.  q-1  «o  _ 

- — 2 2 2 k2w.  . w cos  322 . _ ZL  2 \ 2 k2w,  cos922 

4a  k»l  t=l  q=2  k,t  k»q“t  b 2a  k=l  t*l  q*2  k,t  k,q~t  b 


where  t > q-t 
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V 

4a 


S 

k=l 


COS 


2trry 

b 


where  t « q-t 

then  changing  subscripts  by  the  noted  transformations  in  order  to  clarify 
the  combining  of  terms  the  following  results? 


OO 


IT2  a - 2 

“ 75  \ Vql  cos  „ 

E b q-1  ^ b 


am  m _ 


00  q-1 

S 2 


2a  k=l  t=l  q*2 


k2vk,t''k„-t 008 


am 

b 


where  t > q->t 


oo 


4a  k-1 


„ .2  2 2trry  tt2  !°  ^I1  " . 2 a 

S k w . . cos  + — E Z Ek-w,  w.  cos 

■ »*  b 4a  k-1  s-1  q-2  q k»8  k*“« 


where  s>q-s  let  t - s 


2 «>  q-l  oo 

+ — E E E k2  3rS  w.  w, 
4a  k-1  s=l  q=2  q ic*<*“8  *« 


,s 


cos 


where  s > q-s  let  q-t  * s 


am 

b 


+ 


r2 


4a 


+ 


8a 


oo  -O  vo 

2 2 2 
k*l  t=l  q=l 


go  Oo  p a 

E E kV,  . cos 
k-1  t=l  k,t  b 


,2 

k w,  . w,  . 
k,t  kfq+t 


cos 


am 

b 
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After  noting  that  the  third  and  fourth  summations  of  this  equation  may- 
be combined  as  follows  % 


E \ 2 (k 2 ~ w,  w + k2  ^ wk  w,  ) cos  2^s- 

v k,s  k,q-s  „ k,q-s  k,s' 


k=l  s=l  q=2  q 


q 


where  s >q-s 


c 2 Z Z k'V,  . w,  . cos 
k*l  t=l  q-2  kft  k,q-t 


am 

b 


where  t > q-t 

the  equation  may  be  written 


TT2 


E bf  q2  cos 
E b^  q=l  °>q  b 


2 a.  q-1  <*>  ? 

- — E E E k v,  . v,  . 
4a  k=l  t=l  q=2  k,t  k,q-t 


cos 


where  t 2 q-t 


2 cw 


TT 


oo  oO 


E E E k w,  . ' w. 


/ +%  % k,t  k,q+t 

4a  k=l  t=l  q=l 


cos 


am 


TT2 


* " ,2  2 2trry 

E E k w , . cos  •L 

8a  k=l  t=l  ’ b 


(D4) 


Upon  examination  of  the  second  term  of  eq.  (D3)  and  noting  the  two  cases 
which  may  occur,  i.e»,  n = n*  and  n / n*  the  following  results? 


£ 

8a  m=l 


S.  ^ - cos(nre') 

3=1  n=l  n'=l  ' * b 


ml 

b J 
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if  2 2 n 

2 2 m v + 


8a  n^l  n=l 


m,n 


2 2 2 any 

I If  AAO  ffff.iA 


2 2 mv  _ COS 


8a  idfI  n=l 


m.n 


OO  OO  OO 


— 2 *2  2 a2w  v cos(n-n')  22 

8a  bfI  n=l  n*«l  m»"  b 


n / n' 


oo  oo  oo 


+ — 2 2 2 ®2w  w cos(n+n*) 

8a  dfI  n=l  n'«l  m»n  m'n  b 


n / n* 


Noting  that 

2 oo  oo 

— 22 
8a  bofI  n*=l 


PO  a __  yr  oo  oo 

2 mv  w . cos(n+n»)  22  *=  — 2 2 2 

n'-l  ",n  b 4a  hfI  n-1  n'-l 


m2u  V cos  (n«')  ^ 
m,n  m,n» 


and 


n»>n 


J2  60  Oo  OO  p 

— 2 2 2 mi w w tco8(nHtt»)^ 

8a  »■!  n-1  n'-l  m,n  B»n  b 


— 222 

b »=1  n=l  n*«l 


rn2wm  v_  . cos  (n-n')  2X 
m#n  mfnf  , 


n>n* 
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the  above  may  be  written  as  follows: 

- — 2 2 2 m2v  v t fcos(n-n')  22  _ cos(n+n*)  ^X~L 

8a  n-1  n=l  n'-l  ",n  L b bJ 


L S m2u„_2-—  2 I 2 m2w_  w .cos(n-o')  12 

™ **  mfn  m,n» 


8a  m»l  n=l  m,n  4a  dfI  n=l  n’»l 


n^n* 


_2  oo  oo  oo  o __ 

+ — 2 2 2 m w w . cos(n+n’ ) 22 

4a  nr=l  n=l  n'=l  * * b 


» ? ? 2 2 2n7Ty 

+ ■—  2 2 m w cos  *• 

8a  m=l  n=l  * b 


Or  making  the  indicated  transformations 
2 

- — 2 2 2 m2w  v . cos(n-n* ) 22  _ cos(n+n' ) 22 

8a  dfI  n*l  n1*!  ’ * b b 


where  n>n*  let  n-n*  e q 


TT2  °Z  % 2 2 TT2  °°  ” c*  2 

— 2 2 m w - — 2 2 2 mv  ,w  ,cos 

8a  m=l  n=l  m,n  4a  m=l  q=l  n»=l  m,q+n  m,n 


where  n * > q-n  since  n*  n let  n+n*  * 


_2  OO  Qo  oo  0 

+ — 2 2 2 m w ,w  Bl  cos  222 

4a  m=l  q=l  n»=l  m»q-n  m,n  b 
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~ 2 2 2nwy 

+ — • E S m v cos  ■ ■'* 

8a  idfI  n»l  m,n  b 


(D5) 


Substituting  (D4)  and  (D5)  in  (D3)  the  following  relationship  results 
for  the  end  shortening  of  the  plate  in  the  x-direction. 


similarly 


Z Z 
hfI  n=l 


2 

m w. 


m.n 


Z Z 

m=l  n=l 


„2w 


m,n 


(D6) 


(D7) 
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